Subdiffusion in hair bundle dynamics: The role of protein conformational fluctuations

The detection of sound signals in vertebrates involves a complex network of different mechano-sensory elements in the inner ear. An especially important element in this network is the hair bundle, an antenna-like array of stereocilia containing gated ion channels that operate under the control of one or more adaptation motors. Deflections of the hair bundle by sound vibrations or thermal fluctuations transiently open the ion channels, allowing the flow of ions through them, and producing an electrical signal in the process, eventually causing the sensation of hearing. Recent high frequency (0.1-10 kHz) measurements by Kozlov et al. Proc. Natl. Acad. Sci. U. S. A. 109, 2896 (2012)] of the power spectrum and the mean square displacement of the thermal fluctuations of the hair bundle suggest that in this regime the dynamics of the hair bundle are subdiffusive. This finding has been explained in terms of the simple Brownian motion of a filament connecting neighboring stereocilia (the tip link), which is modeled as a viscoelastic spring. In the present paper, the diffusive anomalies of the hair bundle are ascribed to tip link fluctuations that evolve by fractional Brownian motion, which originates in fractional Gaussian noise and is characterized by a power law memory. The predictions of this model for the power spectrum of the hair bundle and its mean square displacement are consistent with the experimental data and the known properties of the tip link. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4768902]

Duty cycling and power management with a network of energy harvesting sensors

In this paper, we study duty cycling and power management in a network of energy harvesting sensor (EHS) nodes. We consider a one-hop network, where K EHS nodes send data to a destination over a wireless fading channel. The goal is to find the optimum duty cycling and power scheduling across the nodes that maximizes the average sum data rate, subject to energy neutrality at each node. We adopt a two-stage approach to simplify the problem. In the inner stage, we solve the problem of optimal duty cycling of the nodes, subject to the short-term power constraint set by the outer stage. The outer stage sets the short-term power constraints on the inner stage to maximize the long-term expected sum data rate, subject to long-term energy neutrality at each node. Albeit suboptimal, our solutions turn out to have a surprisingly simple form: the duty cycle allotted to each node by the inner stage is simply the fractional allotted power of that node relative to the total allotted power. The sum power allotted is a clipped version of the sum harvested power across all the nodes. The average sum throughput thus ultimately depends only on the sum harvested power and its statistics. We illustrate the performance improvement offered by the proposed solution compared to other naive schemes via Monte-Carlo simulations.

Interference planning for multicell OFDM downlink

In this paper we consider the downlink of an OFDM cellular system. The objective is to maximise the system utility by means of fractional frequency reuse and interference planning. The problem is a joint scheduling and power allocation problem. Using gradient scheduling scheme, the above problem is transformed to a problem of maximising weighted sum-rate at each time slot. At each slot, an iterative scheduling and power allocation algorithm is employed to address the weighted sum-rate maximisation problem. The power allocation problem in the above algorithm is a nonconvex optimisation problem. We study several algorithms that can tackle this part of the problem. We propose two modifications to the above algorithms to address practical and computational feasibility. Finally, we compare the performance of our algorithm with some existing algorithms based on certain achieved system utility metrics. We show that the practical considerations do not affect the system performance adversely.

Solid-liquid transition in polydisperse Lennard-Jones systems

We study melting of a face-centered crystalline solid consisting of polydisperse Lennard-Jones spheres with Gaussian polydispersity in size. The phase diagram reproduces the existence of a nearly temperature invariant terminal polydispersity (delta(t) similar or equal to 0.11), with no signature of reentrant melting. The absence of reentrant melting can be attributed to the influence of the attractive part of the potential upon melting. We find that at terminal polydispersity the fractional density change approaches zero, which seems to arise from vanishingly small compressibility of the disordered phase. At constant temperature and volume fraction the system undergoes a sharp transition from crystalline solid to the disordered amorphous or fluid state with increasing polydispersity. This has been quantified by second- and third-order rotational invariant bond orientational order, as well as by the average inherent structure energy. The translational order parameter also indicates a similar sharp structural change at delta similar or equal to 0.09 in case of T* = 1.0, phi = 0.58. The free energy calculation further supports the sharp nature of the transition. The third-order rotationally invariant bond order shows that with increasing polydispersity, the local cluster favors a more icosahedral arrangement and the system loses its local crystalline symmetry. Interestingly, the value of structure factor S(k) of the amorphous phase at delta similar or equal to 0.10 (just beyond the solid-liquid transition density at T* = 1) becomes 2.75, which is below the value of 2.85 required for freezing given by the empirical Hansen-Verlet rule of crystallization, well known in the theory of freezing.

On Computing Amplitude, Phase, and Frequency Modulations Using a Vector Interpretation of the Analytic Signal

The amplitude-modulation (AM) and phase-modulation (PM) of an amplitude-modulated frequency-modulated (AM-FM) signal are defined as the modulus and phase angle, respectively, of the analytic signal (AS). The FM is defined as the derivative of the PM. However, this standard definition results in a PM with jump discontinuities in cases when the AM index exceeds unity, resulting in an FM that contains impulses. We propose a new approach to define smooth AM, PM, and FM for the AS, where the PM is computed as the solution to an optimization problem based on a vector interpretation of the AS. Our approach is directly linked to the fractional Hilbert transform (FrHT) and leads to an eigenvalue problem. The resulting PM and AM are shown to be smooth, and in particular, the AM turns out to be bipolar. We show an equivalence of the eigenvalue formulation to the square of the AS, and arrive at a simple method to compute the smooth PM. Some examples on synthesized and real signals are provided to validate the theoretical calculations.

A Nearly Exact Reformulation of the Girsanov Linearization for Stochastically Driven Nonlinear Oscillators

The Girsanov linearization method (GLM), proposed earlier in Saha, N., and Roy, D., 2007, ``The Girsanov Linearisation Method for Stochastically Driven Nonlinear Oscillators,'' J. Appl. Mech., 74, pp. 885-897, is reformulated to arrive at a nearly exact, semianalytical, weak and explicit scheme for nonlinear mechanical oscillators under additive stochastic excitations. At the heart of the reformulated linearization is a temporally localized rejection sampling strategy that, combined with a resampling scheme, enables selecting from and appropriately modifying an ensemble of locally linearized trajectories while weakly applying the Girsanov correction (the Radon-Nikodym derivative) for the linearization errors. The semianalyticity is due to an explicit linearization of the nonlinear drift terms and it plays a crucial role in keeping the Radon-Nikodym derivative ``nearly bounded'' above by the inverse of the linearization time step (which means that only a subset of linearized trajectories with low, yet finite, probability exceeds this bound). Drift linearization is conveniently accomplished via the first few (lower order) terms in the associated stochastic (Ito) Taylor expansion to exclude (multiple) stochastic integrals from the numerical treatment. Similarly, the Radon-Nikodym derivative, which is a strictly positive, exponential (super-) martingale, is converted to a canonical form and evaluated over each time step without directly computing the stochastic integrals appearing in its argument. Through their numeric implementations for a few low-dimensional nonlinear oscillators, the proposed variants of the scheme, presently referred to as the Girsanov corrected linearization method (GCLM), are shown to exhibit remarkably higher numerical accuracy over a much larger range of the time step size than is possible with the local drift-linearization schemes on their own.

Insertion Reactions of Six-Membered Cyclopalladated N,N `,N `'-Triarylguanidine, Pd{kappa(2)(C,N)-C6H3Me-3(NHC(NHAr)(=NAr))-2)(mu-Br)](2) (Ar=2-MeC6H4) with PhCC C-C(O)OR (R = Me and Et): A Gateway to Second Orthopalladation through Novel Rearrangements

The reaction of Pd{kappa(2)(C,N)-C6H3Me-3-(NHC(NHAr)(=NAr))-2}(mu-Br)](2) (Ar = 2-MeC6H4; 1) with 4 equiv of PhC C-C(O)OMe in CH2Cl2 afforded Pd{kappa(2)(C,N)-C(Ph)=C(C(O)OMe)C(Ph)=C(C(O)-OMe)C6H3Me-3(N=C(NH Ar)(2))-2}Br] (Ar = 2-MeC6H4; 2) in 70% yield, and the aforementioned reaction carried out with 10 equiv of PhC C-C(O)OR (R = Me, and Et) afforded an admixture of two regioisomers of Pd{kappa(3)(N,C,O)-O=C(OR)-C5Ph3(C(O)OR)C(C(O)OR)C6H3Me-3(N=C(NHAr)( 2))- 2}Br] (Ar = 2-MeC6H4; R = Me (3a/3b), Et (4a/4b)) in 80 and 87% yields, respectively. In one attempt, the minor regioisomer, 4b, was isolated from the mixture in 6% yield by fractional crystallization. Palladacycles 3a/3b and 4a/4b, upon stirring in CH2Cl2/MeCN (1/1, v/v) mixture at ambient condition for S days, afforded Pd{eta(3)-allyl,(KN)-N-1)-C-5(C(O)OR)(2)Ph3C-(C(O)OR)C6H3Me-3(N=C(NH Ar)(2))(-2)}Br] (Ar = 2-MeC6H4; R = Me (5a/5b), Et (6a/6b)) in 94 and 93% yields, respectively. Palladacycles 3a/3b and 4a/4b, upon reaction with AgOTf in CH2CH2/Me2C(O) (1/1, v/v) mixture at ambient temperature for 15 min, afforded Pd{kappa(3)(N,C,O)-O=C(OR)C5Ph3(C(O)OR)C(C(O)OR)C6H3Me-3(N=C(NHAr)(2 ))-2}(OTf)] (Ar = 2-MeC6H4; R = Me (7a/7b), Et (8a/8b)) in 79 and 77% yields, respectively. Palladacycles 7a/7b and 8a/ 8b, upon reflux in PhC1 separately for 6 h, or palladacycles 5a/5b and 6a/6b, upon treatment with AgOTf in CH2Cl2/Me2C(O) (7/3, v/v) mixture for 15 min, afforded Pd{(eta(2)-Ph)C5Ph2(C(O)OR)kappa(2)(C,N)-C(C(O)OR)C6H3Me-3(N=C(NHAr) (2))-2}(OTf)] (Ar = 2-MeC6H4; R = Me (9a/9h), Et (10a/10b)) in >= 87% yields. Palladacycles 9a/9b, upon stirring in MeCN in the presence of excess NaOAc followed by crystallization of the reaction mixture in the same solvent, afforded Pd{kappa(3)(N,C,C)-(C6H4)C5Ph2(C(O)OMe)(2)C(C(O)OMe)(2)C6H3Me-3(N=C( NHAr)(2))-2}(NCMe)] (Ar = 2-MeC6H4; 11a/11b) in 82% yield. The new palladacycles were characterized by analytical, IR, and NMR (H-1 and C-13) spectroscopic techniques, and the molecular structures of 2, 3a, 4a, 4b, 5a, 6a, 7a, 9a, 10a, and 11a-d(3) were determined by single crystal X-ray diffraction. The frameworks in the aforementioned palladacycles, except that present in 2, are unprecedented. Plausible pathways for the formation of new palladacycles and the influence of the guanidine unit in 1, substituents in alkynes, reaction conditions, and electrophilicity of the bromide and the triflate upon the frameworks of the insertion products have been discussed.

A finite element method for the Saint-Venant torsion and bending problems for prismatic beams

In this work, we present a finite element formulation for the Saint-Venant torsion and bending problems for prismatic beams. The torsion problem formulation is based on the warping function, and can handle multiply-connected regions (including thin-walled structures), compound and anisotropic bars. Similarly, the bending formulation, which is based on linearized elasticity theory, can handle multiply-connected domains including thin-walled sections. The torsional rigidity and shear centers can be found as special cases of these formulations. Numerical results are presented to show the good coarse-mesh accuracy of both the formulations for both the displacement and stress fields. The stiffness matrices and load vectors (which are similar to those for a variable body force in a conventional structural mechanics problem) in both formulations involve only domain integrals, which makes them simple to implement and computationally efficient. (C) 2014 Elsevier Ltd. All rights reserved.

Functional Corannulene: Diverse Structures, Enhanced Charge Transport, and Tunable Optoelectronic Properties

Chemical functionalization of various hydrocarbons, such as coronene, corannulene, and so forth, shows good promise in electronics applications because of their tunable optoelectronic properties. By using quantum chemical calculations, we have investigated the changes in the corannulene buckybowl structure, which greatly affect its electronic and optical properties when functionalized with different electron-withdrawing imide groups. We find that the chemical nature and position of functional groups strongly regulate the stacking geometry, -stacking interactions, and electronic structure. Herein, a range of optoelectronic properties and structure-property relationships of various imide-functionalized corannulenes are explored and rationalized in detail. In terms of carrier mobility, we find that the functionalization strongly affects the reorganization energy of corannulene, while the enhanced stacking improves hopping integrals, favoring the carrier mobility of crystals of pentafluorophenylcorannulene-5-monoimide. The study shows a host of emerging optoelectronic properties and enhancements in the charge-transport characteristics of functionalized corannulene, which may find possible semiconductor and electronics applications.

Binaural Signal Processing Motivated Generalized Analytic Signal Construction and AM-FM Demodulation

Binaural hearing studies show that the auditory system uses the phase-difference information in the auditory stimuli for localization of a sound source. Motivated by this finding, we present a method for demodulation of amplitude-modulated-frequency-modulated (AM-FM) signals using a ignal and its arbitrary phase-shifted version. The demodulation is achieved using two allpass filters, whose impulse responses are related through the fractional Hilbert transform (FrHT). The allpass filters are obtained by cosine-modulation of a zero-phase flat-top prototype halfband lowpass filter. The outputs of the filters are combined to construct an analytic signal (AS) from which the AM and FM are estimated. We show that, under certain assumptions on the signal and the filter structures, the AM and FM can be obtained exactly. The AM-FM calculations are based on the quasi-eigenfunction approximation. We then extend the concept to the demodulation of multicomponent signals using uniform and non-uniform cosine-modulated filterbank (FB) structures consisting of flat bandpass filters, including the uniform cosine-modulated, equivalent rectangular bandwidth (ERB), and constant-Q filterbanks. We validate the theoretical calculations by considering application on synthesized AM-FM signals and compare the performance in presence of noise with three other multiband demodulation techniques, namely, the Teager-energy-based approach, the Gabor's AS approach, and the linear transduction filter approach. We also show demodulation results for real signals.

Upper Bounds on the Size of Grain-Correcting Codes

In this paper, we revisit the combinatorial error model of Mazumdar et al. that models errors in high-density magnetic recording caused by lack of knowledge of grain boundaries in the recording medium. We present new upper bounds on the cardinality/rate of binary block codes that correct errors within this model. All our bounds, except for one, are obtained using combinatorial arguments based on hypergraph fractional coverings. The exception is a bound derived via an information-theoretic argument. Our bounds significantly improve upon existing bounds from the prior literature.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Unusual eigenvalue spectrum and relaxation in the Levy-Ornstein-Uhlenbeck process

We consider the rates of relaxation of a particle in a harmonic well, subject to Levy noise characterized by its Levy index mu. Using the propagator for this Levy-Ornstein-Uhlenbeck process (LOUP), we show that the eigenvalue spectrum of the associated Fokker-Planck operator has the form (n + m mu)nu where nu is the force constant characterizing the well, and n, m is an element of N. If mu is irrational, the eigenvalues are all nondegenerate, but rational mu can lead to degeneracy. The maximum degeneracy is shown to be 2. The left eigenfunctions of the fractional Fokker-Planck operator are very simple while the right eigenfunctions may be obtained from the lowest eigenfunction by a combination of two different step-up operators. Further, we find that the acceptable eigenfunctions should have the asymptotic behavior vertical bar x vertical bar(-n1-n2 mu) as vertical bar x vertical bar -> infinity, with n(1) and n(2) being positive integers, though this condition alone is not enough to identify them uniquely. We also assert that the rates of relaxation of LOUP are determined by the eigenvalues of the associated fractional Fokker-Planck operator and do not depend on the initial state if the moments of the initial distribution are all finite. If the initial distribution has fat tails, for which the higher moments diverge, one can have nonspectral relaxation, as pointed out by Toenjes et al. Phys. Rev. Lett. 110, 150602 (2013)].

Generalized Model Predictive Static Programming and Angle-Constrained Guidance of Air-to-Ground Missiles

A new generalized model predictive static programming technique is presented for rapidly solving a class of finite-horizon nonlinear optimal control problems with hard terminal constraints. Two key features for its high computational efficiency include one-time backward integration of a small-dimensional weighting matrix dynamics, followed bya static optimization formulation that requires only a static Lagrange multiplier to update the control history. It turns out that under Euler integration and rectangular approximation of finite integrals it is equivalent to the existing model predictive static programming technique. In addition to the benchmark double integrator problem, usefulness of the proposed technique is demonstrated by solving a three-dimensional angle-constrained guidance problem for an air-to-ground missile, which demands that the missile must meet constraints on both azimuth and elevation angles at the impact point in addition to achieving near-zero miss distance, while minimizing the lateral acceleration demand throughout its flight path. Simulation studies include maneuvering ground targets along with a first-order autopilot lag. Comparison studies with classical augmented proportional navigation guidance and modern general explicit guidance lead to the conclusion that the proposed guidance is superior to both and has a larger capture region as well.

The Tetrablock as a Spectral Set

The tetrablock, roughly speaking, is the set of all linear fractional maps that map the open unit disc to itself. A formal definition of this inhomogeneous domain is given below. This paper considers triples of commuting bounded operators (A,B,P) that have the tetrablock as a spectral set. Such a triple is named a tetrablock contraction. The motivation comes from the success of model theory in another inhomogeneous domain, namely, the symmetrized bidisc F. A pair of commuting bounded operators (S,P) with Gamma as a spectral set is called a Gamma-contraction, and always has a dilation. The two domains are related intricately as the Lemma 3.2 below shows. Given a triple (A, B, P) as above, we associate with it a pair (F-1, F-2), called its fundamental operators. We show that (A,B,P) dilates if the fundamental operators F-1 and F-2 satisfy certain commutativity conditions. Moreover, the dilation space is no bigger than the minimal isometric dilation space of the contraction P. Whether these commutativity conditions are necessary, too, is not known. what we have shown is that if there is a tetrablock isometric dilation on the minimal isometric dilation space of P. then those commutativity conditions necessarily get imposed on the fundamental operators. En route, we decipher the structure of a tetrablock unitary (this is the candidate as the dilation triple) and a tertrablock isometry (the restriction of a tetrablock unitary to a joint invariant sub-space). We derive new results about r-contractions and apply them to tetrablock contractions. The methods applied are motivated by 11]. Although the calculations are lengthy and more complicated, they beautifully reveal that the dilation depends on the mutual relationship of the two fundamental operators, so that certain conditions need to be satisfied. The question of whether all tetrablock contractions dilate or not is unresolved.