Influence of poling methods on the orientational dynamics of 2-methyl-4-nitro-aniline in poly(methyl methacrylate) studied by second harmonic generation

The magnitude and stability of the induced dipolar orientation of 2-methyl-4-nitroaniline (MNA)/poly(methyl methacrylate) (PMMA) guest/host system is investigated. The chromophores are aligned using both the corona discharge and contact electrode poling techniques. The magnitude of order parameter (also an indicator for the second order nonlinear susceptibility) is measured by recording absorbances of the poled (by the two different techniques) and unpoled PMMA films at different concentrations of MNA. Under the same conditions the corona poling technique creates a higher alignment of molecules along the field direction. The time dependence of the second harmonic intensity of the MNA/PMMA film prepared by the two techniques can be described by a Kohlrausch-Williams-Watts stretched exponential. The temperature dependence of the decay time constant is found to generally follow a modified Williams-Landel-Ferry (WLF) or Vogel-Tamann-Fulcher (VTF) equation. The glass transition temperature seems to be the single most important parameter for determining the relaxation time tau(T).

Influence of glycosidic linkage on the nature of carbohydrate binding in beta-prism I fold lectins: An X-ray and molecular dynamics investigation on banana lectin-carbohydrate complexes

The three crystal structures reported here provide details of the interactions of mannose and the mannosyl-alpha-1,3-mannose component of a pentamannose with banana lectin and evidence for the binding of glucosyl-alpha-1,2-glucose to the lectin. The known structures involving the lectin include a complex with glucosyl-beta-1,3-glucose. Modeling studies on the three disaccharide complexes with the reducing end and the nonreducing end at the primary binding site are also provided here. The results of the Xray and modeling studies show that the disaccharides with an alpha-1,3 linkage prefer to have the nonreducing end at the primary binding site, whereas the reducing end is preferred at the site when the linkage is beta-1,3 in mannose/glucose-specific beta-prism I fold lectins. In the corresponding galactose-specific lectins, however, alpha-1,3-linked disaccharides cannot bind the lectin with the nonreducing end at the primary binding site on account of steric clashes with an aromatic residue that occurs only when the lectin is galactose-specific. Molecular dynamics simulations based on the known structures involving banana lectin enrich the information on lectin-carbohydrate interactions obtained from crystal structures. They demonstrate that conformational selection as well as induced fit operate when carbohydrates bind to banana lectin.

Crystal structures, dynamics and functional implications of molybdenum-cofactor biosynthesis protein MogA from two thermophilic organisms

Molybdenum-cofactor (Moco) biosynthesis is an evolutionarily conserved pathway in almost all kingdoms of life, including humans. Two proteins, MogA and MoeA, catalyze the last step of this pathway in bacteria, whereas a single two-domain protein carries out catalysis in eukaryotes. Here, three crystal structures of the Moco-biosynthesis protein MogA from the two thermophilic organisms Thermus thermophilus (TtMogA; 1.64 angstrom resolution, space group P2(1)) and Aquifex aeolicus (AaMogA; 1.70 angstrom resolution, space group P2(1) and 1.90 angstrom resolution, space group P1) have been determined. The functional roles and the residues involved in oligomerization of the protein molecules have been identified based on a comparative analysis of these structures with those of homologous proteins. Furthermore, functional roles have been proposed for the N- and C-terminal residues. In addition, a possible protein-protein complex of MogA and MoeA has been proposed and the residues involved in protein-protein interactions are discussed. Several invariant water molecules and those present at the subunit interfaces have been identified and their possible structural and/or functional roles are described in brief. In addition, molecular-dynamics and docking studies with several small molecules (including the substrate and the product) have been carried out in order to estimate their binding affinities towards AaMogA and TtMogA. The results obtained are further compared with those obtained for homologous eukaryotic proteins.

Nonlinear dynamics of a rigid block on a rigid base

The planar rocking of a prismatic rectangular rigid block about either of its corners is considered. The problem of homoclinic intersections of the stable and unstable manifolds of the perturbed separatrix is addressed to and the corresponding Melnikov functions are derived. Inclusion of the vertical forcing in the Hamiltonian permits the construction of a three-dimensional separatrix. The corresponding modified Melnikov function of Wiggins for homoclinic intersections is derived. Further, the 1-period symmetric orbits are predicted analytically using the method of averaging and compared with the simulation results. The stability boundary for such orbits is also established.

Critical earthquake input power spectral density function models for engineering structures

The problem of determining optimal power spectral density models for earthquake excitation which satisfy constraints on total average power, zero crossing rate and which produce the highest response variance in a given linear system is considered. The solution to this problem is obtained using linear programming methods. The resulting solutions are shown to display a highly deterministic structure and, therefore, fail to capture the stochastic nature of the input. A modification to the definition of critical excitation is proposed which takes into account the entropy rate as a measure of uncertainty in the earthquake loads. The resulting problem is solved using calculus of variations and also within linear programming framework. Illustrative examples on specifying seismic inputs for a nuclear power plant and a tall earth dam are considered and the resulting solutions are shown to be realistic.

Quantum and non-Markovian effects in the electron transfer reaction dynamics in the Marcus inverted region

Electron transfer reactions in large molecules may often be coupled to both the polar solvent modes and the intramolecular vibrational modes of the molecule. This can give rise to a complex dynamics which may in some systems, like betaine, be controlled more by vibrational rather than by solvent effects. Additionally, a significant contribution from an ultrafast relaxation component in the solvation dynamics may enhance the complexity. To explain the wide range of behavior that has been observed experimentally, Barbara et al. recently proposed that a model of an electron transfer reaction should minimally consist of a low-frequency classical solvent mode (X), a low-frequency vibrational mode (Q), and a high-frequency quantum mode (q) (J. Phys. Chem. 1991, 96, 3728). In the present work, a theoretical study of this model is described. This study generalizes earlier work by including the biphasic solvent response and the dynamics of the low-frequency vibrational mode in the presence of a delocalized, extended reaction zone. A novel Green's function technique has been developed which allowed us to study the non-Markovian dynamics on a multidimensional surface. The contributions from the high-frequency vibrational mode and the ultrafast component in the non-Markovian solvent dynamics are found to be primarily responsible for the dramatic increase in charge transfer rate over the prediction of the classical theories that neglect both these factors. These, along with a large coupling between the reactant and the product states, may combine to render the electron transfer rate both very large and constant over a wide range of solvent relaxation rates. A study on the free energy gap dependence of the electron transfer rate reveals that the rates are sensitive to changes in the quantum frequency particularly when the free energy gap is very large.

Self-consistent microscopic treatment of the effects of self-motion of the probe on ionic and dipolar solvation dynamics

A simple but self-consistent microscopic theory for the time dependent solvation energy of both ions and dipoles is presented which includes, for the first time, the details of the self-motion of the probe on its own solvation dynamics. The theory leads to several interesting predictions. The most important of them is that, for dipolar solvation, both the rotational and the translational motions of the dipolar solute probe can significantly accelerate the rate of solvation. In addition, the rotational self-motion of the solute can also give rise to an additional mechanism of nonexponentiality in solvation time correlation functions in otherwise slow liquids. A comparison between the present theoretical predictions and the recent experimental studies of Maroncelli et al. on solvation dynamics of aniline in l-propanol seems to indicate that the said experiments have missed the initial solvent response up to about 45 ps. After mapping the experimental results on the redefined time scale, the theoretical results can explain the experimental results for solvation of aniline in 1-propanol very well. For ionic solvation, the translational motion is significant for light solutes only. For example, for Li+ in water, translational motion speeds up the solvation by about 20%. The present theory demonstrates that in dipolar solvation the partial quenching of the self-motion due to the presence of specific solute-solvent interactions (such as H-bonding) may lead to a much slower solvation than that when the self-motion is present. This point has been discussed. In addition, we present the theoretical results for solvation of aniline in propylene carbonate, Here, the solvation is predicted to be complete within 15-20 ps.

Solvation dynamics in slow, viscous liquids: Application to amides

As the viscosity of a liquid increases rapidly in the supercooled regime, the nature of molecular relaxation can exhibit dynamics rather different from the fast dynamics observed in the normal regime. In this article, we present theoretical studies of solvation dynamics and orientational relaxation in slow liquids. As the local short-range correlations are important in the slow liquids, we have extended our previous theory to take into account the shea-range pair correlations between the polar solute and the dipolar solvent molecules. Application of the generalized theory To the study of solvation dynamics of amide systems gives nice agreement with the experimental results of Maroncelli and co-workers (J. Phys. Chem. 1990, 94, 4929). The theory also provides valuable insight into the orientational relaxation precesses in the viscous liquids.

Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Stability of stochastic Leipholz column with stochastic loading

The Leipholz column which is having the Young modulus and mass per unit length as stochastic processes and also the distributed tangential follower load behaving stochastically is considered. The non self-adjoint differential equation and boundary conditions are considered to have random field coefficients. The standard perturbation method is employed. The non self-adjoint operators are used within the regularity domain. Full covariance structure of the free vibration eigenvalues and critical loads is derived in terms of second order properties of input random fields characterizing the system parameter fluctuations. The mean value of critical load is calculated using the averaged problem and the corresponding eigenvalue statistics are sought. Through the frequency equation a transformation is done to yield load parameter statistics. A numerical study incorporating commonly observed correlation models is reported which illustrates the full potentials of the derived expressions.

Stochastic differential games: Occupation measure based approach

A new approach based on occupation measures is introduced for studying stochastic differential games. For two-person zero-sum games, the existence of values and optimal strategies for both players is established for various payoff criteria. ForN-person games, the existence of equilibria in Markov strategies is established for various cases.

Stochastic fluctuations of the solar dynamo

Attempts in the past to model the irregularities of the solar cycle (such as the Maunder minimum) were based on studies of the nonlinear feedback of magnetic fields on the dynamo source terms. Since the alpha-coefficient is obtained by averaging over the turbulence, it is expected to have stochastic fluctuations, and we show that these fluctuations can explain the irregularities of the solar cycle in a more satisfactory way. We solve the dynamo equations in a slab with a single mode, taking the alpha-coefficient to be constant in space but fluctuating stochastically in time with some given amplitude and given correlation time. The same level of percentile fluctuations (about 10 %) produces no effect on an alpha-omega dynamo, but makes an alpha-2 dynamo completely chaotic. The level of irregularities in an alpha-2-omega dynamo qualitatively agrees with the solar behavior, reinforcing the conclusion of Choudhuri (1990a) that the solar dynamo is of the alpha-2-omega-type. The irregularities are found to increase on increasing either the amplitude or the correlation time of the stochastic fluctuations. The alpha-quenching mechanism tends to make the system stable against the irregularities and hence it is inferred that the alpha-quenching should not be too strong so that the irregularities are not completely suppressed. We also present a simple-minded analysis to understand why the stochastic fluctuations in the alpha-omega, alpha-2-omega and alpha-2 regimes have such different outcomes.

Stochastic modelling of flexible manufacturing systems

Mathematical modelling plays a vital role in the design, planning and operation of flexible manufacturing systems (FMSs). In this paper, attention is focused on stochastic modelling of FMSs using Markov chains, queueing networks, and stochastic Petri nets. We bring out the role of these modelling tools in FMS performance evaluation through several illustrative examples and provide a critical comparative evaluation. We also include a discussion on the modelling of deadlocks which constitute an important source of performance degradation in fully automated FMSs.