Dynamical symmetry breaking with optimal control: Reducing the number of pieces

We analyze the production of defects during the dynamical crossing of a mean-field phase transition with a real order parameter. When the parameter that brings the system across the critical point changes in time according to a power-law schedule, we recover the predictions dictated by the well-known Kibble-Zurek theory. For a fixed duration of the evolution, we show that the average number of defects can be drastically reduced for a very large but finite system, by optimizing the time dependence of the driving using optimal control techniques. Furthermore, the optimized protocol is robust against small fluctuations.

On the maximum energy of shock-accelerated cosmic rays at ultra-relativistic shocks

The maximum energy to which cosmic rays can be accelerated at weakly magnetised ultra-relativistic shocks is investigated. We demonstrate that for such shocks, in which the scattering of energetic particles is mediated exclusively by ion skin-depth scale structures, as might be expected for a Weibel-mediated shock, there is an intrinsic limit on the maximum energy to which particles can be accelerated. This maximum energy is determined from the requirement that particles must be isotropized in the downstream plasma frame before the mean field transports them far downstream, and falls considerably short of what is required to produce ultra-high-energy cosmic rays. To circumvent this limit, a highly disorganized field is required on larger scales. The growth of cosmic ray-induced instabilities on wavelengths much longer than the ion-plasma skin depth, both upstream and downstream of the shock, is considered. While these instabilities may play an important role in magnetic field amplification at relativistic shocks, on scales comparable to the gyroradius of the most energetic particles, the calculated growth rates have insufficient time to modify the scattering. Since strong modification is a necessary condition for particles in the downstream region to re-cross the shock, in the absence of an alternative scattering mechanism, these results imply that acceleration to higher energies is ruled out. If weakly magnetized ultra-relativistic shocks are disfavoured as high-energy particle accelerators in general, the search for potential sources of ultra-high-energy cosmic rays can be narrowed.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Robust nonadiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of nonadiabatic electronic transitions and test it against exact integration of the time-dependent Schrodinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the nonadiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.

Case study of the uniaxial anisotropic spin-1 bilinear-biquadratic Heisenberg model on a triangular lattice

We study the spin-1 model on a triangular lattice in the presence of a uniaxial anisotropy field using a cluster mean-field (CMF) approach. The interplay among antiferromagnetic exchange, lattice geometry, and anisotropy forces Gutzwiller mean-field approaches to fail in a certain region of the phase diagram. There, the CMF method yields two supersolid phases compatible with those present in the spin-1/2 XXZ model onto which the spin-1 system maps. Between these two supersolid phases, the three-sublattice order is broken and the results of the CMF approach depend heavily on the geometry and size of the cluster. We discuss the possible presence of a spin liquid in this region.

Positron scattering and annihilation on noble-gas atoms

Positron scattering and annihilation on noble-gas atoms is studied ab initio using many-body theory methods for positron energies below the positronium formation threshold. We show that in this energy range, the many-body theory yields accurate numerical results and provides a near-complete understanding of the positron–noble-gas atom system. It accounts for positron-atom and electron-positron correlations, including the polarization of the atom by the positron and the nonperturbative effect of virtual positronium formation. These correlations have a large influence on the scattering dynamics and result in a strong enhancement of the annihilation rates compared to the independent-particle mean-field description. Computed elastic scattering cross sections are found to be in good agreement with recent experimental results and Kohn variational and convergent close-coupling calculations. The calculated values of the annihilation rate parameter Zeff (effective number of electrons participating in annihilation) rise steeply along the sequence of noble-gas atoms due to the increasing strength of the correlation effects, and agree well with experimental data.

Detection of visual-field deterioration by Glaucoma Progression Analysis and Threshold Noiseless Trend programs

BACKGROUND: To compare the ability of Glaucoma Progression Analysis (GPA) and Threshold Noiseless Trend (TNT) programs to detect visual-field deterioration.

METHODS: Patients with open-angle glaucoma followed for a minimum of 2 years and a minimum of seven reliable visual fields were included. Progression was assessed subjectively by four masked glaucoma experts, and compared with GPA and TNT results. Each case was judged to be stable, deteriorated or suspicious of deterioration

RESULTS: A total of 56 eyes of 42 patients were followed with a mean of 7.8 (SD 1.0) tests over an average of 5.5 (1.04) years. Interobserver agreement to detect progression was good (mean kappa = 0.57). Progression was detected in 10-19 eyes by the experts, in six by GPA and in 24 by TNT. Using the consensus expert opinion as the gold standard (four clinicians detected progression), the GPA sensitivity and specificity were 75% and 83%, respectively, while the TNT sensitivity and specificity was 100% and 77%, respectively.

CONCLUSION: TNT showed greater concordance with the experts than GPA in the detection of visual-field deterioration. GPA showed a high specificity but lower sensitivity, mainly detecting cases of high focality and pronounced mean defect slopes.

Simulation of quantum random walks using the interference of a classical field

We suggest a theoretical scheme for the simulation of quantum random walks on a line using beam splitters, phase shifters, and photodetectors. Our model enables us to simulate a quantum random walk using of the wave nature of classical light fields. Furthermore, the proposed setup allows the analysis of the effects of decoherence. The transition from a pure mean-photon-number distribution to a classical one is studied varying the decoherence parameters.

An Investigation of Compressor Map Width Enhancement and a Detailed Analysis of Inducer Flow Field

This paper describes an investigation of various shroud bleed slot configurations of a centrifugal compressor using CFD with a manual multi-block structured grid generation method. The compressor under investigation is used in a turbocharger application for a heavy duty diesel engine of approximately 400 hp. The baseline numerical model has been developed and validated against experimental performance measurements. The influence of the bleed slot flow field on a range of operating conditions between surge and choke has been analysed in detail.
The impact of the returning bleed flow on the incidence at the impeller blade leading edge due to its mixing with the main through-flow has also been studied. From the baseline geometry, a number of modifications to the bleed slot width have been proposed, and a detailed comparison of the flow characteristics performed. The impact of slot variations on the inlet incidence angle has been investigated, highlighting the improvement in surge and choked flow capability. Along with this, the influence of the bleed slot on stabilising the blade passage flow by the suction of the tip and over-tip vortex flow by the slot has been considered near surge.

A Monte Carlo model of DNA double-strand break clustering and rejoining kinetics for the analysis of pulsed-field gel electrophoresis data

In studies of radiation-induced DNA fragmentation and repair, analytical models may provide rapid and easy-to-use methods to test simple hypotheses regarding the breakage and rejoining mechanisms involved. The random breakage model, according to which lesions are distributed uniformly and independently of each other along the DNA, has been the model most used to describe spatial distribution of radiation-induced DNA damage. Recently several mechanistic approaches have been proposed that model clustered damage to DNA. In general, such approaches focus on the study of initial radiation-induced DNA damage and repair, without considering the effects of additional (unwanted and unavoidable) fragmentation that may take place during the experimental procedures. While most approaches, including measurement of total DNA mass below a specified value, allow for the occurrence of background experimental damage by means of simple subtractive procedures, a more detailed analysis of DNA fragmentation necessitates a more accurate treatment. We have developed a new, relatively simple model of DNA breakage and the resulting rejoining kinetics of broken fragments. Initial radiation-induced DNA damage is simulated using a clustered breakage approach, with three free parameters: the number of independently located clusters, each containing several DNA double-strand breaks (DSBs), the average number of DSBs within a cluster (multiplicity of the cluster), and the maximum allowed radius within which DSBs belonging to the same cluster are distributed. Random breakage is simulated as a special case of the DSB clustering procedure. When the model is applied to the analysis of DNA fragmentation as measured with pulsed-field gel electrophoresis (PFGE), the hypothesis that DSBs in proximity rejoin at a different rate from that of sparse isolated breaks can be tested, since the kinetics of rejoining of fragments of varying size may be followed by means of computer simulations. The problem of how to account for background damage from experimental handling is also carefully considered. We have shown that the conventional procedure of subtracting the background damage from the experimental data may lead to erroneous conclusions during the analysis of both initial fragmentation and DSB rejoining. Despite its relative simplicity, the method presented allows both the quantitative and qualitative description of radiation-induced DNA fragmentation and subsequent rejoining of double-stranded DNA fragments. (C) 2004 by Radiation Research Society.

Perimetric progression in open angle glaucoma and the Visual Field Index (VFI)

PURPOSE: To evaluate the changes in the Visual Field Index (VFI) in eyes with perimetric glaucomatous progression, and to compare these against stable glaucoma patients.

PATIENTS AND METHODS: Consecutive patients with open angle glaucoma with a minimum of 6 reliable visual fields and 2 years of follow-up were identified. Perimetric progression was assessed by 4 masked glaucoma experts from different units, and classified into 3 categories: "definite progression," "suspected progression," or "no progression." This was compared with the Glaucoma Progression Analysis (GPA) II and VFI linear regression analysis, where progression was defined as a negative slope with significance of <5%.

RESULTS: Three hundred ninety-seven visual fields from 51 eyes of 39 patients were assessed. The mean number of visual fields was 7.8 (SD 1.1) per eye, and the mean follow-up duration was 63.7 (SD 13.4) months. The mean VFI linear regression slope showed an overall statistically significant difference (P<0.001, analysis of variance) for each category of progression. Using expert consensus opinion as the reference standard, both VFI analysis and GPA II had high specificity (0.93 and 0.90, respectively), but relatively low sensitivity (0.45 and 0.41, respectively).

CONCLUSIONS: The mean VFI regression slope in our cohort of eyes without perimetric progression showed a statistically significant difference compared with those with suspected and definite progression. VFI analysis and GPA II both had similarly high specificity but low sensitivity when compared with expert consensus opinion.

Three dimensional analysis of the lamina cribrosa in glaucoma

Background/aim: Structural changes in the lamina cribrosa have been implicated in the pathogenesis of glaucomatous optic atrophy. The aim of this study was to determine a measure the surface variability of the cup floor in normal subjects and patients with glaucoma. Methods: A sample of age matched normal subjects (NN), patients with low tension glaucoma (LTG), and primary open angle glaucoma (POAG) were included in the study. The glaucoma groups were matched for the severity of the visual field loss. Mean 10 degree topographic images of normal and glaucomatous eyes from the Heidelberg retina tomograph were imported into ERDAS image processing software where topographic analysis of the cup floor could be assessed. Each image was processed using customised spatial filters that calculated the surface depth variation in localised neighbourhood areas across each image. The local change in depth across the cup floor surface was determined and compared between the three clinical groups. Results: The depth variation in the cup floor was largest in normal subjects followed by LTG and POAG. Highly statistically significant differences in surface depth variability of the cup floor existed between normal and LTG (p=0.005), between normal and POAG (p