First-order hyperpolarizabilities and pKa of weak organic acids in protic solvents are lineraly related

Second-order nonlinearities (beta) of five weak organic acids in protic solvents have been measured by the double-quantum Rayleigh scattering (DRS) technique. beta is found to bear a linear relationship to the pK(a) of these compounds in those solvents. A direct implication of this observation is that the DRS technique can be used to determine the pK(a) of weak organic acids in any solvent.

First-Order Hyperpolarizabilities of Sulfophthalein Dyes

In this paper we report the first hyperpolarizabilities (beta) of 12, sulfophthalein dyes. Since these dyes are ionic in nature, their second-order nonlinearities were measured by the hyper-Rayleigh scattering technique in solution. The measured beta values are large and highly solvent dependent. Inclusion of solvent polarity in ab initio estimates of static second-order polarizability does not fully account for the experimental beta values. Contributions from the dissociated forms of the dye in different solvents seem to play an important role in enhancing beta in these systems.

Synthesis and esterolytic chemistry of some (dialkylamino) pyridine-functionalized micellar aggregates. Evidence of catalytic turnover

Four new (dialkylamino)pyridine-functionalized surfactants have been synthesized. Micelles were generated either from the surfactant alone in aqueous buffer (pH 8.5 or 9.0) or by comicellization in 1 x 10(-3)-1 x 10(-4) M aqueous micellar cetyltrimethylammonium bromide (CTABr) solution at pH 8.5 or 9.0. Such aggregates were used to cleave p-nitrophenyl alkanoates or p-nitrophenyl diphenylphosphate. The nucleophilic reagents and the second-order ''catalytic'' rate constants toward esterolysis of the substrate p-nitrophenyl octanoate (at 25 degrees C, pH 9.0) were [cat.] = 1 x 10(-4) M, [CTABr] = 1 x 10(-3) M, and k(cat.) = 440.13 M(-1) s(-1) for 1b, [cat.] = 5 x 10(-4) M, [CTABr] = 5 x 10(-4) M, and k(cat.) = 30.8 M(-1) s(-1) for 1c, [cat.] = 5 x 10(-4) M, [CTABr] = 5 x 10(-3) M, and k(cat.) = 183.64 M(-1) s(-1) for 2a, and [cat.] = 3 x 10(-4) M and k(cat.) = 54.1 M(-1) s(-1) for 2b. The catalytic systems, especially 1b/CTABr and 2a/CTABr, also conferred significantly greater reactivity toward the esters derived from alkanoic acids of moderate chain length (C-6-C-10) during hydrolytic cleavages relative to their shorter and longer counterparts. Importantly, the catalytic systems comprising the coaggregates of either neutral 1b and CTABr (1:10) or anionic 2a and CTABr (1:10) conformed to the Michaelis-Menten kinetic scheme and demonstrated turnover behavior in the presence of excess substrate.

Robust Aeroelastic Optimization Of Composite Helicopter Rotor

A robust aeroelastic optimization is performed to minimize helicopter vibration with uncertainties in the design variables. Polynomial response surfaces and space-¯lling experimental designs are used to generate the surrogate model of aeroelastic analysis code. Aeroelastic simulations are performed at the sample inputs generated by Latin hypercube sampling. The response values which does not satisfy the frequency constraints are eliminated from the data for model ¯tting. This step increased the accuracy of response surface models in the feasible design space. It is found that the response surface models are able to capture the robust optimal regions of design space. The optimal designs show a reduction of 10 percent in the objective function comprising six vibratory hub loads and 1.5 to 80 percent reduction for the individual vibratory forces and moments. This study demonstrates that the second-order response surface models with space ¯lling-designs can be a favorable choice for computationally intensive robust aeroelastic optimization.

The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

Electrochemical, phosphate hydrolysis, DNA binding and DNA cleavage properties of new polyaza macrobicyclic dinickel(II) complexes

A new class of macrobicyclic dinickel(II) complexes Ni2L1,2 B](ClO4)(4) (1-6), where L-1,L-2 are polyaza macrobicyclic binucleating ligands, and B is a N,N-donor heterocyclic base (viz. 2,2'-bipyridine (bipy) and 1,10-phenanthroline (phen)) are synthesized and characterized. The redox, catalytic, DNA binding and DNA cleavage properties were studied. They exhibit two irreversible waves in the cathodic region around E-pc = -0.95 V and E-pa = -0.85 V vs. Ag/Ag+ in CH3CN-0.1 M TBAP, respectively. The first order rate constants for the hydrolysis of 4-nitrophenylphosphate to 4-nitrophenolate by the dinickel(II) complexes 1-6 are in the range from 3.36 x 10(-5) to 10.83 x 10(-5) Ms-1. The complexes 3 and 6 show good binding propensity to calf thymus DNA giving binding constant values (K-b) in the range from 3.08 x 10(5) to 5.37 x 10(5) M-1. The binding site sizes and viscosity data suggest the DNA intercalative and/or groove binding nature of the complexes. The complexes display significant hydrolytic cleavage of supercoiled pBR322DNA at pH 7.2 and 37 degrees C. The hydrolytic cleavage of DNA by the complexes is supported by the evidence from free radical quenching and T4 ligase ligation. The pseudo Michaelis-Menten kinetic parameters k(cat) = 5.44 x 10(-2) h(-1) and K-M = 6.23 x 10(-3) M for complex 3 were obtained. Complex 3 also shows an enormous enhancement of the cleavage rate, of 1.5 x 10(6), in comparison to the uncatalysed hydrolysis rate (k = 3.6 x 10(-8) h(-1)) of ds-DNA.

Thermal degradation kinetics of thermoresponsive poly(N-isopropylacrylamide-co-N,N-dimethylacrylamide) copolymers prepared via RAFT polymerization

Reversible addition-fragmentation chain transfer polymerization at 70 A degrees C in N,N-dimethylformamide was used to prepare poly(N-isopropylacrylamide-co-N,N-dimethylacrylamide) copolymers in various compositions to afford well-defined polymers with pre-determined molecular weight, narrow molecular weight distribution, and precise chain end structure. The copolymer compositions were determined by H-1 NMR spectroscopy. The reactivity ratios of N-isopropylacrylamide (NIPAM) and N,N-dimethylacrylamide (DMA) were calculated as r (NIPAM) = 0.838 and r (DMA) = 1.105, respectively, by the extended Kelen-Tudos method at high conversions. The lower critical solution temperature of PNIPAM can be altered by changing the DMA content in the copolymer chain. Differential scanning calorimetry and thermogravimetric analysis at different heating rates were carried out on these copolymers to understand the nature of thermal degradation and to determine its kinetics. Different kinetic models were applied to estimate various parameters like the activation energy, the order, and the frequency factor. These studies are important to understand the solid state polymer degradation of N-alkyl substituted polymers, which show great potential in the preparation of miscible polymer blends due to their ability to interact through hydrogen bonding.

Correspondence between neutron depolarization and higher order magnetic susceptibility to investigate ferromagnetic clusters in phase separated systems

It is a tough task to distinguish a short-range ferromagnetically correlated cluster-glass phase from a canonical spin-glass-like phase in many magnetic oxide systems using conventional magnetometry measurements. As a case study, we investigate the magnetic ground state of La0.85Sr0.15CoO3, which is often debated based on phase separation issues. We report the results of two samples of La0.85Sr0.15CoO3 (S-1 and S-2) prepared under different conditions. Neutron depolarization, higher harmonic ac susceptibility and magnetic relaxation studies were carried out along with conventional magnetometry measurements to differentiate subtle changes at the microscopic level. There is no evidence of ferromagnetic correlation in the sample S-2 attributed to a spin-glass phase, and this is compounded by the lack of existence of a second order component of higher harmonic ac susceptibility and neutron depolarization. A magnetic relaxation experiment at different temperatures complements the spin glass characteristic in S-2. All these signal a sharp variance when we consider the cluster-glass-like phase (phase separated) in S-1, especially when prepared from an improper chemical synthesis process. This shows that the nonlinear ac susceptibility is a viable tool to detect ferromagnetic clusters such as those the neutron depolarization study can reveal.

Lysis dynamics and membrane oligomerization pathways for Cytolysin A (ClyA) pore-forming toxin

Pore-forming toxins are known for their ability to efficiently form transmembrane pores which eventually leads to cell lysis. The dynamics of lysis and underlying self-assembly or oligomerization pathways leading to pore formation are incompletely understood. In this manuscript the pore-forming kinetics and lysis dynamics of Cytolysin-A (ClyA) toxins on red blood cells (RBCs) are quantified and compared with experimental lysis data. Lysis experiments are carried out on a fixed mass of RBCs, under isotonic conditions in phosphate-buffered saline, for different initial toxin concentrations ranging from 2.94-14.7 nM. Kinetic models which account for monomer binding, conformation and oligomerization to form the dodecameric ClyA pore complex are developed and lysis is assumed to occur when the number of pores per RBC (n(p)) exceeds a critical number, n(pc). By analysing the model in a sublytic regime (n(p) < n(pc)) the number of pores per RBC to initiate lysis is found to lie between 392 and 768 for the sequential oligomerization mechanism and between 5300 and 6300 for the non-sequential mechanism. Rupture rates which are first order in the number of RBCs are seen to provide the best agreement with the lysis experiments. The time constants for pore formation are estimated to lie between 1 and 20 s and monomer conformation time scales were found to be 2-4 times greater than the oligomerization times. Cell rupture takes places in 100s of seconds, and occurs predominantly with a steady number of pores ranging from 515 to 11 000 on the RBC surface for the sequential mechanism. Both the sequential irreversible and non-sequential kinetics provide similar predictions of the hemoglobin release dynamics, however the hemoglobin released as a function of the toxin concentration was accurately captured only with the sequential model. Each mechanism develops a distinct distribution of mers on the surface, providing a unique experimentally observable fingerprint to identify the underlying oligomerization pathways. Our study offers a method to quantify the extent and dynamics of lysis which is an important aspect of developing novel drug and gene delivery strategies based on pore-forming toxins.

Zero miss distance guidance using feedforward and periodic control

There have been attempts at obtaining robust guidance laws to ensure zero miss distance (ZMD) for interceptors with parametric uncertainties. All these laws require the plant to be of minimum phase type to enable the overall guidance loop transfer function to satisfy strict positive realness (SPR). The SPR property implies absolute stability of the closed loop system, and has been shown in the literature to lead to ZMD because it avoids saturation of lateral acceleration. In these works higher order interceptors are reduced to lower order equivalent models for which control laws are designed to ensure ZMD. However, it has also been shown that when the original system with right half plane (RHP) zeros is considered, the resulting miss distances, using such strategies, can be quite high. In this paper, an alternative approach using the circle criterion establishes the conditions for absolute stability of the guidance loop and relaxes the conservative nature of some earlier results arising from assumption of in�nite engagement time. Further, a feedforward scheme in conjunction with a lead-lag compensator is used as one control strategy while a generalized sampled hold function is used as a second strategy, to shift the RHP transmission zeros, thereby achieving ZMD. It is observed that merely shifting the RHP zero(s) to the left half plane reduces miss distances signi�cantly even when no additional controllers are used to ensure SPR conditions.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Faster computation of the Karhunen-Loeve expansion using its domain independence property

The goal of this work is to reduce the cost of computing the coefficients in the Karhunen-Loeve (KL) expansion. The KL expansion serves as a useful and efficient tool for discretizing second-order stochastic processes with known covariance function. Its applications in engineering mechanics include discretizing random field models for elastic moduli, fluid properties, and structural response. The main computational cost of finding the coefficients of this expansion arises from numerically solving an integral eigenvalue problem with the covariance function as the integration kernel. Mathematically this is a homogeneous Fredholm equation of second type. One widely used method for solving this integral eigenvalue problem is to use finite element (FE) bases for discretizing the eigenfunctions, followed by a Galerkin projection. This method is computationally expensive. In the current work it is first shown that the shape of the physical domain in a random field does not affect the realizations of the field estimated using KL expansion, although the individual KL terms are affected. Based on this domain independence property, a numerical integration based scheme accompanied by a modification of the domain, is proposed. In addition to presenting mathematical arguments to establish the domain independence, numerical studies are also conducted to demonstrate and test the proposed method. Numerically it is demonstrated that compared to the Galerkin method the computational speed gain in the proposed method is of three to four orders of magnitude for a two dimensional example, and of one to two orders of magnitude for a three dimensional example, while retaining the same level of accuracy. It is also shown that for separable covariance kernels a further cost reduction of three to four orders of magnitude can be achieved. Both normal and lognormal fields are considered in the numerical studies. (c) 2014 Elsevier B.V. All rights reserved.

Global response sensitivity analysis using probability distance measures and generalization of Sobol's analysis

The study introduces two new alternatives for global response sensitivity analysis based on the application of the L-2-norm and Hellinger's metric for measuring distance between two probabilistic models. Both the procedures are shown to be capable of treating dependent non-Gaussian random variable models for the input variables. The sensitivity indices obtained based on the L2-norm involve second order moments of the response, and, when applied for the case of independent and identically distributed sequence of input random variables, it is shown to be related to the classical Sobol's response sensitivity indices. The analysis based on Hellinger's metric addresses variability across entire range or segments of the response probability density function. The measure is shown to be conceptually a more satisfying alternative to the Kullback-Leibler divergence based analysis which has been reported in the existing literature. Other issues addressed in the study cover Monte Carlo simulation based methods for computing the sensitivity indices and sensitivity analysis with respect to grouped variables. Illustrative examples consist of studies on global sensitivity analysis of natural frequencies of a random multi-degree of freedom system, response of a nonlinear frame, and safety margin associated with a nonlinear performance function. (C) 2015 Elsevier Ltd. All rights reserved.

Correlation equations for average deposition rate coefficients of nanoparticles in a cylindrical pore

Nanoparticle deposition behavior observed at the Darcy scale represents an average of the processes occurring at the pore scale. Hence, the effect of various pore-scale parameters on nanoparticle deposition can be understood by studying nanoparticle transport at pore scale and upscaling the results to the Darcy scale. In this work, correlation equations for the deposition rate coefficients of nanoparticles in a cylindrical pore are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Based on dominant processes, the pore space is divided into three different regions, namely, bulk, diffusion, and potential regions. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile. Among other effects, we account for nanoparticle-collector interaction forces on nanoparticle deposition. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. In this study, pore-scale simulations are performed for three values of Peclet number, Pe = 0.05, 5, and 50. We find that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Peclet numbers (Pe = 0.05) and by a kinetic model at high Peclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relation with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with the colloid filtration theory. These equations can be incorporated into pore network models to study the effect of pore-scale parameters on nanoparticle deposition at larger length scales such as Darcy scale.

Structural evolution and second harmonic properties of lithium niobate-tantalate nanocrystals embedded in a borate glass

Nano-crystals of LiNbxTa1 (-) O-x(3) were evolved by subjecting melt-quenched 1.5Li(2)O-2B(2)O(3)-xNb(2)O(5)-(1 - x)Ta2O5 glasses (where x = 0, 0.25, 0.5, 0.75 and 1.00) to a controlled 3-h isothermal heat treatment between 530 and 560 degrees C. Detailed X-ray diffraction and Raman spectral studies confirmed the formation of nano-crystalline LiNbxTa1 (-) O-x(3) along with a minor phase of ferroelectric and non-linear optic Li2B4O7. The sizes of the nanocrystals evolved in the glass were in the range of 19-37 nm for x = 0-0.75 and 23-45 nm for x = 1.00. Electron microscopic studies confirmed a transformation of the morphology of the nano-crystallites from dendritic star-shaped spherulites for x = 0 to rod-shaped structures for x = 1.00 brought about by a coalescence of crystallites. Broad Maker-fringe patterns (recorded at 532 nm) were obtained by subjecting the heat-treated glass plates to 1064 nm fundamental radiation. However, an effective second order non-linear optic coefficient, d(eff), of 0.45 pm/V, which is nearly 1.2 times the d(36) of KDP single crystal, was obtained for a 560 degrees C/3 h heat-treated glass of the representative composition x = 0.50 comprising 37 nm sized crystallites. (C) 2015 Elsevier B.V. All rights reserved.