Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Structural investigation of resorcinol based symmetrical banana mesogens by XRD, NMR and polarization measurements

Synthesis and structural characterization of two novel symmetrical banana mesogens built from resorcinol with seven phenyl rings linked by ester and imine with a terminal dodecyl/dodecyloxy chain has been carried out. Density functional theory (DFT) has been employed for obtaining the geometry optimized structures, the dipole moments and C-13 NMR chemical shifts. The HOPM and DSC studies revealed enantiotropic B-2 and B-7 phases for the dodecyl and dodecyloxy homologs respectively. The powder X-ray studies of both the mesogens indicate the presence of layer ordering. The polarization measurements reveal an anti-ferroelectric switching for the B-2 phase of the dodecyl homolog whose structure has been identified as SmCSPA. The B-7 phase of the dodecyloxy homolog was found to be non-switchable. High resolution C-13 NMR study of the dodecyl homolog in its mesophase has been carried out. C-13-H-1 dipolar couplings obtained from the 2-dimensional separated local field spectroscopy experiment were used to obtain the orientational order parameters of the different segments of the mesogen. Very large C-13-H-1 dipolar couplings observed for the carbons of the central phenyl ring (9.7-12.3 kHz) in comparison to the dipolar couplings of those of the side arm phenyl rings (less than 3 kHz) are a direct consequence of the ordering in the banana phase and the shape of the molecule. From the ratio of the local order parameter values, the bent-angle of the mesogen could be determined in a straight forward manner to be 120.5 degrees.

Investigation of dielectric relaxation, Jahn-Teller distortion, and magnetic ordering in Y substituted Pr1-xYxMnO3 (0.1 <= x <= 0.4)

We report structural, magnetic, and dielectric properties of the perovskite compound Pr1-xYxMnO3 (0.1 <= x <= 0.4) studied using dc magnetization, ac susceptibility, neutron powder diffraction, and dielectric techniques. These compounds crystallize in orthorhombic space group (Pnma) in the temperature range 5-300 K. The Mn-O-Mn bond angle decreases with the Y substitution along with an increase in the Jahn-Teller distortion. The Jahn-Teller distortion for Pr0.9Y0.1MnO3 shows an anomalous change near 50 K, below which it falls sharply. Neutron powder diffraction patterns of all reported compositions at low temperature constitute additional magnetic Bragg peaks that suggest magnetic ordering. Magnetic reflections were indexed in the nuclear lattice with the propagation vector k = (0, 0, 0). Rietveld refinement of powder patterns conform to A type antiferromagnetic ordering where moments are aligned ferromagnetically in a-c plane and coupled nearly antiferromagnetically along b-axis resulting in a net ferromagnetic component along the b-direction. The antiferromagnetic transition temperature was deduced from dc magnetization and ac susceptibility data. The transition temperature decreases by nearly 22 K (from 81 K to 59 K) as yttrium content (x) increases from 0.1 to 0.4. Measurements reveal strong frequency dispersion in dielectric constant and dielectric loss. Activation energy and relaxation time are estimated from the Arrhenius plot. It is further shown that relaxation behaviour is altered with yttrium doping concentration. (C) 2015 AIP Publishing LLC.

Adaptive Mesh Refinement for Fast Convergence of EFIE-Based 3-D Extraction

3-D full-wave method of moments (MoM) based electromagnetic analysis is a popular means toward accurate solution of Maxwell's equations. The time and memory bottlenecks associated with such a solution have been addressed over the last two decades by linear complexity fast solver algorithms. However, the accurate solution of 3-D full-wave MoM on an arbitrary mesh of a package-board structure does not guarantee accuracy, since the discretization may not be fine enough to capture spatial changes in the solution variable. At the same time, uniform over-meshing on the entire structure generates a large number of solution variables and therefore requires an unnecessarily large matrix solution. In this paper, different refinement criteria are studied in an adaptive mesh refinement platform. Consequently, the most suitable conductor mesh refinement criterion for MoM-based electromagnetic package-board extraction is identified and the advantages of this adaptive strategy are demonstrated from both accuracy and speed perspectives. The results are also compared with those of the recently reported integral equation-based h-refinement strategy. Finally, a new methodology to expedite each adaptive refinement pass is proposed.

Third post-Newtonian gravitational waveforms for compact binary systems in general orbits: Instantaneous terms

We compute the instantaneous contributions to the spherical harmonic modes of gravitational waveforms from compact binary systems in general orbits up to the third post-Newtonian (PN) order. We further extend these results for compact binaries in quasielliptical orbits using the 3PN quasi-Keplerian representation of the conserved dynamics of compact binaries in eccentric orbits. Using the multipolar post-Minkowskian formalism, starting from the different mass and current-type multipole moments, we compute the spin-weighted spherical harmonic decomposition of the instantaneous part of the gravitational waveform. These are terms which are functions of the retarded time and do not depend on the history of the binary evolution. Together with the hereditary part, which depends on the binary's dynamical history, these waveforms form the basis for construction of accurate templates for the detection of gravitational wave signals from binaries moving in quasielliptical orbits.

COMPUTATIONAL MODELING OF A MULTI-LAYERED PIEZO-COMPOSITE BEAM MADE UP OF MFC

The Variational Asymptotic Method (VAM) is used for modeling a coupled non-linear electromechanical problem finding applications in aircrafts and Micro Aerial Vehicle (MAV) development. VAM coupled with geometrically exact kinematics forms a powerful tool for analyzing a complex nonlinear phenomena as shown previously by many in the literature 3 - 7] for various challenging problems like modeling of an initially twisted helicopter rotor blades, matrix crack propagation in a composite, modeling of hyper elastic plates and various multi-physics problems. The problem consists of design and analysis of a piezocomposite laminate applied with electrical voltage(s) which can induce direct and planar distributed shear stresses and strains in the structure. The deformations are large and conventional beam theories are inappropriate for the analysis. The behavior of an elastic body is completely understood by its energy. This energy must be integrated over the cross-sectional area to obtain the 1-D behavior as is typical in a beam analysis. VAM can be used efficiently to approximate 3-D strain energy as closely as possible. To perform this simplification, VAM makes use of thickness to width, width to length, width multiplied by initial twist and strain as small parameters embedded in the problem definition and provides a way to approach the exact solution asymptotically. In this work, above mentioned electromechanical problem is modeled using VAM which breaks down the 3-D elasticity problem into two parts, namely a 2-D non-linear cross-sectional analysis and a 1-D non-linear analysis, along the reference curve. The recovery relations obtained as a by-product in the cross-sectional analysis earlier are used to obtain 3-D stresses, displacements and velocity contours. The piezo-composite laminate which is chosen for an initial phase of computational modeling is made up of commercially available Macro Fiber Composites (MFCs) stacked together in an arbitrary lay-up and applied with electrical voltages for actuation. The expressions of sectional forces and moments as obtained from cross-sectional analysis in closed-form show the electro-mechanical coupling and relative contribution of electric field in individual layers of the piezo-composite laminate. The spatial and temporal constitutive law as obtained from the cross-sectional analysis are substituted into 1-D fully intrinsic, geometrically exact equilibrium equations of motion and 1-D intrinsic kinematical equations to solve for all 1-D generalized variables as function of time and an along the reference curve co-ordinate, x(1).

Exploring CP violation in the MSSM

We explore the prospects for observing CP violation in the minimal supersymmetric extension of the Standard Model (MSSM) with six CP-violating parameters, three gaugino mass phases and three phases in trilinear soft supersymmetry-breaking parameters, using the CPsuperH code combined with a geometric approach to maximise CP-violating observables subject to the experimental upper bounds on electric dipole moments. We also implement CP-conserving constraints from Higgs physics, flavour physics and the upper limits on the cosmological dark matter density and spin-independent scattering. We study possible values of observables within the constrained MSSM (CMSSM), the non-universal Higgs model (NUHM), the CPX scenario and a variant of the phenomenological MSSM (pMSSM). We find values of the CP-violating asymmetry A(CP) in b -> s gamma decay that may be as large as 3 %, so future measurements of ACP may provide independent information about CP violation in the MSSM. We find that CP-violating MSSM contributions to the B-s meson mass mixing term Delta M-Bs are in general below the present upper limit, which is dominated by theoretical uncertainties. If these could be reduced, Delta M-Bs could also provide an interesting and complementary constraint on the six CP-violating MSSM phases, enabling them all to be determined experimentally, in principle. We also find that CP violation in the h(2,3)tau(+)tau(-) and h(2,3) (t) over bart couplings can be quite large, and so may offer interesting prospects for future pp, e(+) e(-), mu(+) mu(-) and gamma gamma colliders.

Effect of Manganese (II) Oxide on microstructure and ionic transport properties of nanostructured cubic zirconia

Effect of MnO addition on microstructure and ionic transport properties of nanocrystalline cubic(c)-ZrO2 is reported. Monoclinic (m) ZrO2 powders with 10-30 mol% MnO powder are mechanically alloyed in a planetary ball mill at room temperature for 10 h and annealed at 550 degrees C for 6 h. In all compositions m-ZrO2 transforms completely to nanocrystalline c-ZrO2 phase and MnO is fully incorporated into c-ZrO2 lattice. Rietveld's refinement technique is employed for detailed microstructure analysis by analyzing XRD patterns. High resolution transmission electron microscopy (HRTEM) analysis confirms the complete formation of c-ZrO2 phase. Presence of stoichiometric Mn in c-ZrO2 powder is confirmed by Electron Probe Microscopy analysis. XPS analysis reveals that Mn is mostly in Mn2+ oxidation state. A correlation between lattice parameter and oxygen vacancy is established. A detailed ionic conductivity measurement in the 250 degrees-575 degrees C temperature range describes the effect of MnO on conductivity of c-ZrO2. The ionic conductivity (s) of 30 mol% MnO alloyed ZrO2 at 550 degrees C is 0.04 s cm(-1). Electrical relaxation studies are carried out by impedance and modulus spectroscopy. Relaxation frequency is found to increase with temperature and MnO mol fraction. Electrical characterization predicts that these compounds have potentials for use as solid oxide fuel cell electrolyte material. (C) 2015 Elsevier Ltd. All rights reserved.

Structure of an amidohydrolase, SACOL0085, from methicillin-resistant Staphylococcus aureus COL

Staphylococcus aureus is an opportunistic pathogen that rapidly acquires resistance to frontline antibiotics. The characterization of novel protein targets from this bacterium is thus an important step towards future therapeutic strategies. Here, the crystal structure of an amidohydrolase, SACOL0085, from S. aureus COL is described. SACOL0085 is a member of the M20D family of peptidases. Unlike other M20D peptidases, which are either monomers or dimers, SACOL0085 adopts a butterfly-shaped homotetrameric arrangement with extensive intersubunit interactions. Each subunit of SACOL0085 contains two Mn2+ ions at the active site. A conserved cysteine residue at the active site distinguishes M20D peptidases from other M20 family members. This cysteine, Cys103, serves as bidentate ligand coordinating both Mn2+ ions in SACOL0085.

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.

A recursive multi-scaling approach to regional flood frequency analysis

Scaling approaches are widely used by hydrologists for Regional Frequency Analysis (RFA) of floods at ungauged/sparsely gauged site(s) in river basins. This paper proposes a Recursive Multi-scaling (RMS) approach to RFA that overcomes limitations of conventional simple- and multi-scaling approaches. The approach involves identification of a separate set of attributes corresponding to each of the sites (being considered in the study area/region) in a recursive manner according to their importance, and utilizing those attributes to construct effective regional regression relationships to estimate statistical raw moments (SMs) of peak flows. The SMs are then utilized to arrive at parameters of flood frequency distribution and quantile estimate(s) corresponding to target return period(s). Effectiveness of the RMS approach in arriving at flood quantile estimates for ungauged sites is demonstrated through leave-one-out cross-validation experiment on watersheds in Indiana State, USA. Results indicate that the approach outperforms index-flood based Region-of-Influence approach, simple- and multi-scaling approaches and a multiple linear regression method. (C) 2015 Elsevier B.V. All rights reserved.

Evaluation of the Index-Flood Approach Related Regional Frequency Analysis Procedures

Index-flood related regional frequency analysis (RFA) procedures are in use by hydrologists to estimate design quantiles of hydrological extreme events at data sparse/ungauged locations in river basins. There is a dearth of attempts to establish which among those procedures is better for RFA in the L-moment framework. This paper evaluates the performance of the conventional index flood (CIF), the logarithmic index flood (LIF), and two variants of the population index flood (PIF) procedures in estimating flood quantiles for ungauged locations by Monte Carlo simulation experiments and a case study on watersheds in Indiana in the U.S. To evaluate the PIF procedure, L-moment formulations are developed for implementing the procedure in situations where the regional frequency distribution (RFD) is the generalized logistic (GLO), generalized Pareto (GPA), generalized normal (GNO) or Pearson type III (PE3), as those formulations are unavailable. Results indicate that one of the variants of the PIF procedure, which utilizes the regional information on the first two L-moments is more effective than the CIF and LIF procedures. The improvement in quantile estimation using the variant of PIF procedure as compared with the CIF procedure is significant when the RFD is a generalized extreme value, GLO, GNO, or PE3, and marginal when it is GPA. (C) 2015 American Society of Civil Engineers.

Thermal, magnetic and electrical properties of Tb1-xDyxMnO3 multiferroics

A series of multiferroic materials with the compositional formula, Tb1 - xDyxMnO3 (where x=0, 0.1, 0.2, 0.3 and 0.4) were prepared by the sol gel method. After characterizing the samples structurally, a systematic investigation of specific heat, magnetization and dielectric properties over the temperature range, 4-300 K, was undertaken. Based on these studies, it was found that all the samples exhibit a transition at 40 K and the observed behavior may be attributed to the ordering of Mn3+ ions. Further, all the five samples are found to exhibit a ferroelectric transition in the temperature range 20-24 K. Finally, yet another transition was also exhibited by all the samples at temperatures below 10 K and is attributed to the antiferromagnetic (AF) ordering of rare-earth ionic moments. The magnetic entropy of all the samples was also computed with the help of their heat capacity data. (C) 2015 Elsevier B.V. All rights reserved.

Learning-Based Fast Iterative Convergence of 3-D MoM via Eigen-AGMRES Method

Three-dimensional (3-D) full-wave electromagnetic simulation using method of moments (MoM) under the framework of fast solver algorithms like fast multipole method (FMM) is often bottlenecked by the speed of convergence of the Krylov-subspace-based iterative process. This is primarily because the electric field integral equation (EFIE) matrix, even with cutting-edge preconditioning techniques, often exhibits bad spectral properties arising from frequency or geometry-based ill-conditioning, which render iterative solvers slow to converge or stagnate occasionally. In this communication, a novel technique to expedite the convergence of MoMmatrix solution at a specific frequency is proposed, by extracting and applying Eigen-vectors from a previously solved neighboring frequency in an augmented generalized minimum residual (AGMRES) iterative framework. This technique can be applied in unison with any preconditioner. Numerical results demonstrate up to 40% speed-up in convergence using the proposed Eigen-AGMRES method.

The role of intra- and inter-site exchange correlations in the extended Falicov-Kimball model on a triangular lattice

Ground state magnetic properties of the spin-dependent Falicov-Kimball model (FKM) are studied by incorporating the intrasite exchange correlation J (between itinerant d- and localized f-electrons) and intersite (superexchange) correlation J (between localized f-electrons) on a triangular lattice for two different fillings. Numerical diagonalization and Monte-Carlo techniques are used to determine the ground state magnetic properties. Transitions from antiferromagnetic to ferromagnetic and again to re-entrant antiferromagnetic phase is observed in a wide range of parameter space. The magnetic moments of d- and f-electrons are observed to depend strongly on the value off, J and also on the total number of d-electrons (N-d). (C) 2015 Elsevier Ltd. All rights reserved.