GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

An efficient hole transport material composite based on poly(3-hexylthiophene) and bamboo-structured carbon nanotubes for high performance perovskite solar cells

In this work, we have developed a new efficient hole transport material (HTM) composite based on poly(3- hexylthiophene) (P3HT) and bamboo-structured carbon nanotubes (BCNs) for CH3NH3PbI3 (MAPbI3) based perovskite solar cells. Compared to pristine P3HT, it is found that the crystallinity of P3HT was significantly improved by addition of BCNs, which led to over one order of magnitude higher conductivity for the composite containing 1–2 wt% BCNs in P3HT. In the meantime, the interfacial charge transfer between the MAPbI3 light absorbing layer and the HTM composite layer based on P3HT/BCNs was two-fold faster than pristine P3HT. More importantly, the HTM film with a superior morphological structure consisting of closely compact large grains was achieved with the composite containing 1 wt% BCNs in P3HT. The study by electrochemical impedance spectroscopy has confirmed that the electron recombination in the solar cells was reduced nearly ten-fold with the addition of 1 wt% carbon nanotubes in the HTM composite. Owing to the superior HTM film morphology and the significantly reduced charge recombination, the energy conversion efficiency of the perovskite solar cells increased from 3.6% for pristine P3HT to 8.3% for P3HT/(1 wt% BCNs) with a significantly enhanced open circuit voltage (Voc) and fill factor (FF). The findings of this work are important for development of new HTM for high performance perovskite solar cells.

Conductivity Scale In Disordered-Systems

A reanalysis of the correction to the Boltzmann conductivity due to maximally crossed graphs for degenerate bands explains why the conductivity scale in many-valley semiconductors is an order of magnitude higher than Mott's "minimum metallic conductivity." With the use of a reasonable assumption for the Boltzmann mean free path, the lowest-order perturbation theory is seen to give a remarkably good, semiquantitative, description of the conductivity variation in both uncompensated doped semiconductors and amorphous alloys.

Good science for improving policy: greenhouse gas emissions from agricultural manures

Australia’s and New Zealand’s major agricultural manure management emission sources are reported to be, in descending order of magnitude: (1) methane (CH4) from dairy farms in both countries; (2) CH4 from pig farms in Australia; and nitrous oxide (N2O) from (3) beef feedlots and (4) poultry sheds in Australia. We used literature to critically review these inventory estimates. Alarmingly for dairy farm CH4 (1), our review revealed assumptions and omissions that when addressed could dramatically increase this emission estimate. The estimate of CH4 from Australian pig farms (2) appears to be accurate, according to industry data and field measurements. The N2O emission estimates for beef feedlots (3) and poultry sheds (4) are based on northern hemisphere default factors whose appropriateness for Australia is questionable and unverified. Therefore, most of Australasia’s key livestock manure management greenhouse gas (GHG) emission profiles are either questionable or are unsubstantiated by region-specific research. Encouragingly, GHG from dairy shed manure are relatively easy to mitigate because they are a point source which can be managed by several ‘close-to-market’ abatement solutions. Reducing these manure emissions therefore constitutes an opportunity for meaningful action sooner compared with the more difficult-to-implement and long-term strategies that currently dominate agricultural GHG mitigation research. At an international level, our review highlights the critical need to carefully reassess GHG emission profiles, particularly if such assessments have not been made since the compilation of original inventories. Failure to act in this regard presents the very real risk of missing the ‘low hanging fruit’ in the rush towards a meaningful response to climate change

Strong gravity and the Yukawa field

A study of strong gravity field coupled to the Yukawa field is carried out for a conformally flat space-time. A quantitative relation between the strong interaction coupling constantg 2/hstrokc and the strong gravity constants (Lambda f~1028 cm–2,G f ~6.6×1030 C.G.S. units) is obtained givingg 2/hstroksim17, which is of the right order of magnitude. This justifies the contention that strong gravity is relevant for elementary particles (e.g., hadrons).

Problems and prospects in neurochemistry

The importance of neurochemistry in understanding the functional basis of the nervous system was emphasized. Attention was drawn to the role of lipids, particularly the sphingolipids,whose metabolic abnormalities lead to 'sphingolipidosis' In the brain and to gangliosides, which show growth-promoting and neuritogenic properties. Several questions that remain to be answered in this area were enumerated. It was pointed out that neurons make a large number of proteins, an order of magnitude higher than other cells, and several of these are yet to be characterized and their functional significance established. Myelination and synapto-genesis are two fundamental processes in brain development. Although much is known about myelin lipids and proteins, it is not known what signals the glial cell receives to initiate myelin synthesis around the axon, In fact, the process of myelination provides an excellent system for studying membrane biogenesis and cell-sell interaction. Great strides were made in the understanding of neurotransmitter receptors and their function in synaptic transmission, but how neurons make synapses with other specific neurons in a preprogrammed manner is not known and requires immediate study. In this context, it was stressed that developmental neurobiology of the human brain could be most profitably done in India. The importance and complexity of signal transduction mechanisms in the brain was explained and many fundamental questions that remain to be answered were discussed. In conclusion, several other areas of contemporary research interest in the nervous system were mentioned and it was suggested that a 'National Committee for Brain Research' be constituted to identify and intensify research programmes in this vital field.

Dielectric relaxation and electrical conductivity in Bi5NbO10 oxygen ion conductors prepared by a modified sol–gel process

Crystalline Bi5NbO10 nanoparticles have been achieved through a modified sol–gel process using a mixture of ethylenediamine and ethanolamine as a solvent. The Bi5NbO10 nanoparticles were characterized by X-ray diffraction (XRD), differential scanning calorimetry/thermogravimetry (DSC/TG), Fourier transform infrared spectroscopy (FT-IR), transmission electron microscopy (TEM) and Raman spectroscopy. The results showed that well-dispersed 5–60 nm Bi5NbO10 nanoparticles were prepared through heat-treating the precursor at 650 °C and the high density pellets were obtained at temperatures lower than those commonly employed. The frequency and temperature dependence of the dielectric constant and the electrical conductivity of the Bi5NbO10 solid solutions were investigated in the 0.1 Hz to 1 MHz frequency range. Two distinct relaxation mechanisms were observed in the plots of dielectric loss and the imaginary part of impedance (Z″) versus frequency in the temperature range of 200–350 °C. The dielectric constant and the loss in the low frequency regime were electrode dependent. The ionic conductivity of Bi5NbO10 solid solutions at 700 °C is 2.86 Ω−1 m−1 which is in same order of magnitude for Y2O3-stabilized ZrO2 ceramics at same temperature. These results suggest that Bi5NbO10 is a promising material for an oxygen ion conductor.

Prospects for wind energy utilisation in Karnataka State

An examination of the data available at 22 meteorological stations in Karnataka State shows that wind velocities in the State as a whole are neither spectacularly high nor negligibly low. The highest winds (annual mean of around 13 km/hr) are experienced in parts of the northern maidan region of the State (Gulbarga, Raichur and Bidar districts) and in Bangalore. The winds are strongly seasonal: typically, the five monsoon months May-September account for about 80% of the annual wind energy flux. Although the data available are inadequate to make precise estimates, they indicate that the total wind energy potential of the State is about an order of magnitude higher than the current electrical energy consumption. The possible exploitation of wind energy for applications in rural areas therefore requires serious consideration, but it is argued that to be successful it is essential to formulate an integrated and carefully planned programme. The output of current windpumps needs to be increased; a doubling should be feasible by the design of suitable load-matching devices. The first cost has to be reduced by careful design, by the use of local materials and skills and by employing a labour-intensive technology. A consideration of the agricultural factors in the northern maidan region of the State shows that there is likely to be a strong need for mechanical assistance in supplemental and life-saving irrigation for the dry crops characteristic of the area. A technological target for a windmill that could find applications in this area would be one with a rotor diameter of about 10 m that can lift about 10,000 litres of water per hour in winds of 10 km/hr (2.8 m/s) hourly average speed and costs less than about Rs 10,000. Although no such windmills exist as of today, the authors believe that achievement of this target is feasible. An examination of various possible scenarios for the use of windmills in this area suggests that with a windpump costing about Rs 12,000, a three hectare farm growing two dry crops a year can expect an annual return of about 150% from an initial investment of about Rs 15,000. It is concluded that it should be highly worthwhile to undertake a coordinated programme for wind energy development that will include more detailed wind surveys in the northern maidan area (as well as some others, such as the Western Ghats), the development of suitable windmill designs and a study of their applications to agriculture as well as to other fields.

An Accurate Model For Interference From Spatially Distributed Shadowed Users in CDMA Uplinks

A detailed characterization of interference power statistics in CDMA systems is of considerable practical and theoretical interest. Such a characterization for uplink inter-cell interference has been difficult because of transmit power control, randomness in the number of interfering mobile stations, and randomness in their locations. We develop a new method to model the uplink inter-cell interference power as a lognormal distribution, and show that it is an order of magnitude more accurate than the conventional Gaussian approximation even when the average number of mobile stations per cell is relatively large and even outperforms the moment-matched lognormal approximation considered in the literature. The proposed method determines the lognormal parameters by matching its moment generating function with a new approximation of the moment generating function for the inter-cell interference. The method is tractable and exploits the elegant spatial Poisson process theory. Using several numerical examples, the accuracy of the proposed method in modeling the probability distribution of inter-cell interference is verified for both small and large values of interference.

Optical Tweezers for Single Molecule Biology

Molecular machinery on the micro-scale, believed to be the fundamental building blocks of life, involve forces of 1-100 pN and movements of nanometers to micrometers. Micromechanical single-molecule experiments seek to understand the physics of nucleic acids, molecular motors, and other biological systems through direct measurement of forces and displacements. Optical tweezers are a popular choice among several complementary techniques for sensitive force-spectroscopy in the field of single molecule biology. The main objective of this thesis was to design and construct an optical tweezers instrument capable of investigating the physics of molecular motors and mechanisms of protein/nucleic-acid interactions on the single-molecule level. A double-trap optical tweezers instrument incorporating acousto-optic trap-steering, two independent detection channels, and a real-time digital controller was built. A numerical simulation and a theoretical study was performed to assess the signal-to-noise ratio in a constant-force molecular motor stepping experiment. Real-time feedback control of optical tweezers was explored in three studies. Position-clamping was implemented and compared to theoretical models using both proportional and predictive control. A force-clamp was implemented and tested with a DNA-tether in presence of the enzyme lambda exonuclease. The results of the study indicate that the presented models describing signal-to-noise ratio in constant-force experiments and feedback control experiments in optical tweezers agree well with experimental data. The effective trap stiffness can be increased by an order of magnitude using the presented position-clamping method. The force-clamp can be used for constant-force experiments, and the results from a proof-of-principle experiment, in which the enzyme lambda exonuclease converts double-stranded DNA to single-stranded DNA, agree with previous research. The main objective of the thesis was thus achieved. The developed instrument and presented results on feedback control serve as a stepping stone for future contributions to the growing field of single molecule biology.

Parallel min-cut placement on reduced hardware SIMD architecture.

Massively parallel SIMD computing is applied to obtain an order of magnitude improvement in the executional speed of an important algorithm in VLSI design automation. The physical design of a VLSI circuit involves logic module placement as a subtask. The paper is concerned with accelerating the well known Min-cut placement technique for logic cell placement. The inherent parallelism of the Min-cut algorithm is identified, and it is shown that a parallel machine based on the efficient execution of the placement procedure.

Photoconductivity in β-AgI

Excitation spectra and transient and steady-state photoconductivity have been studied in undoped and 0.8-mole% Cu-doped single-crystal β-AgI between 150 and 260°K. A single peak in the spectral response was found to occur in each case, at 2.88 eV for undoped and at 2.81 eV for copper-doped specimens at 260 K, the difference being due to a decrease in band gap. The anisotropy due to polarization of incident radiation parallel or perpendicular to the c direction, which is a measure of the energy difference between the Γ9 and Γ7 levels in the valence band, was found to be 0.010 eV. Transient-photoconductivity experiments showed that the hole lifetime was 6 μ sec at 300°K, an order of magnitude larger than the electron lifetime. The hole drift mobility was found to be 12±2 cm2/ V sec at 300°K and limited by traps at a depth of 0.51±0.01 eV with concentration (3-5)×109/cm3 and capture cross section 10-11 cm2. The study of photoconductivity decay versus temperature revealed the presence of shallow hole traps at 0.14±0.02 eV with concentration greater than 1016/cm3 and capture cross section 10-19 cm2. The steady-state photoconductivity was determined by the deep hole traps at 0.51 eV, and showed the presence of shallow electron traps at a depth of 0.28 eV. The trap distribution was found to be substantially the same in pure and copper-doped specimens, indicating the monovalent substitutional role of copper. The effects of iodine annealing, cadmium doping, and heating above the transition temperature were also studied. The possible nature of the traps is discussed.

Impurity diffusion in bismuth single-crystals

Measurements of impurity diffusion of 86Rb, 90Sr, 133Ba, and 137Cs in single crystal Bi were carried out. Diffusion samples were prepared from single crystal Bi by ion implantation. About 1012-1013 ions were implanted, resulting in surface activities approx =104 cpm. After implantation, specimens were annealed for specified times at 220-265 deg C, and tracer penetration profiles were determined by an electrolytic method. A typical penetration profile for 137Cs in Bi showed a linear relationship for log C vs x in with Fick's law for volume diffusion. Laws of grain boundary diffusion were not obeyed and the order of magnitude of the penetration distances was much less than on a grain boundary mechanism. Results were interpreted in terms of a modified Fischer analysis using a kinetic trapping term. Effective half lengths for trapping at a twin boundary were determined for each impurity.

Fatigue in Ti-6Al-4V-B alloys

The addition of small amounts of B to Ti-6Al-4V alloy reduces the as-cast grain size by an order of magnitude and introduces TiB phase into the microstructure. The effects of these microstructural modifications on both the high cycle fatigue and cyclic stress-strain response were investigated. Experimental results show that B addition markedly enhances the fatigue strength of the alloy; however, the influence of prior-beta grain size was found to be only marginal. The presence of TiB particles in the matrix appears to be beneficial with the addition of 0.55 wt.% B to Ti-6Al-4V enhancing the fatigue strength by more than 50%. Strain-controlled fatigue experiments reveal softening in the cyclic stress-strain response, which increases with the B content in the alloy. Transmission electron microscopy of the fatigued specimens indicates that generation of dislocations during cyclic loading and creation of twins due to strain incompatibility between the matrix and the TiB phase are possible reasons for the observed softening. (c) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Is tetramethyleneethane a ground state triplet?

We have examined a number of possible ways by which tetramethyleneethane (TME) can be a ground state triplet, as claimed by experimental studies, in violation of Ovchinnikov’s theorem for alternant hydrocarbons of equal bond lengths. Model exact π calculations of the low-lying states of TME, 3,4-dimethylenefuran and 3,4-dimethylenepyrrole were carried out using a diagrammatic valence bond approach. The calculations failed to yield a triplet ground state even after (a) tuning of electron correlation, (b) breaking alternancy symmetry, and (c) allowing for geometric distortions. In contrast to earlier studies of fine structure constants in other conjugated systems, the computedD andE values of all the low-lying triplet states of TME for various geometries are at least an order of magnitude different from the experimentally reported values. Incorporation of σ-π mixing by means of UHF MNDO calculations is found to favour a singlet ground state even further. A reinterpretation of the experimental results of TME is therefore suggested to resolve the conflict.