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.

Nearly Constant Switching Frequency Hysteresis Current Controller with Fast Online Computation of Boundary for a 2-Level Induction Motor Drive

A nearly constant switching frequency current hysteresis controller for a 2-level inverter fed induction motor drive is proposed in this paper: The salient features of this controller are fast dynamics for the current, inherent protection against overloads and less switching frequency variation. The large variation of switching frequency as in the conventional hysteresis controller is avoided by defining a current-error boundary which is obtained from the current-error trajectory of the standard space vector PWM. The current-error boundary is computed at every sampling interval based on the induction machine parameters and from the estimated fundamental stator voltage. The stator currents are always monitored and when the current-error exceeds the boundary, voltage space vector is switched to reduce the current-error. The proposed boundary computation algorithm is applicable in linear and over-modulation region and it is simple to implement in any standard digital signal processor: Detailed experimental verification is done using a 7.5 kW induction motor and the results are given to show the performance of the drive at various operating conditions and validate the proposed advantages.

Hysteresis current controller for a general n-level inverter fed drive with online current error boundary computation and nearly constant switching frequency

A space vector-based hysteresis current controller for any general n-level three phase inverter fed induction motor drive is proposed in this study. It offers fast dynamics, inherent overload protection and low harmonic distortion for the phase voltages and currents. The controller performs online current error boundary calculations and a nearly constant switching frequency is obtained throughout the linear modulation range. The proposed scheme uses only the adjacent voltage vectors of the present sector, similar to space vector pulse-width modulation and exhibits fast dynamic behaviour under different transient conditions. The steps involved in the boundary calculation include the estimation of phase voltages from the current ripple, computation of switching time and voltage error vectors. Experimental results are given to show the performance of the drive at various speeds, effect of sudden change of the load, acceleration, speed reversal and validate the proposed advantages.

A Nearly Constant Switching Frequency Current Error Space Vector Based Hysteresis Controller for an IM Drive with 12-sided Polygonal Voltage Space Vectors

In this paper, a current error space vector (CESV) based hysteresis controller for a 12-sided polygonal voltage space vector inverter fed induction motor (IM) drive is proposed, for the first time. An open-end winding configuration is used for the induction motor. The proposed controller uses parabolic boundary with generalized vector selection logic for all sectors. The drive scheme is first studied with a space vector based PWM (SVPWM) control and from this the current error space phasor boundary is obtained. This current error space phasor boundary is approximated with four parabolas and then the system is run with space phasor based hysteresis PWM controller by limiting the CESV within the parabolic boundary. The proposed controller has increased modulation range, absence of 5th and 7th order harmonics for the entire modulation range, nearly constant switching frequency, fast dynamic response with smooth transition to the over modulation region and a simple controller implementation.

First-Principles Study of Structure, Vibrational, and Elastic Properties of Stoichiometric and Calcium-Deficient Hydroxyapatite

Hydroxyapatite (HAp), a primary constituent of human bone, is usually nonstoichiometric with varying Ca/P molar ratios, with the well-known fact that Ca deficiency can cause marked reductions in its mechanical properties. To gain insights into the mechanism of this degradation, we employ first-principles calculations based on density functional theory and determine the effects of Ca deficiency on structure, vibrational, and elastic properties of HAp. Our simulation results confirm a considerable reduction in the elastic constants of HAp due to Ca deficiency, which was experimentally reported earlier. Stress-induced transformation of the Ca-deficient defected structure into a metastable state upon the application of stress could be a reason for this. Local structural stability of HAp and Ca-deficient HAp structures is assessed with full phonon dispersion studies. Further, specific signatures in the computed vibrational spectra for Ca deficiency in HAp can be utilized in experimental characterization of different types of defected HAp.

Pyroelectric and second harmonic responses from LiTaO3 nanocrystallites evolved in a Li2O-B2O3-Ta2O5 glass system

Dendritic growth of trigonal and square bipyramidal structures of LiTaO3 nanocrystallites, of 19-30 nm size, was observed when 1.5Li(2)O-2B(2)O(3)-0.5Ta(2)O(5) glasses were subjected to controlled heat treatment between 530 degrees C and 560 degrees C/3 h. X-ray diffraction and Raman spectral studies carried out on the heat-treated samples confirmed the formation of a LiTaO3 phase along with a minor phase of ferroelectric Li2B4O7. The sample that was heat-treated at 550 degrees C/3 h was found to possess similar to 26 nm sized crystallites which exhibited a pyroelectric coefficient as high as 15 nC cm(-2) K-1 which is in the same range (23 nC cm(-2) K-1) as that of single crystalline LiTaO3 at room temperature. The corresponding figures of merit that were calculated for the fast pulse detector (F-i), the large area pyroelectric detector (F-v) and the pyroelectric point detector (F-D) were 0.517 x 10(-10) m V-1, 0.244 m(2) C-1 and 1.437 x 10(-5) Pa-1/2, respectively. Glass nanocrystal composites comprising similar to 30 nm sized crystallites exhibited broad Maker fringes and the second harmonic intensity emanated from these was 0.5 times that of KDP single crystals.

Rapid and facile one-step synthesis of LiTaO3 nanorods

Single crystalline LiTaO3 nanorods (length of 2.5-6 mu m and diameter of 200-500 nm) were synthesized via a facile molten-salt technique. An individual single crystalline nanorod exhibited a piezoelectric coefficient of 8 pm V-1. An improved optical frequency-doubling efficiency was observed in the case of LiTaO3 nanorods as compared to that of cubic crystallites of similar size.

Origin of the Spin-Orbital Liquid State in a Nearly J=0 Iridate Ba3ZnIr2O9

We show using detailed magnetic and thermodynamic studies and theoretical calculations that the ground state of Ba3ZnIr2O9 is a realization of a novel spin-orbital liquid state. Our results reveal that Ba3ZnIr2O9 with Ir5+ (5d(4)) ions and strong spin-orbit coupling (SOC) arrives very close to the elusive J = 0 state but each Ir ion still possesses a weak moment. Ab initio density functional calculations indicate that this moment is developed due to superexchange, mediated by a strong intradimer hopping mechanism. While the Ir spins within the structural Ir2O9 dimer are expected to form a spin-orbit singlet state (SOS) with no resultant moment, substantial frustration arising from interdimer exchange interactions induce quantum fluctuations in these possible SOS states favoring a spin-orbital liquid phase down to at least 100 mK.

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.

Universality in the fast orientational relaxation near isotropic-nematic transition

Detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out to investigate the emergence of criticality in the single-particle orientational relaxation near the isotropic-nematic (IN) phase transition. The simulations show a sudden appearance of a power-law behavior in the decay of the second-rank orientational relaxation as the IN transition is approached. The simulated value of the power-law exponent is 0.56, which is larger than the mean-field value (0.5) but less than the observed value (0.63) and may be due to the finite size of the simulated system. The decay of the first-rank orientational time correlation function, on the other hand, is nearly exponential but its decay becomes very slow near the isotropic-nematic transition, The zero-frequency rotational friction, calculated from the simulated angular Velocity correlation function, shows a marked increase near the IN transition.

Phase conversions in calcium manganites with changing Ca/Mn ratios and their influence on the electrical transport properties

The phase-interconversions between the spinel-, brownmillerite-, defect rocksalt and perovskite-type structures have been investigated by way of (i) introducing deficiency in A-sites in CaxMn2-xO3 (0.05 <= x <= 1) i.e., by varying Ca/Mn ratio from 0.025 to 1 and (ii) nonstoichiometric CaMnO3-delta (CMO) with 0.02 <= delta <= 1. The temperature dependence of resistivity (rho-T) have been investigated on nonstoichiometric CaMnO3-delta (undoped) as well as the CMO substituted with donor impurities such as La3+, Y3+, Bi3+ or acceptor such as Na1+ ion at the Ca-site. The rho-T characteristics of nonstoichiometric CaMnO3-delta is strongly influenced by oxygen deficiency, which controls the concentration of Mn3+ ions and, in turn, affects the resistivity, rho. The results indicated that the substitution of aliovalent impurities at Ca-site in CaMnO3 has similar effects as of CaMnO3-delta ( undoped) annealed in atmospheres of varying partial pressures whereby electron or hole concentration can be altered, yet the doped samples can be processed in air or atmospheres of higher P-O2. The charge transport mechanisms of nonstoichiometric CaMnO3-delta as against the donor or acceptor doped CaMnO3 (sintered in air, P-O2 similar to 0.2 atm) have been predicted. The rho (T) curves of both donor doped CaMnO3 as well as non-stoichiometric CaMnO3-delta, is predictable by the small polaron hopping (SPH) model, which changes to the variable range hopping (VRH) at low temperatures whereas the acceptor doped CaMnO3 exhibited an activated semiconducting hopping ( ASH) throughout the measured range of temperature (10-500 K).

A comparative social network analysis of wasp colonies and classrooms: Linking network structure to functioning

A major question in current network science is how to understand the relationship between structure and functioning of real networks. Here we present a comparative network analysis of 48 wasp and 36 human social networks. We have compared the centralisation and small world character of these interaction networks and have studied how these properties change over time. We compared the interaction networks of (1) two congeneric wasp species (Ropalidia marginata and Ropalidia cyathiformis), (2) the queen-right (with the queen) and queen-less (without the queen) networks of wasps, (3) the four network types obtained by combining (1) and (2) above, and (4) wasp networks with the social networks of children in 36 classrooms. We have found perfect (100%) centralisation in a queen-less wasp colony and nearly perfect centralisation in several other queen-less wasp colonies. Note that the perfectly centralised interaction network is quite unique in the literature of real-world networks. Differences between the interaction networks of the two wasp species are smaller than differences between the networks describing their different colony conditions. Also, the differences between different colony conditions are larger than the differences between wasp and children networks. For example, the structure of queen-right R. marginata colonies is more similar to children social networks than to that of their queen-less colonies. We conclude that network architecture depends more on the functioning of the particular community than on taxonomic differences (either between two wasp species or between wasps and humans).

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

Gel-coated ion-exchange resin: A new kinetic model

Regenerable 'gel-coated' cationic resins with fast sorption kinetics and high sorption capacity have application potential for removal of trace metal ions even in large-scale operations. Poly(acrylic acid) has been gel-coated on high-surface area silica (pre-coated with ethylene-vinyl acetate copolymer providing a thin barrier layer) and insolubilized by crosslinking with a low-molecular-weight diepoxide (epoxy equivalent 180 g) in the presence of benzyl dimethylamine catalyst at 70 degrees C, In experiments performed for Ca2+ sorption from dilute aqueous solutions of Ca(NO,),, the gel-coated acrylic resin is found to have nearly 40% higher sorption capacity than the bead-form commercial methacrylic resin Amberlite IRC-50 and also several limes higher rate of sorption. The sorption on the gel-coated sorbent under vigorous agitation has the characteristics of particle diffusion control with homogeneous (gel) diffusion in resin phase. A new mathematical model is proposed for such sorption on gel-coated ion-exchange resin in finite bath and solved by applying operator-theoretic methods. The analytical solution so obtained shows goad agreement with experimental sorption kinetics at relatively low levels (< 70%) of resin conversion.

Structure and magnetotransport properties of the layered manganites REl.2Srl.8Mn,07(R E = La, Pr, Nd)

We have examined the magnetotransport properties and the structure, by Rietveld refinement of powder X-ray data, of the phases RE(1.2)Sr(1.8)Mn(2)O(7) (RE = La, Pr, Nd). We find that on cooling, La1.2Sr1.8Mn2O7 undergoes a transition to a nearly perfect ferromagnet with 90% magnetization at 1.45 T, as reported by earlier workers, but the Pr and Nd phases show only a small magnetization that grows gradually as the temperature is decreased. There seems to be significant correlation between electrical transport and the Jahn-Teller elongation of the apical Mn-O bonds in these systems. The elongation of the apical Mn-O bonds forces the nine-coordinate rock-salt site to be occupied preferentially by the smaller rare-earth-metal cations. This preferential occupation is reliably obtained from the X-ray refinement. All three title phases show a magnetoresistance ratio of about 4(corresponding to a magnetoresistance, [R(0)-R(H)]/R(0), of about 75%) at a field of 7 T and temperatures around 100 K.