Quantitative photoacoustic tomography from boundary pressure measurements: noniterative recovery of optical absorption coefficient from the reconstructed absorbed energy map

We describe a noniterative method for recovering optical absorption coefficient distribution from the absorbed energy map reconstructed using simulated and noisy boundary pressure measurements. The source reconstruction problem is first solved for the absorbed energy map corresponding to single- and multiple-source illuminations from the side of the imaging plane. It is shown that the absorbed energy map and the absorption coefficient distribution, recovered from the single-source illumination with a large variation in photon flux distribution, have signal-to-noise ratios comparable to those of the reconstructed parameters from a more uniform photon density distribution corresponding to multiple-source illuminations. The absorbed energy map is input as absorption coefficient times photon flux in the time-independent diffusion equation (DE) governing photon transport to recover the photon flux in a single step. The recovered photon flux is used to compute the optical absorption coefficient distribution from the absorbed energy map. In the absence of experimental data, we obtain the boundary measurements through Monte Carlo simulations, and we attempt to address the possible limitations of the DE model in the overall reconstruction procedure.

Outward spherical solidification of a superheated melt with time dependent boundary flux

Short-time analytical solutions of solid and liquid temperatures and freezing front have been obtained for the outward radially symmetric spherical solidification of a superheated melt. Although results are presented here only for time dependent boundary flux, the method of solution can be used for other kinds of boundary conditions also. Later, the analytical solution has been compared with the numerical solution obtained with the help of a finite difference numerical scheme in which the grid points change with the freezing front position. An efficient method of execution of the numerical scheme has been discussed in details. Graphs have been drawn for the total solidification times and temperature distributions in the solid.

Direct stability evaluation of power systems with detailed generator models using structure-preserving energy functions

Energy-based direct methods for transient stability analysis are potentially useful both as offline tools for planning purposes as well as for online security assessment. In this paper, a novel structure-preserving energy function (SPEF) is developed using the philosophy of structure-preserving model for the system and detailed generator model including flux decay, transient saliency, automatic voltage regulator (AVR), exciter and damper winding. A simpler and yet general expression for the SPEF is also derived which can simplify the computation of the energy function. The system equations and the energy function are derived using the centre-of-inertia (COI) formulation and the system loads are modelled as arbitrary functions of the respective bus voltages. Application of the proposed SPEF to transient stability evaluation of power systems is illustrated with numerical examples.

The energy balance in the Ramdas layer

On calm clear nights, air at a height of a few decimetres above bare soil can be cooler than the surface by several degrees in what we shall call the Ramdas layer (Ramdas and Atmanathan, 1932). The authors have recently offered a logical explanation for such a lifted temperature minimum, together with a detailed numerical model. In this paper, we provide physical insight into the phenomenon by a detailed discussion of the energy budget in four typical cases, including one with a lifted minimum. It is shown that the net cooling rate near ground is the small difference between two dominant terms, representing respectively radiative upflux from the ground and from the air layers just above ground. The delicate energy balance that leads to the lifted minimum is upset by turbulent transport, by surface emissivity approaching unity, or by high ground cooling rates. The rapid variation of the flux emissivity of humid air is shown to dominate radiative transport near the ground.

Measurement of Gibbs energy of formation of Ca2PbO4 using a solid-state cell with three electrodes

Phase relations in the system Ca-Pb-O at 1100 K have been determined by equilibrating 18 compositions in the ternary and identifying the phases present in quenched samples by X-ray diffraction and energy dispersive X-ray analysis (EDX). Only one ternary compound Ca2PbO4 was found to be present. The compound coexists with CaO and PbO. The intermetallic compounds Ca2Pb, Ca5Pb3 and CaPb and liquid alloys are in equilibrium with CaO. The standard Gibbs energies of formation of Ca2PbO4 (880 - 1100 K) and Pb3O4 (770 - 910 K) were determined using solid-state cells based on yttria-stabilized zirconia as the solid electrolyte. Pure oxygen gas at 0.1 MPa was used as the reference electrode. For measurements on Ca2PbO4, a novel cell design with three electrodes in series, separated by solid electrolyte membranes, was used to avoid polarization of the electrode containing three solid phases. Two three-phase electrodes were used. The first absorbs the electrochemical flux of oxygen from the reference electrode to the measuring electrode. The other three-phase electrode, which is unaffected by the oxygen flux through the solid electrolyte, is used for electromotive force (EMF) measurement. The results from EMF studies were cross-checked using thermogravimetry (TG) under controlled oxygen partial pressures. The stability of Pb3O4 was investigated using a conventional solid-state cell with RuO2 electrodes. The results can be summarized by the following equations: 2CaO + PbO +1/2O(2) --> Ca2PbO4 Delta(r)G degrees/J mol(-1) = (- 128340 + 93.21 T/K) +/- 200 3PbO + 1/2O(2) --> Pb3O4 Delta(r)G degrees/J mol(-1) = (- 70060 + 77.5 T/K) +/- 150

A method of tracking the peak power points for a variable speed wind energy conversion system

In this paper, a method of tracking the peak power in a wind energy conversion system (WECS) is proposed, which is independent of the turbine parameters and air density. The algorithm searches for the peak power by varying the speed in the desired direction. The generator is operated in the speed control mode with the speed reference being dynamically modified in accordance with the magnitude and direction of change of active power. The peak power points in the P-omega curve correspond to dP/domega = 0. This fact is made use of in the optimum point search algorithm. The generator considered is a wound rotor induction machine whose stator is connected directly to the grid and the rotor is fed through back-to-back pulse-width-modulation (PWM) converters. Stator flux-oriented vector control is applied to control the active and reactive current loops independently. The turbine characteristics are generated by a dc motor fed from a commercial dc drive. All of the control loops are executed by a single-chip digital signal processor (DSP) controller TMS320F240. Experimental results show that the performance of the control algorithm compares well with the conventional torque control method.

Coupled electro-mechanical response of an electroactive polymer cantilever structure and its application in energy harvesting

Ionic polymer-metal composites (IPMC), piezoelectric polymer composites and nematic elastomer composites are materials, which exhibit characteristics of both sensors and actuators. Large deformation and curvature are observed in these systems when electric potential is applied. Effects of geometric non-linearity due to the chargeinduced motion in these materials are poorly understood. In this paper, a coupled model for understanding the behavior of an ionic polymer beam undergoing large deformation and large curvature is presented. Maxwell's equations and charge transport equations are considered which couple the distribution of the ion concentration and the pressure gradient along length of a cantilever beam with interdigital electrodes. A nonlinear constitutive model is derived accounting for the visco-elasto-plastic behavior of these polymers and based on the hypothesis that the presence of electrical charge stretches/contracts bonds, which give rise to electrical field dependent softening/hardening. Polymer chain orientation in statistical sense plays a role on such softening or hardening. Elementary beam kinematics with large curvature is considered. A model for understanding the deformation due to electrostatic repulsion between asymmetrical charge distributions across the cross-sections is presented. Experimental evidence that Silver(Ag) nanoparticle coated IPMCs can be used for energy harvesting is reported. An IPMC strip is vibrated in different environments and the electric power against a resistive load is measured. The electrical power generated was observed to vary with the environment with maximum power being generated when the strip is in wet state. IPMC based energy harvesting systems have potential applications in tidal wave energy harvesting, residual environmental energy harvesting to power MEMS and NEMS devices.

Effect of N/Ga flux ratio on transport behavior of Pt/GaN Schottky diodes

GaN films were grown on c-plane sapphire by plasma-assisted molecular beam epitaxy (PAMBE). The effect of N/Ga flux ratio on structural, morphological, and optical properties was studied. The dislocation density found to increase with increasing the N/Ga ratio. The surface morphology of the films as seen by scanning electron microscopy shows pits on the surface and found that the pit density on the surface increases with N/Ga ratio. The room temperature photoluminescence study reveals the shift in band-edge emission toward the lower energy with increase in N/Ga ratio. This is believed to arise from the reduction in compressive stress in the films as is evidenced by room temperature Raman study. The transport studied on the Pt/GaN Schottky diodes showed a significant increase in leakage current with an increase in N/Ga ratio and was found to be caused by the increase in pit density as well as increase in dislocation density in the GaN films. (C) 2011 American Institute of Physics. [doi:10.1063/1.3634116]

Indium flux, growth temperature and RF power induced effects in InN layers grown on GaN/Si substrate by plasma-assisted MBE

In the present work, we report the growth of wurtzite InN epilayers on GaN/Si (1 1 1) substrate by plasma-assisted molecular beam epitaxy (PAMBE). The growth parameters such as indium flux, substrate temperature and RF power affect the crystallographic and morphological properties of InN layers, which were evaluated using high resolution X-ray diffraction (HRXRD) analysis and atomic force microscopy (AFM). It is found that excess indium (In) concentrations and surface roughness were increased with increase in In flux and growth temperature. The intensity of HRXRD (0 0 0 2) peak, corresponding to c-axis orientation has been increased and full width at half maxima (FWHM) has decreased with increase in RF power. It was found that highly c-axis oriented InN epilayers can be grown at 450 degrees C growth temperature, 450 W RF power and 1.30 x 10(-7) mbar In beam equivalent pressure (BEP). The energy gap of InN layers grown by optimizing growth conditions was determined by photoluminescence and optical absorption measurement. (C) 2011 Elsevier B.V. All rights reserved.

Sensitivity of terrestrial water and energy budgets to CO(2)-physiological forcing: an investigation using an offline land model

Increasing concentrations of atmospheric carbon dioxide (CO(2)) influence climate by suppressing canopy transpiration in addition to its well- known greenhouse gas effect. The decrease in plant transpiration is due to changes in plant physiology (reduced opening of plant stomata). Here, we quantify such changes in water flux for various levels of CO(2) concentrations using the National Center for Atmospheric Research's (NCAR) Community Land Model. We find that photosynthesis saturates after 800 ppmv (parts per million, by volume) in this model. However, unlike photosynthesis, canopy transpiration continues to decline at about 5.1% per 100 ppmv increase in CO(2) levels. We also find that the associated reduction in latent heat flux is primarily compensated by increased sensible heat flux. The continued decline in canopy transpiration and subsequent increase in sensible heat flux at elevated CO(2) levels implies that incremental warming associated with the physiological effect of CO(2) will not abate at higher CO(2) concentrations, indicating important consequences for the global water and carbon cycles from anthropogenic CO(2) emissions.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Flux tube model signals for baryon correlations in heavy ion collisions

The flux tube model offers a pictorial description of what happens during the deconfinement phase transition in QCD. The three-point vertices of a flux tube network lead to formation of baryons upon hadronization. Therefore, correlations in the baryon number distribution at the last scattering surface are related to the preceding pattern of the flux tube vertices in the quark-gluon plasma, and provide a signature of the nearby deconfinement phase transition. I discuss the nature of the expected signal, and how to extract it from the experimental data for heavy ion collisions at RHIC and LHC.

Probing the indefinite CP nature of the Higgs boson through decay distributions in the process e(+)e(-) -> t(t)over-bar Phi

The recently discovered scalar resonance at the Large Hadron Collider is now almost confirmed to be a Higgs boson, whose CP properties are yet to be established. At the International Linear Collider with and without polarized beams, it may be possible to probe these properties at high precision. In this work, we study the possibility of probing departures from the pure CP-even case, by using the decay distributions in the process e(+)e(-) -> t (t) over bar Phi, with Phi mainly decaying into a b (b) over bar pair. We have compared the case of a minimal extension of the Standard Model case (model I) with an additional pseudoscalar degree of freedom, with a more realistic case namely the CP-violating two-Higgs doublet model (model II) that permits a more general description of the couplings. We have considered the International Linear Collider with root s = 800 GeV and integrated luminosity of 300 fb(-1). Our main findings are that even in the case of small departures from the CP-even case, the decay distributions are sensitive to the presence of a CP-odd component in model II, while it is difficult to probe these departures in model I unless the pseudoscalar component is very large. Noting that the proposed degrees of beam polarization increase the statistics, the process demonstrates the effective role of beam polarization in studies beyond the Standard Model. Further, our study shows that an indefinite CP Higgs would be a sensitive laboratory to physics beyond the Standard Model.

Studies of grand minima in sunspot cycles by using a flux transport solar dynamo model

We propose that grand minima in solar activity are caused by simultaneous fluctuations in the meridional circulation and the Babcock-Leighton mechanism for the poloidal field generation in the flux transport dynamo model. We present the following results: (a) fluctuations in the meridional circulation are more effective in producing grand minima; (b) both sudden and gradual initiations of grand minima are possible; (c) distributions of durations and waiting times between grand minima seem to be exponential; (d) the coherence time of the meridional circulation has an effect on the number and the average duration of grand minima, with a coherence time of about 30 yr being consistent with observational data. We also study the occurrence of grand maxima and find that the distributions of durations and waiting times between grand maxima are also exponential, like the grand minima. Finally we address the question of whether the Babcock-Leighton mechanism can be operative during grand minima when there are no sunspots. We show that an alpha-effect restricted to the upper portions of the convection zone can pull the dynamo out of the grand minima and can match various observational requirements if the amplitude of this alpha-effect is suitably fine-tuned.

Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse 10(d) x 10(d) matrix, with d > 6. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in nonequilibrium setups without mass transport.