Polarization Caging in Diffusion-Controlled Electron Transfer Reactions in Solution

In some bimolecular diffusion-controlled electron transfer (ET) reactions such as ion recombination (IR), both solvent polarization relaxation and the mutual diffusion of the reacting ion pair may determine the rate and even the yield of the reaction. However, a full treatment with these two reaction coordinates is a challenging task and has been left mostly unsolved. In this work, we address this problem by developing a dynamic theory by combining the ideas from ET reaction literature and barrierless chemical reactions. Two-dimensional coupled Smoluchowski equations are employed to compute the time evolution of joint probability distribution for the reactant (P-(1)(X,R,t)) and the product (p((2))(X,R,t)), where X, as is usual in ET reactions, describes the solvent polarization coordinate and R is the distance between the reacting ion pair. The reaction is described by a reaction line (sink) which is a function of X and R obtained by imposing a condition of equal energy on the initial and final states of a reacting ion pair. The resulting two-dimensional coupled equations of motion have been solved numerically using an alternate direction implicit (ADI) scheme (Peaceman and Rachford, J. Soc. Ind. Appl. Math. 1955, 3, 28). The results reveal interesting interplay between polarization relaxation and translational dynamics. The following new results have been obtained. (i) For solvents with slow longitudinal polarization relaxation, the escape probability decreases drastically as the polarization relaxation time increases. We attribute this to caging by polarization of the surrounding solvent, As expected, for the solvents having fast polarization relaxation, the escape probability is independent of the polarization relaxation time. (ii) In the slow relaxation limit, there is a significant dependence of escape probability and average rate on the initial solvent polarization, again displaying the effects of polarization caging. Escape probability increases, and the average rate decreases on increasing the initial polarization. Again, in the fast polarization relaxation limit, there is no effect of initial polarization on the escape probability and the average rate of IR. (iii) For normal and barrierless regions the dependence of escape probability and the rate of IR on initial polarization is stronger than in the inverted region. (iv) Because of the involvement of dynamics along R coordinate, the asymmetrical parabolic (that is, non-Marcus) energy gap dependence of the rate is observed.

Kinematic analysis and design of articulated manipulators with joint motion constraints

For an articulated manipulator with joint rotation constraints, we show that the maximum workspace is not necessarily obtained for equal link lengths but is also determined by the range and mean positions of the joint motions. We present expressions for sectional area, workspace volume, overlap volume and work area in terms of link ratios, mean positions and ranges of joint motion. We present a numerical procedure to obtain the maximum rectangular area that can be embedded in the workspace of an articulated manipulator with joint motion constraints. We demonstrate the use of analytical expressions and the numerical plots in the kinematic design of an articulated manipulator with joint rotation constraints.

Breakage of a drop of inviscid fluid due to a pressure fluctuation at its surface

In this work, an attempt is made to gain a better understanding of the breakage of low-viscosity drops in turbulent flows by determining the dynamics of deformation of an inviscid drop in response to a pressure variation acting on the drop surface. Known scaling relationships between wavenumbers and frequencies, and between pressure fluctuations and velocity fluctuations in the inertial subrange are used in characterizing the pressure fluctuation. The existence of a maximum stable drop diameter d(max) follows once scaling laws of turbulent flow are used to correlate the magnitude of the disruptive forces with the duration for which they act. Two undetermined dimensionless quantities, both of order unity, appear in the equations of continuity, motion, and the boundary conditions in terms of pressure fluctuations applied on the surface. One is a constant of proportionality relating root-mean-square values of pressure and velocity differences between two points separated by a distance l. The other is a Weber number based on turbulent stresses acting on the drop and the resisting stresses in the drop due to interfacial tension. The former is set equal to 1, and the latter is determined by studying the interaction of a drop of diameter equal to d(max) with a pressure fluctuation of length scale equal to the drop diameter. The model is then used to study the breakage of drops of diameter greater than d(max) and those with densities different from that of the suspending fluid. It is found that, at least during breakage of a drop of diameter greater than d(max) by interaction with a fluctuation of equal length scale, a satellite drop is always formed between two larger drops. When very large drops are broken by smaller-length-scale fluctuations, highly deformed shapes are produced suggesting the possibility of further fragmentation due to instabilities. The model predicts that as the dispersed-phase density increases, d(max) decreases.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Throughput Maximization of Delay-Constrained Traffic in Wireless Energy Harvesting Sensors

A wireless Energy Harvesting Sensor (EHS) needs to send data packets arriving in its queue over a fading channel at maximum possible throughput while ensuring acceptable packet delays. At the same time, it needs to ensure that energy neutrality is satisfied, i.e., the average energy drawn from a battery should equal the amount of energy deposited in it minus the energy lost due to the inefficiency of the battery. In this work, a framework is developed under which a system designer can optimize the performance of the EHS node using power control based on the current channel state information, when the EHS node employs a single modulation and coding scheme and the channel is Rayleigh fading. Optimal system parameters for throughput optimal, delay optimal and delay-constrained throughput optimal policies that ensure energy neutrality are derived. It is seen that a throughput optimal (maximal) policy is packet delay-unbounded and an average delay optimal (minimal) policy achieves negligibly small throughput. Finally, the influence of the harvested energy profile on the performance of the EHS is illustrated through the example of solar energy harvesting.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

Steady compressible flow of cohesionless granular materials through a wedge-shaped bunker

A continuum model based on the critical-state theory of soil mechanics is used to generate stress, density, and velocity profiles, and to compute discharge rates for the flow of granular material in a mass flow bunker. The bin–hopper transition region is idealized as a shock across which all the variables change discontinuously. Comparison with the work of Michalowski (1987) shows that his experimentally determined rupture layer lies between his prediction and that of the present theory. However, it resembles the former more closely. The conventional condition involving a traction-free surface at the hopper exit is abandoned in favour of an exit shock below which the material falls vertically with zero frictional stress. The basic equations, which are not classifiable under any of the standard types, require excessive computational time. This problem is alleviated by the introduction of the Mohr–Coulomb approximation (MCA). The stress, density, and velocity profiles obtained by integration of the MCA converge to asymptotic fields on moving down the hopper. Expressions for these fields are derived by a perturbation method. Computational difficulties are encountered for bunkers with wall angles θw [gt-or-equal, slanted] 15° these are overcome by altering the initial conditions. Predicted discharge rates lie significantly below the measured values of Nguyen et al. (1980), ranging from 38% at θw = 15° to 59% at θw = 32°. The poor prediction appears to be largely due to the exit condition used here. Paradoxically, incompressible discharge rates lie closer to the measured values. An approximate semi-analytical expression for the discharge rate is obtained, which predicts values within 9% of the exact (numerical) ones in the compressible case, and 11% in the incompressible case. The approximate analysis also suggests that inclusion of density variation decreases the discharge rate. This is borne out by the exact (numerical) results – for the parameter values investigated, the compressible discharge rate is about 10% lower than the incompressible value. A preliminary comparison of the predicted density profiles with the measurements of Fickie et al. (1989) shows that the material within the hopper dilates more strongly than predicted. Surprisingly, just below the exit slot, there is good agreement between theory and experiment.

Novel ABO(3) oxides related to perovskite and YAlO3 structure types in the La-B-V-O (B = Ni, Cu) systems

The synthesis, structure and magnetic properties of mixed-metal oxides of ABO(3) composition in the La-B-V-O (B = Ni, Cu) systems are described in the present paper. While the B = Ni oxides adopt GdFeO3-like perovskite structure containing disordered nickel and vanadium at the octahedral B site, La3Cu2VO9 crystallizes in a YAlO3-type structure. A detailed investigation of the superstructure of nominal La3Cu2VO9 by WDS analysis and Rietveld refinement of powder XRD data reveal that the likely composition of the phase is La13Cu9V4O38.5, where the Cu and V atoms are ordered in a root13a(h) (a(h) = hexagonal a parameter of YAlO3-like subcell) superstructure. Magnetic susceptibility data support the proposed superstructure consisting of triangular Cu-3 clusters. At low temperatures, the magnetic moment corresponds to S = 1/2 per Cu-3 cluster, while at high temperatures the behavior is Curie-Weiss like, showing S = 1/2 per copper. The present work reveals the contrasting behavior of La-Cu-V-O and La-Ni-V-O systems: while a unique line-phase related to YAlO3 structure is formed around La3Cu2VO9 Composition in the copper system, a continuous series of perovskite-GdFeO3 solid solutions, LaNi1-xVxO3 for 0 less than or equal to x less than or equal to 1/3 seems to be obtained in the nickel system, where the oxidation state of nickel varies from 3+ to 2+.

Analysis for the Reflected Waves Transmitted on to the Channel by Tall Grounded Objects Using a Special Case

During lightning strike to a tall grounded object (TGO), reflected current waves from TGO are transmitted on to the channel. With regard to these transmitted waves, there seems to be some uncertainties like: 1) will they get reflected at the main wavefront; and 2) if so, what would be their final status. This study makes an attempt to address these issues considering a special case of strike to a TGO involving equal channel core and TGO radii. A macroscopic physical model for the lightning return stroke is adopted for the intended work. Analysis showed that the waves transmitted on to the channel merges with the main wavefront without any sign of reflection. Investigation revealed that: 1) the nonlinear spatio-temporal resistance profile of the channel at the wavefront is mainly responsible for the same; and 2) the distributed source provides additional support. The earlier findings are not limited to the special case of TGO considered. In spite of considering equal TGO and channel core radii, salient features of the model predicted remote electromagnetic fields agree well with the measured data reported in literature.

Comparative statistical aging studies on an oil-pressboard insulation model

Accelerated aging experiments have been conducted on a representative oil-pressboard insulation model to investigate the effect of constant and sequential stresses on the PD behavior using a built-in phase resolved partial discharge analyzer. A cycle of the applied voltage starting from the zero of the positive half cycle was divided into 16 equal phase windows (Φ1 to Φ16) and partial discharge (PD) magnitude distribution in each phase was determined. Based on the experimental results, three stages of aging mechanism were identified. Gumbel's extreme value distribution of the largest element was used to model the first stage of aging process. Second and subsequent stages were modeled using two-parameter Weibull distribution. Spearman's non-parametric rank correlation test statistic and Kolmogrov-Smirnov two sample test were used to relate the aging process of each phase with the corresponding process of the full cycle. To bring out clearly the effect of stress level, its duration and test procedure on the distribution parameters and hence of the aging process, non-parametric ANOVA techniques like Kruskal-Wallis and Fisher's LSD multiple comparison tests were used. Results of the analysis show that two phases (Φ13 and Φ14) near the vicinity of the negative voltage peak were found to contribute significantly to the aging process and their aging mechanism also correlated well with that of the corresponding full cycle mechanism. Attempts have been made to relate these results with the published work of other workers

Analysis for the Reflected Waves Transmitted on to the Channel by Tall Grounded Objects Using a Special Case

During lightning strike to a tall grounded object (TGO), reflected current waves from TGO are transmitted on to the channel. With regard to these transmitted waves, there seems to be some uncertainties like: 1) will they get reflected at the main wavefront; and 2) if so, what would be their final status. This study makes an attempt to address these issues considering a special case of strike to a TGO involving equal channel core and TGO radii. A macroscopic physical model for the lightning return stroke is adopted for the intended work. Analysis showed that the waves transmitted on to the channel merges with the main wavefront without any sign of reflection. Investigation revealed that: 1) the nonlinear spatio-temporal resistance profile of the channel at the wavefront is mainly responsible for the same; and 2) the distributed source provides additional support. The earlier findings are not limited to the special case of TGO considered. In spite of considering equal TGO and channel core radii, salient features of the model predicted remote electromagnetic fields agree well with the measured data reported in literature.

A Simple Charge Model for Symmetric Double-Gate MOSFETs Adapted to Gate-Oxide-Thickness Asymmetry

Surface-potential-based compact charge models for symmetric double-gate metal-oxide-semiconductor field-effect transistors (SDG-MOSFETs) are based on the fundamental assumption of having equal oxide thicknesses for both gates. However, for practical devices, there will always be some amount of asymmetry between the gate oxide thicknesses due to process variations and uncertainties, which can affect device performance significantly. In this paper, we propose a simple surface-potential-based charge model, which is applicable for tied double-gate MOSFETs having same gate work function but could have any difference in gate oxide thickness. The proposed model utilizes the unique so-far-unexplored quasi-linear relationship between the surface potentials along the channel. In this model, the terminal charges could be computed by basic arithmetic operations from the surface potentials and applied biases, and thus, it could be implemented in any circuit simulator very easily and extendable to short-channel devices. We also propose a simple physics-based perturbation technique by which the surface potentials of an asymmetric device could be obtained just by solving the input voltage equation of SDG devices for small asymmetry cases. The proposed model, which shows excellent agreement with numerical and TCAD simulations, is implemented in a professional circuit simulator through the Verilog-A interface and demonstrated for a 101-stage ring oscillator simulation. It is also shown that the proposed model preserves the source/drain symmetry, which is essential for RF circuit design.

Secure Computation in a Bidirectional Relay

Bidirectional relaying, where a relay helps two user nodes to exchange equal length binary messages, has been an active area of recent research. A popular strategy involves a modified Gaussian MAC, where the relay decodes the XOR of the two messages using the naturally-occurring sum of symbols simultaneously transmitted by user nodes. In this work, we consider the Gaussian MAC in bidirectional relaying with an additional secrecy constraint for protection against a honest but curious relay. The constraint is that, while the relay should decode the XOR, it should be fully ignorant of the individual messages of the users. We exploit the symbol addition that occurs in a Gaussian MAC to design explicit strategies that achieve perfect independence between the received symbols and individual transmitted messages. Our results actually hold for a more general scenario where the messages at the two user nodes come from a finite Abelian group G, and the relay must decode the sum within G of the two messages. We provide a lattice coding strategy and study optimal rate versus average power trade-offs for asymptotically large dimensions.

On the generalized degrees of freedom of the K-user symmetric MIMO Gaussian interference channel

This work derives inner and outer bounds on the generalized degrees of freedom (GDOF) of the K-user symmetric MIMO Gaussian interference channel. For the inner bound, an achievable GDOF is derived by employing a combination of treating interference as noise, zero-forcing at the receivers, interference alignment (IA), and extending the Han-Kobayashi (HK) scheme to K users, depending on the number of antennas and the INR/SNR level. An outer bound on the GDOF is derived, using a combination of the notion of cooperation and providing side information to the receivers. Several interesting conclusions are drawn from the bounds. For example, in terms of the achievable GDOF in the weak interference regime, when the number of transmit antennas (M) is equal to the number of receive antennas (N), treating interference as noise performs the same as the HK scheme and is GDOF optimal. For K >; N/M+1, a combination of the HK and IA schemes performs the best among the schemes considered. However, for N/M <; K ≤ N/M+1, the HK scheme is found to be GDOF optimal.

Channel training signal design for reciprocal multiple antenna systems with beamforming

Fast and efficient channel estimation is key to achieving high data rate performance in mobile and vehicular communication systems, where the channel is fast time-varying. To this end, this work proposes and optimizes channel-dependent training schemes for reciprocal Multiple-Input Multiple-Output (MIMO) channels with beamforming (BF) at the transmitter and receiver. First, assuming that Channel State Information (CSI) is available at the receiver, a channel-dependent Reverse Channel Training (RCT) signal is proposed that enables efficient estimation of the BF vector at the transmitter with a minimum training duration of only one symbol. In contrast, conventional orthogonal training requires a minimum training duration equal to the number of receive antennas. A tight approximation to the capacity lower bound on the system is derived, which is used as a performance metric to optimize the parameters of the RCT. Next, assuming that CSI is available at the transmitter, a channel-dependent forward-link training signal is proposed and its power and duration are optimized with respect to an approximate capacity lower bound. Monte Carlo simulations illustrate the significant performance improvement offered by the proposed channel-dependent training schemes over the existing channel-agnostic orthogonal training schemes.