Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Decoupling of diffusion from viscosity: Difference scenario for translational and rotational motions

Recent experiments have indicated a dramatically different viscosity dependence of the translational and the rotational diffusion coefficients in a supercooled liquid as the glass transition temperature is approached from above. While the translational motion seems to be decoupled from the rising viscosity (eta), the rotational motion seems to remain firmly coupled to eta. In order to understand the microscopic origin of this behavior, we have carried nut detailed theoretical calculations of both the quantities by using a self-consistent mode-coupling theory (MCT). it is found that when the size of the solute is same as that of the solvent molecules, the conventional MCT fails to predict the observed decoupling. The solvent inhomogeneity is found to play a decisive role in determining the decoupling. The difference in the viscosity dependence between rotation and translational diffusion coefficient is discussed.

A Correlation between the Dynamic Viscosity and Crystallisation of Felodipine from its Amorphous Polymer Solid Dispersion

Purpose Previously, it has been reported that molecular mobility determines the rate of molecular approach to crystal surfaces, while entropy relates to the probability of that approaching molecule having the desirable configuration for further growth of the existing crystal; and the free energy dictates the probability of that molecule not returning to the liquid phase1. If we plot the crystal growth rate and viscosity of a supercooled liquid in a log-log format, the relationship between the two is linear, indicating the influence viscosity has upon crystal growth rate. However, such approximation has been derived from pure drug compounds and it is apparent that further understanding of crystallization from drug-polymer solid dispersion is required in order to stabilise drugs embedded within amorphous polymeric solid dispersions. Methods Mixtures of felodipine and polymer (HPMCAS-HF, PVPK15 and Soluplus®) at specified compositions were prepared using a Restch MM200 ball mill. To examine crystal growth within amorphous solid dispersions, samples were prepared by melting 5-10 mg of ball milled mixture at 150°C for 3-5 minutes on a glass slip pre-cleaned with methanol and acetone. All prepared samples were confirmed to be crystal free by visual observation using a polarised light microscope (Olympus BX50). Prepared samples were stored at 0% RH (P2O5), inside desiccators, maintained in ovens at 80°C. For the dynamic viscosity measurement, approximately 100-200mg ball milled mixture was heated on the base plate of a rotational rheometer at 150°C for 5 minutes and the top plate was lowered to a defined gap to form a good contact with the material. The sandwiched amorphous material was heated to 80°C and the viscosity was measured. Results The equation was used to probe the correlation of viscosity to crystal growth rate. In comparison to the value of xi in log-log equation reported from pure drug compound, a value of 1.63 was obtained for FD-polymer solid dispersions irrespective of the polymer involved. &#8733 Conclusion The high xi value suggests stronger viscosity dependence may exist for amorphous FD once incorporated with amorphous polymer.

Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Energy transfer efficiency distributions in polymers in solution during folding and unfolding

Distribution of fluorescence resonance energy transfer (FRET) efficiency between the two ends of a Lennard-Jones polymer chain both at equilibrium and during folding and unfolding has been calculated, for the first time, by Brownian dynamics simulations. The distribution of FRET efficiency becomes bimodal during folding of the extended state subsequent to a temperature quench, with the width of the distribution for the extended state broader than that for the folded state. The reverse process of unfolding subsequent to a upward temperature jump shows different characteristics. The distributions show significant viscosity dependence which can be tested against experiments.

Experimental Evidence of the Occurrence of Intramolecular Electron Transfer in Catalyzed 1,2-Dioxetane Decomposition

The activation parameters for the thermal decomposition of 13 acridinium-substituted 1,2-dioxetanes, bearing an aromatic moiety, were determined and their chemiluminescence emission quantum yields estimated, utilizing in situ photosensitized 1,2-dioxetane generation and observation of its thermal decomposition kinetics, without isolation of these highly unstable cyclic peroxides. Decomposition rate constants show linear free-energy correlation for electron-withdrawing substituents, with a Hammett reaction constant of rho = 1.3 +/- 0.1, indicating the occurrence of an intramolecular electron transfer from the acridinium moiety to the 1,2-dioxetane ring, as postulated by the intramolecular chemically initiated electron exchange luminescence (CIEEL) mechanism. Emission quantum yield behavior can also be rationalized on the basis of the intramolecular CIEEL mechanism, additionally evidencing its occurrence in this transformation. Both relations constitute the first experimental evidence for the occurrence of the postulated intramolecular electron transfer in the catalyzed and induced decomposition of properly substituted 1,2-dioxetanes.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Barriers to protein folding: formation of buried polar interactions is a slow step in acquisition of structure.

In the MYL mutant of the Arc repressor dimer, sets of partially buried salt-bridge and hydrogen-bond interactions mediated by Arg-31, Glu-36, and Arg-40 in each subunit are replaced by hydrophobic interactions between Met-31, Tyr-36, and Leu-40. The MYL refolding/dimerization reaction differs from that of wild type in being 10- to 1250-fold faster, having an earlier transition state, and depending upon viscosity but not ionic strength. Formation of the wild-type salt bridges in a hydrophobic environment clearly imposes a kinetic barrier to folding, which can be lowered by high salt concentrations. The changes in the position of the transition state and viscosity dependence can be explained if denatured monomers interact to form a partially folded dimeric intermediate, which then continues folding to form the native dimer. The second step is postulated to be rate limiting for wild type. Replacing the salt bridge with hydrophobic interactions lowers this barrier for MYL. This makes the first kinetic barrier rate limiting for MYL refolding and creates a downhill free-energy landscape in which most molecules which reach the intermediate state continue to form native dimers.

Viscosity of aged bio-oils from fast pyrolysis of beech wood and miscanthus:shear rate and temperature dependence

The viscosity of four aged bio-oil samples was measured experimentally at various shear rates and temperatures using a rotational viscometer. The experimental bio-oils were derived from fast pyrolysis of beech wood at 450, 500, and 550 °C and Miscanthus at 500 °C (in this work, they were named as BW1, BW2, BW3, and MXG) in a bubbling fluidized bed reactor. The viscosity of all bio-oils was kept constant at various shear rates at the same temperature, which indicated that they were Newtonian fluids. The viscosity of bio-oils was strongly dependent upon the temperature, and with the increase of the temperature from 30 to 80 °C, the viscosity of BW1, BW2, BW3, and MXG decreased by 90.7, 93.3, 92.6, and 90.2%, respectively. The Arrhenius viscosity model, which has been commonly used to represent the temperature dependence of the viscosity of many fluids, did not fit the viscosity-temperature experimental data of all bio-oils very well, especially in the low- and high-temperature regions. For comparison, the Williams-Landel-Ferry (WLF) model was also used. The results showed that the WLF model gave a very good description of the viscosity-temperature relationship of each bio-oil with very small residuals and the BW3 bio-oil had the strongest viscosity-temperature dependence.

An Artificial Neutral Network (ANN) model for predicting biodiesel kinetic viscosity as a function of temperature and chemical compositions

An Artificial Neural Network (ANN) is a computational modeling tool which has found extensive acceptance in many disciplines for modeling complex real world problems. An ANN can model problems through learning by example, rather than by fully understanding the detailed characteristics and physics of the system. In the present study, the accuracy and predictive power of an ANN was evaluated in predicting kinetic viscosity of biodiesels over a wide range of temperatures typically encountered in diesel engine operation. In this model, temperature and chemical composition of biodiesel were used as input variables. In order to obtain the necessary data for model development, the chemical composition and temperature dependent fuel properties of ten different types of biodiesels were measured experimentally using laboratory standard testing equipments following internationally recognized testing procedures. The Neural Networks Toolbox of MatLab R2012a software was used to train, validate and simulate the ANN model on a personal computer. The network architecture was optimised following a trial and error method to obtain the best prediction of the kinematic viscosity. The predictive performance of the model was determined by calculating the absolute fraction of variance (R2), root mean squared (RMS) and maximum average error percentage (MAEP) between predicted and experimental results. This study found that ANN is highly accurate in predicting the viscosity of biodiesel and demonstrates the ability of the ANN model to find a meaningful relationship between biodiesel chemical composition and fuel properties at different temperature levels. Therefore the model developed in this study can be a useful tool in accurately predict biodiesel fuel properties instead of undertaking costly and time consuming experimental tests.

Effect of dynamical asymmetry on the viscosity of a random copolymer melt

The variation of the viscosity as a function of the sequence distribution in an A-B random copolymer melt is determined. The parameters that characterize the random copolymer are the fraction of A monomers f, the parameter lambda which determines the correlation in the monomer identities along a chain and the Flory chi parameter chi(F) which determines the strength of the enthalpic repulsion between monomers of type A and B. For lambda>0, there is a greater probability of finding like monomers at adjacent positions along the chain, and for lambda<0 unlike monomers are more likely to be adjacent to each other. The traditional Markov model for the random copolymer melt is altered to remove ultraviolet divergences in the equations for the renormalized viscosity, and the phase diagram for the modified model has a binary fluid type transition for lambda>0 and does not exhibit a phase transition for lambda<0. A mode coupling analysis is used to determine the renormalization of the viscosity due to the dependence of the bare viscosity on the local concentration field. Due to the dissipative nature of the coupling. there are nonlinearities both in the transport equation and in the noise correlation. The concentration dependence of the transport coefficient presents additional difficulties in the formulation due to the Ito-Stratonovich dilemma, and there is some ambiguity about the choice of the concentration to be used while calculating the noise correlation. In the Appendix, it is shown using a diagrammatic perturbation analysis that the Ito prescription for the calculation of the transport coefficient, when coupled with a causal discretization scheme, provides a consistent formulation that satisfies stationarity and the fluctuation dissipation theorem. This functional integral formalism is used in the present analysis, and consistency is verified for the present problem as well. The upper critical dimension for this type of renormaliaation is 2, and so there is no divergence in the viscosity in the vicinity of a critical point. The results indicate that there is a systematic dependence of the viscosity on lambda and chi(F). The fluctuations tend to increase the viscosity for lambda<0, and decrease the viscosity for lambda>0, and an increase in chi(F) tends to decrease the viscosity. (C) 1996 American Institute of Physics.

Saturated liquid viscosity correlations for alternative refrigerants

This article presents dimensionless equations for the temperature dependence of the saturated liquid viscosity of R32, R123, R124, R125, R134a, R141b, and R152a valid over a temperature range of engineering interest. The correlation has the form Phi(D)(n)=A+BTD where Phi(D) is the dimensionless fluidity (1/eta(D)) and T-D is a dimensionless temperature. n, A, and B are evaluated for each of the above refrigerants based on a least-squares fit to experimental data. This equation is found to provide an improved fit over those existing in the literature up to T-D=0.8.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.