Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

Field observations of instantaneous water slopes and horizontal pressure gradients in the swash-zone

Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within +/- 0.001, or 1/20th -1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05-0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly Suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., delta u/delta t > 0 similar to delta p/delta x < 0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.

Online approximations for wind-field models

We study online approximations to Gaussian process models for spatially distributed systems. We apply our method to the prediction of wind fields over the ocean surface from scatterometer data. Our approach combines a sequential update of a Gaussian approximation to the posterior with a sparse representation that allows to treat problems with a large number of observations.

Developments and extensions of Roth's method of double Fourier series with applications to the solution of electromagnetic field problems

The first part of the thesis compares Roth's method with other methods, in particular the method of separation of variables and the finite cosine transform method, for solving certain elliptic partial differential equations arising in practice. In particular we consider the solution of steady state problems associated with insulated conductors in rectangular slots. Roth's method has two main disadvantages namely the slow rate of con­vergence of the double Fourier series and the restrictive form of the allowable boundary conditions. A combined Roth-separation of variables method is derived to remove the restrictions on the form of the boundary conditions and various Chebyshev approximations are used to try to improve the rate of convergence of the series. All the techniques are then applied to the Neumann problem arising from balanced rectangular windings in a transformer window. Roth's method is then extended to deal with problems other than those resulting from static fields. First we consider a rectangular insulated conductor in a rectangular slot when the current is varying sinusoidally with time. An approximate method is also developed and compared with the exact method.The approximation is then used to consider the problem of an insulated conductor in a slot facing an air gap. We also consider the exact method applied to the determination of the eddy-current loss produced in an isolated rectangular conductor by a transverse magnetic field varying sinusoidally with time. The results obtained using Roth's method are critically compared with those obtained by other authors using different methods. The final part of the thesis investigates further the application of Chebyshdev methods to the solution of elliptic partial differential equations; an area where Chebyshev approximations have rarely been used. A poisson equation with a polynomial term is treated first followed by a slot problem in cylindrical geometry.

Approximation of photonic crystal fibres with large air holes by the step index fibre model

An equivalent step index fibre with a silica core and air cladding is used to model photonic crystal fibres with large air holes. We model this fibre for linear polarisation (we focus on the lowest few transverse modes of the electromagnetic field). The equivalent step index radius is obtained by equating the lowest two eigenvalues of the model to those calculated numerically for the photonic crystal fibres. The step index parameters thus obtained can then be used to calculate nonlinear parameters like the nonlinear effective area of a photonic crystal fibre or to model nonlinear few-mode interactions using an existing model.

Analysis of radicals and radical reaction products in cell signalling and biomolecular damage:the long hard road to gold-standard measures

The field of free radical biology and medicine continues to move at a tremendous pace, with a constant flow of ground-breaking discoveries. The following collection of papers in this issue of Biochemical Society Transactions highlights several key areas of topical interest, including the crucial role of validated measurements of radicals and reactive oxygen species in underpinning nearly all research in the field, the important advances being made as a result of the overlap of free radical research with the reinvigorated field of lipidomics (driven in part by innovations in MS-based analysis), the acceleration of new insights into the role of oxidative protein modifications (particularly to cysteine residues) in modulating cell signalling, and the effects of free radicals on the functions of mitochondria, extracellular matrix and the immune system. In the present article, we provide a brief overview of these research areas, but, throughout this discussion, it must be remembered that it is the availability of reliable analytical methodologies that will be a key factor in facilitating continuing developments in this exciting research area.

Demonstration of 10 Gbit/s transmission over 900 km SMF with <400 ns adaptation time using full-field EDC

This reported work significantly extends the reach of 10Gbit/s on-off keying singlemode fibre (SMF) transmission using full-field based electronic dispersion compensation (EDC) to 900 km. In addition, the EDC balances the complexity and the adaptation capability by employing a simple dispersive transmission line with static parameters for coarse dispersion compensation and 16-state maximum likelihood sequence estimation with Gaussian approximation based channel training for adaptive impairment trimming. Improved adaptation times of less than 400 ns for a bit error rate target of 10-3 over distances ranging from 0 to 900 km are reported.

Magnetic resonance imaging of the manipulation of a chemical wave using an inhomogeneous magnetic field

The effects of applied magnetic fields on the traveling wave formed by the reaction of (ethylenediaminetetraacetato)cobalt(II) (Co(II)EDTA2-) and hydrogen peroxide have been studied using magnetic resonance imaging (MRI). It was found that the wave could be manipulated by applying pulsed magnetic field gradients to a sample contained in a vertical cylindrical tube in the 7.0 T magnetic field of the spectrometer. Transverse field gradients decelerated the propagation of the wave down the high-field side of the tube and accelerated it down the low-field side. This control of the wave propagation eventually promoted the formation of a finger on the low-field side of the tube and allowed the wave to be maneuvered within the sample tube. The origin of these effects is rationalized by considering the Maxwell stress arising from the combined homogeneous and inhomogeneous magnetic fields and the magnetic susceptibility gradient across the wave front.

Magnetic resonance imaging of a magnetic field-dependent chemical wave

The magnetic field dependence of the travelling wave formed during the reaction of (ethylenediaminetetraacetato)cobalt (II) (Co(II)EDTA2- ) and hydrogen peroxide was studied using magnetic resonance imaging (MRI). The reaction was investigated in a vertical tube, in which the wave was initiated from above. The wave propagated downwards, initially with a flat wavefront before forming a finger. Magnetic field effects were observed only once the finger had formed. The wave propagation was accelerated by a magnetic field with a negative gradient (i.e., when the field was stronger at the top of the tube than at the bottom) and slightly decelerated by positive field gradients.

Multiscale modelling of heterogeneously catalysed transesterification reaction process:an overview

Biodiesel is fast becoming one of the key transport fuels as the world endeavours to reduce its carbon footprint and find viable alternatives to oil derived fuels. Research in the field is currently focusing on more efficient ways to produce biodiesel, with the most promising avenue of research looking into the use of heterogeneous catalysis. This article presents a framework for kinetic reaction and diffusive transport modelling of the heterogeneously catalysed transesterification of triglycerides into fatty acid methyl esters (FAMEs), unveiled by a model system of tributyrin transesterification in the presence of MgO catalysts. In particular, the paper makes recommendations on multicomponent diffusion calculations such as the diffusion coefficients and molar fluxes from infinite dilution diffusion coefficients using the Wilke and Chang correlation, intrinsic reaction kinetic studies using the Eley-Rideal kinetic mechanism with methanol adsorption as the rate determining steps and multiscale reaction-diffusion process simulation between catalytic porous and bulk reactor scales. © 2013 The Royal Society of Chemistry.

Demonstration of 10 Gbit/s transmission over 900 km SMF with <400 ns adaptation time using full-field EDC

This reported work significantly extends the reach of 10Gbit/s on-off keying singlemode fibre (SMF) transmission using full-field based electronic dispersion compensation (EDC) to 900 km. In addition, the EDC balances the complexity and the adaptation capability by employing a simple dispersive transmission line with static parameters for coarse dispersion compensation and 16-state maximum likelihood sequence estimation with Gaussian approximation based channel training for adaptive impairment trimming. Improved adaptation times of less than 400 ns for a bit error rate target of 10-3 over distances ranging from 0 to 900 km are reported.

Mean field evolution in an open quantum system

In this thesis, we consider N quantum particles coupled to collective thermal quantum environments. The coupling is energy conserving and scaled in the mean field way. There is no direct interaction between the particles, they only interact via the common reservoir. It is well known that an initially disentangled state of the N particles will remain disentangled at times in the limit N -> [infinity]. In this thesis, we evaluate the η-body reduced density matrix (tracing over the reservoirs and the N - η remaining particles). We identify the main disentangled part of the reduced density matrix and obtain the first order correction term in 1/N. We show that this correction term is entangled. We also estimate the speed of convergence of the reduced density matrix as N -> [infinity]. Our model is exactly solvable and it is not based on numerical approximation.

Habitat Map for Lizard Island reef, Australia derived from a photo-transect survey field data collected in December 2011 and September/October 2012

An object based image analysis approach (OBIA) was used to create a habitat map of the Lizard Reef. Briefly, georeferenced dive and snorkel photo-transect surveys were conducted at different locations surrounding Lizard Island, Australia. For the surveys, a snorkeler or diver swam over the bottom at a depth of 1-2m in the lagoon, One Tree Beach and Research Station areas, and 7m depth in Watson's Bay, while taking photos of the benthos at a set height using a standard digital camera and towing a surface float GPS which was logging its track every five seconds. The camera lens provided a 1.0 m x 1.0 m footprint, at 0.5 m height above the benthos. Horizontal distance between photos was estimated by fin kicks, and corresponded to a surface distance of approximately 2.0 - 4.0 m. Approximation of coordinates of each benthic photo was done based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software (www.geospatialexperts.com). Coordinates of each photo were interpolated by finding the gps coordinates that were logged at a set time before and after the photo was captured. Dominant benthic or substrate cover type was assigned to each photo by placing 24 points random over each image using the Coral Point Count excel program (Kohler and Gill, 2006). Each point was then assigned a dominant cover type using a benthic cover type classification scheme containing nine first-level categories - seagrass high (>=70%), seagrass moderate (40-70%), seagrass low (<= 30%), coral, reef matrix, algae, rubble, rock and sand. Benthic cover composition summaries of each photo were generated automatically in CPCe. The resulting benthic cover data for each photo was linked to GPS coordinates, saved as an ArcMap point shapefile, and projected to Universal Transverse Mercator WGS84 Zone 56 South. The OBIA class assignment followed a hierarchical assignment based on membership rules with levels for "reef", "geomorphic zone" and "benthic community" (above).

Near-field diffraction-based focal length determination technique

An accurate and simple technique for determining the focal length of a lens is presented. It consists of measuring the period of the fringes produced by a diffraction grating at the near field when it is illuminated with a beam focused by the unknown lens. In paraxial approximation, the period of the fringes varies linearly with the distance. After some calculations, a simple extrapolation of data is performed to obtain the locations of the principal plane and the focal plane of the lens. Thus, the focal length is obtained as the distance between the two mentioned planes. The accuracy of the method is limited by the collimation degree of the incident beam and by the algorithm used to obtain the period of the fringes. We have checked the technique with two commercial lenses, one convergent and one divergent, with nominal focal lengths (+100±1) mm and (−100±1) mm respectively. We have experimentally obtained the focal lengths resulting into the interval given by the manufacturer but with an uncertainty of 0.1%, one order of magnitude lesser than the uncertainty given by the manufacturer.

Tracing the Chromospheric and Coronal Magnetic Field with AIA, IRIS, IBIS, and ROSA Data

The aim of this study is to explore the suitability of chromospheric images for magnetic modeling of active regions. We use high-resolutionimages (≈0.2"-0.3"), from the Interferometric Bidimensional Spectrometer in the Ca II 8542 Å line, the Rapid Oscillations in the Solar Atmosphere instrument in the Hα 6563Å line, the Interface Region Imaging Spectrograph in the 2796Å line, and compare non-potential magnetic field models obtainedfrom those chromospheric images with those obtained from images of the Atmospheric Imaging Assembly in coronal (171 Å, etc.) and inchromospheric (304 Å) wavelengths. Curvi-linear structures are automatically traced in those images with the OCCULT-2 code, to which we forward-fitted magnetic field lines computed with the Vertical-current Approximation Nonlinear Force Free Field code. We find that the chromospheric images: (1) reveal crisp curvi-linear structures (fibrils, loop segments, spicules) that are extremely well-suited for constraining magnetic modeling; (2) that these curvi-linear structures arefield-aligned with the best-fit solution by a median misalignment angle of μ2 ≈ 4°–7° (3) the free energy computed from coronal data may underestimate that obtained from chromospheric data by a factor of ≈2–4, (4) the height range of chromospheric features is confined to h≲4000 km, while coronal features are detected up to h = 35,000 km; and (5) the plasma-β parameter is β ≈ 10^-5 - 10^-1 for all traced features. We conclude that chromospheric images reveal important magnetic structures that are complementary to coronal images and need to be included in comprehensive magnetic field models, something that is currently not accomodated in standard NLFFF codes.