An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. (C) 2012 Elsevier B.V. All rights reserved.

Influence of Side-Chain on Structural Order and Photophysical Properties in Thiophene Based Diketopyrrolopyrroles: A Systematic Study

In this work, we have synthesized a series of TDPP derivatives with different alkyl groups such as n-hexyl (-C6H13) 3a, 2-ethylhexyl (-(2-C2H5)C6H12) 3b, triethylene glycol mono methyl ether (-(CH2CH2O)(3c)H-3, TEG) 3c, and octadodecyl (-(8-C8H17)C12H22) 3d. N,N dialkylation of Othiophene-diketopyrrolopyrrole (TDPP, 1) strongly influences its solubility, solid state packing, and structural order. These materials allow us to explicitly study the influence of alkyl chain on solid state packing and photophysical properties. TDPP moiety containing two different alkyl groups 3e (TEG and 2-ethylhexyl) and 3f (TEG and n-hexyl) were synthesized for the first time. The absorption spectra of all derivatives exhibited a red shift in solid state when compared to their solution spectra. The type of alkyl chains leads to change in the optical band gaps in solid state. The fluorescence study reveals that TDPP derivatives have strong pi-pi interaction in the solid state and the extent of bathochromic shift is due to combination of intramolecular interaction and formation of aggregates in solid state. This behavior strongly depends on the nature of alkyl chain. The presence of strong C-H center dot center dot center dot O inter chain interactions and CH-pi interactions in solid state exhibits strong influence on the photophysical properties of TDPP chromophore.

Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium

We derive exact expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electromagnetic field. We also evaluate these expressions for the homogeneous case in equilibrium, where all time dependence and external fields vanish. Unlike similar sum rules for the Fermi-Hubbard model, in the Bose-Hubbard model case, the sum rules often depend on expectation values that cannot be determined simply from parameters in the Hamiltonian like the interaction strength and chemical potential but require knowledge of equal-time many-body expectation values from some other source. We show how one can approximately evaluate these expectation values for the Mott-insulating phase in a systematic strong-coupling expansion in powers of the hopping divided by the interaction. We compare the exact moment relations to the calculated moments of spectral functions determined from a variety of different numerical approximations and use them to benchmark their accuracy. DOI: 10.1103/PhysRevA.87.013628

Implications of program phase behavior on timing analysis

Knowledge about program worst case execution time (WCET) is essential in validating real-time systems and helps in effective scheduling. One popular approach used in industry is to measure execution time of program components on the target architecture and combine them using static analysis of the program. Measurements need to be taken in the least intrusive way in order to avoid affecting accuracy of estimated WCET. Several programs exhibit phase behavior, wherein program dynamic execution is observed to be composed of phases. Each phase being distinct from the other, exhibits homogeneous behavior with respect to cycles per instruction (CPI), data cache misses etc. In this paper, we show that phase behavior has important implications on timing analysis. We make use of the homogeneity of a phase to reduce instrumentation overhead at the same time ensuring that accuracy of WCET is not largely affected. We propose a model for estimating WCET using static worst case instruction counts of individual phases and a function of measured average CPI. We describe a WCET analyzer built on this model which targets two different architectures. The WCET analyzer is observed to give safe estimates for most benchmarks considered in this paper. The tightness of the WCET estimates are observed to be improved for most benchmarks compared to Chronos, a well known static WCET analyzer.

A weighted average based external clock synchronization protocol for wireless sensor networks

Clock synchronization is an extremely important requirement of wireless sensor networks(WSNs). There are many application scenarios such as weather monitoring and forecasting etc. where external clock synchronization may be required because WSN itself may consists of components which are not connected to each other. A usual approach for external clock synchronization in WSNs is to synchronize the clock of a reference node with an external source such as UTC, and the remaining nodes synchronize with the reference node using an internal clock synchronization protocol. In order to provide highly accurate time, both the offset and the drift rate of each clock with respect to reference node are estimated from time to time, and these are used for getting correct time from local clock reading. A problem with this approach is that it is difficult to estimate the offset of a clock with respect to the reference node when drift rate of clocks varies over a period of time. In this paper, we first propose a novel internal clock synchronization protocol based on weighted averaging technique, which synchronizes all the clocks of a WSN to a reference node periodically. We call this protocol weighted average based internal clock synchronization(WICS) protocol. Based on this protocol, we then propose our weighted average based external clock synchronization(WECS) protocol. We have analyzed the proposed protocols for maximum synchronization error and shown that it is always upper bounded. Extensive simulation studies of the proposed protocols have been carried out using Castalia simulator. Simulation results validate our theoretical claim that the maximum synchronization error is always upper bounded and also show that the proposed protocols perform better in comparison to other protocols in terms of synchronization accuracy. A prototype implementation of the proposed internal clock synchronization protocol using a few TelosB motes also validates our claim.

Infinite push: a new support vector ranking algorithm that directly optimizes accuracy at the absolute top of the list

Ranking problems have become increasingly important in machine learning and data mining in recent years, with applications ranging from information retrieval and recommender systems to computational biology and drug discovery. In this paper, we describe a new ranking algorithm that directly maximizes the number of relevant objects retrieved at the absolute top of the list. The algorithm is a support vector style algorithm, but due to the different objective, it no longer leads to a quadratic programming problem. Instead, the dual optimization problem involves l1, ∞ constraints; we solve this dual problem using the recent l1, ∞ projection method of Quattoni et al (2009). Our algorithm can be viewed as an l∞-norm extreme of the lp-norm based algorithm of Rudin (2009) (albeit in a support vector setting rather than a boosting setting); thus we refer to the algorithm as the ‘Infinite Push’. Experiments on real-world data sets confirm the algorithm’s focus on accuracy at the absolute top of the list.

A Savitzky-Golay Filtering Perspective of Dynamic Feature Computation

We address the classical problem of delta feature computation, and interpret the operation involved in terms of Savitzky- Golay (SG) filtering. Features such as themel-frequency cepstral coefficients (MFCCs), obtained based on short-time spectra of the speech signal, are commonly used in speech recognition tasks. In order to incorporate the dynamics of speech, auxiliary delta and delta-delta features, which are computed as temporal derivatives of the original features, are used. Typically, the delta features are computed in a smooth fashion using local least-squares (LS) polynomial fitting on each feature vector component trajectory. In the light of the original work of Savitzky and Golay, and a recent article by Schafer in IEEE Signal Processing Magazine, we interpret the dynamic feature vector computation for arbitrary derivative orders as SG filtering with a fixed impulse response. This filtering equivalence brings in significantly lower latency with no loss in accuracy, as validated by results on a TIMIT phoneme recognition task. The SG filters involved in dynamic parameter computation can be viewed as modulation filters, proposed by Hermansky.

Multiuser detection in large-dimension multicode MIMO-CDMA systems with higher-order modulation

In this paper, we are interested in high spectral efficiency multicode CDMA systems with large number of users employing single/multiple transmit antennas and higher-order modulation. In particular, we consider a local neighborhood search based multiuser detection algorithm which offers very good performance and complexity, suited for systems with large number of users employing M-QAM/M-PSK. We apply the algorithm on the chip matched filter output vector. We demonstrate near-single user (SU) performance of the algorithm in CDMA systems with large number of users using 4-QAM/16-QAM/64-QAM/8-PSK on AWGN, frequency-flat, and frequency-selective fading channels. We further show that the algorithm performs very well in multicode multiple-input multiple-output (MIMO) CDMA systems as well, outperforming other linear detectors and interference cancelers reported in the literature for such systems. The per-symbol complexity of the search algorithm is O(K2n2tn2cM), K: number of users, nt: number of transmit antennas at each user, nc: number of spreading codes multiplexed on each transmit antenna, M: modulation alphabet size, making the algorithm attractive for multiuser detection in large-dimension multicode MIMO-CDMA systems with M-QAM.

Electronically tunable zero order resonator based on CRLH-TLs

Composite Right/Left Handed (CRLH) transmission line (TL) based electronically tunable 1.5 cell zero order resonator (ZOR) is demonstrated with microstrip technology by use of varactors. A novel mechanism for DC bias for the varactor is proposed. This is achieved by patterning the ground plane of microstrip thereby reducing the complexity of DC feed mechanism. This approach also mitigates the effect of parasitics arising from DC feed choke appearing in the RF signal path.

Live Soap: Stability, Order, and Fluctuations in Apolar Active Smectics

We construct a hydrodynamic theory of noisy, apolar active smectics in bulk suspension or on a substrate. Unlike purely orientationally ordered active fluids, active apolar smectics can be dynamically stable in Stokesian bulk suspensions. Smectic order in these systems is quasilong ranged in dimension d = 2 and long ranged in d = 3. We predict reentrant Kosterlitz-Thouless melting to an active nematic in our simplest model in d = 2, a nonzero second-sound speed parallel to the layers in bulk suspensions, and that there are no giant number fluctuations in either case. We also briefly discuss possible instabilities in these systems. DOI: 10.1103/PhysRevLett.110.118102

Reconstructing the properties of dark energy using standard sirens

Future space-based gravity wave (GW) experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of GW sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source in only 68% of the cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini et al. T.D. Saini, S.K. Sethi, and V. Sahni, Phys. Rev. D 81, 103009 (2010)] argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as type Ia supernovae, baryon acoustic oscillations, cosmic microwave background, etc. DOI:10.1103/PhysRevD.87.083001

A unified model for the dynamics of driven ribbon with strain and magnetic order parameters

We develop a unified model to explain the dynamics of driven one dimensional ribbon for materials with strain and magnetic order parameters. We show that the model equations in their most general form explain several results on driven magnetostrictive metallic glass ribbons such as the period doubling route to chaos as a function of a dc magnetic field in the presence of a sinusoidal field, the quasiperiodic route to chaos as a function of the sinusoidal field for a fixed dc field, and induced and suppressed chaos in the presence of an additional low amplitude near resonant sinusoidal field. We also investigate the influence of a low amplitude near resonant field on the period doubling route. The model equations also exhibit symmetry restoring crisis with an exponent close to unity. The model can be adopted to explain certain results on magnetoelastic beam and martensitic ribbon under sinusoidal driving conditions. In the latter case, we find interesting dynamics of a periodic one orbit switching between two equivalent wells as a function of an ac magnetic field that eventually makes a direct transition to chaos under resonant driving condition. The model is also applicable to magnetomartensites and materials with two order parameters. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4790845]

Anisotropy induced crossover from weakly to strongly first order melting of two dimensional solids

Melting and freezing transitions in two dimensional (2D) systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems: the melting of 2D molecular solids is still largely unexplored. In order to understand the role of anisotropy as well as multiple energy and length scales present in molecular systems, here we report computer simulation studies of melting of 2D molecular systems. We computed a limited portion of the solid-liquid phase diagram. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram consisting of isotropic liquid and two crystalline phases-honeycomb and oblique. The nature of the transition depends on the relative strength of the anisotropic interaction and a strongly first order melting turns into a weakly first order transition on increasing the strength of the isotropic interaction. This crossover can be attributed to an increase in stiffness of the solid phase free energy minimum on increasing the strength of the anisotropic interaction. The defects involved in melting of molecular systems are quite different from those known for the atomic systems.

Structural characterization and molecular order of rodlike mesogens with three- and four-ring core by XRD and C-13 NMR spectroscopy

Structural characterizations using XRD and C-13 NMR spectroscopy of two rodlike mesogens consisting of (i) three phenyl ring core with a polar cyano terminal and (ii) four phenyl ring core with flexible dodecyl terminal chain are presented. The three-ring-core mesogen with cyano terminal exhibits enantiotropic smectic A phase while the four-ring mesogen reveals polymesomorphism and shows enantiotropic nematic, smectic C, and tilted hexatic phases. The molecular organization in the three-ring mesogen is found to be partial bilayer smectic Ad type, and the interdigitation of the molecules in the neighboring layers is attributed to the presence of the polar terminal group. For the four-ring mesogen, the XRD results confirm the existence of the smectic C and the tilted hexatic mesophases. A thermal variation of the layer spacing across the smectic C phase followed by a discrete jump at the transition to the tilted hexatic phase is also observed. The tilt angles have been estimated to be about 45 degrees in the smectic C phase and about 40 degrees in tilted hexatic phase. C-13 NMR results indicate that in the mesophase the molecules are aligned parallel to the magnetic field. From the C-13-H-1 dipolar couplings determined from the 2D experiments, the overall order parameter for the three-ring mesogen in its smectic A phase has been estimated to be 0.72 while values ranging from 0.88 to 0.44 have been obtained for the four-ring mesogen as it passes from the tilted hexatic to the nematic phase. The orientations of the different rings of the core unit with respect to each other and also with respect to the long axis of the molecule have also been obtained.

The Indian Ocean forecast system

In order to meet the ever growing demand for the prediction of oceanographic parametres in the Indian Ocean for a variety of applications, the Indian National Centre for Ocean Information Services (INCOIS) has recently set-up an operational ocean forecast system, viz. the Indian Ocean Forecast System (INDOFOS). This fully automated system, based on a state-of-the-art ocean general circulation model issues six-hourly forecasts of the sea-surface temperature, surface currents and depths of the mixed layer and the thermocline up to five-days of lead time. A brief account of INDOFOS and a statistical validation of the forecasts of these parametres using in situ and remote sensing data are presented in this article. The accuracy of the sea-surface temperature forecasts by the system is high in the Bay of Bengal and the Arabian Sea, whereas it is moderate in the equatorial Indian Ocean. On the other hand, the accuracy of the depth of the thermocline and the isothermal layers and surface current forecasts are higher near the equatorial region, while it is relatively lower in the Bay of Bengal.