Molecular computations for reactions and phase transitions: applications to protein stabilization, hydrates and catalysis

In this work we have made significant contributions in three different areas of interest: therapeutic protein stabilization, thermodynamics of natural gas clathrate-hydrates, and zeolite catalysis. In all three fields, using our various computational techniques, we have been able to elucidate phenomena that are difficult or impossible to explain experimentally. More specifically, in mixed solvent systems for proteins we developed a statistical-mechanical method to model the thermodynamic effects of additives in molecular-level detail. It was the first method demonstrated to have truly predictive (no adjustable parameters) capability for real protein systems. We also describe a novel mechanism that slows protein association reactions, called the “gap effect.” We developed a comprehensive picture of methioine oxidation by hydrogen peroxide that allows for accurate prediction of protein oxidation and provides a rationale for developing strategies to control oxidation. The method of solvent accessible area (SAA) was shown not to correlate well with oxidation rates. A new property, averaged two-shell water coordination number (2SWCN) was identified and shown to correlate well with oxidation rates. Reference parameters for the van der Waals Platteeuw model of clathrate-hydrates were found for structure I and structure II. These reference parameters are independent of the potential form (unlike the commonly used parameters) and have been validated by calculating phase behavior and structural transitions for mixed hydrate systems. These calculations are validated with experimental data for both structures and for systems that undergo transitions from one structure to another. This is the first method of calculating hydrate thermodynamics to demonstrate predictive capability for phase equilibria, structural changes, and occupancy in pure and mixed hydrate systems. We have computed a new mechanism for the methanol coupling reaction to form ethanol and water in the zeolite chabazite. The mechanism at 400°C proceeds via stable intermediates of water, methane, and protonated formaldehyde.

Quantitative model for the kinetics of lyotropic phase transitions involving changes in monolayer curvature

We present a kinetic model for transformations between different self-assembled lipid structures. The model shows how data on the rates of phase transitions between mesophases of different geometries can be used to provide information on the mechanisms of the transformations and the transition states involved. This can be used, for example, to gain an insight into intermediate structures in cell membrane fission or fusion. In cases where the monolayer curvature changes on going from the initial to the final mesophase, we consider the phase transition to be driven primarily by the change in the relaxed curvature with pressure or temperature, which alters the relative curvature elastic energies of the two mesophase structures. Using this model, we have analyzed previously published kinetic data on the inter-conversion of inverse bicontinuous cubic phases in the 1-monoolein-30 wt% water system. The data are for a transition between QII(G) and QII(D) phases, and our analysis indicates that the transition state more closely resembles the QII(D) than the QII(G) phase. Using estimated values for the monolayer mean curvatures of the QII(G) and QII(D) phases of -0.123 nm(-1) and -0.133 nm(-1), respectively, gives values for the monolayer mean curvature of the transition state of between -0.131 nm(-1) and -0.132 nm(-1). Furthermore, we estimate that several thousand molecules undergo the phase transition cooperatively within one "cooperative unit", equivalent to 1-2 unit cells of QII(G) or 4-10 unit cells of QII(D).

Assessment of the RE(OH)3 Ising magnetic materials as possible candidates for the study of transverse-field-induced quantum phase transitions

The LiHoxY1−xF4 Ising magnetic material subject to a magnetic field perpendicular to the Ho3+ Ising direction has shown over the past 20 years to be a host of very interesting thermodynamic and magnetic phenomena. Unfortunately, the availability of other magnetic materials other than LiHoxY1−xF4 that may be described by a transverse-field Ising model remains very much limited. It is in this context that we use here a mean-field theory to investigate the suitability of the Ho(OH)3, Dy(OH)3, and Tb(OH)3 insulating hexagonal dipolar Ising-type ferromagnets for the study of the quantum phase transition induced by a magnetic field, Bx, applied perpendicular to the Ising spin direction. Experimentally, the zero-field critical (Curie) temperatures are known to be Tc≈2.54, 3.48, and 3.72 K, for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. From our calculations we estimate the critical transverse field, Bxc, to destroy ferromagnetic order at zero temperature to be Bxc=4.35, 5.03, and 54.81 T for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. We find that Ho(OH)3, similarly to LiHoF4, can be quantitatively described by an effective S=1/2 transverse-field Ising model. This is not the case for Dy(OH)3 due to the strong admixing between the ground doublet and first excited doublet induced by the dipolar interactions. Furthermore, we find that the paramagnetic (PM) to ferromagnetic (FM) transition in Dy(OH)3 becomes first order for strong Bx and low temperatures. Hence, the PM to FM zero-temperature transition in Dy(OH)3 may be first order and not quantum critical. We investigate the effect of competing antiferromagnetic nearest-neighbor exchange and applied magnetic field, Bz, along the Ising spin direction ẑ on the first-order transition in Dy(OH)3. We conclude from these preliminary calculations that Ho(OH)3 and Dy(OH)3 and their Y3+ diamagnetically diluted variants, HoxY1−x(OH)3 and DyxY1−x(OH)3, are potentially interesting systems to study transverse-field-induced quantum fluctuations effects in hard axis (Ising-type) magnetic materials.

Photo-induced phase transitions in azobenzene-doped liquid crystals

The effect of irradiation (UV-visible light) on a nematic liquid crystal doped with a photoactive azobenzene derivative was investigated. The selective irradiation results in either an E implies Z or Z implies E isomerization of the azobenzene unit. The effect of the isomerization is to cause a reversible depression of the liquid crystal to isotropic (LC implies l) phase transition temperature of the doped mixture, which can be monitored optically as an isothermal phase transition. This depression also results in a biphasic liquid crystal+isotropic region which is discussed. The authors investigate the cause and magnitude of the phase depression as a function of the amount of doped 4-butyl-4'-methoxyazobenzene (photoactive unit) in 4-cyano-4'-n-pentylbiphenyl (liquid crystal unit), and as a function of the percentage conversion of E implies Z (caused by isomerization) in the azobenzene. The photostationary state of the doped mixtures achieved by Z implies E isomerization is considered and its effect upon the transition temperature of the mixture and response time of the system is discussed. They discuss the implications of the photostationary state with regards to the reversibility of the photo-induced phase transition and hence potential applications.

Spin-canted antiferromagnetic phase transitions in alternating phenoxo- and carboxylato-bridged Mn(III)-salen complexes

Three new MnIII complexes, {[Mn-2(salen)(2)(OCn)](ClO4)}(n) (1), {[Mn-2(salen)(2)(OPh)](ClO4)}(n) (2) and {[Mn-2(salen)(2)(OBz)](ClO4)}(2) (3) (where salen = N,N'-bis(salicylidene)-1,2-diaminoethane dianion, OCn = cinnamate, OPh = phenylacetate and OBz = benzoate), have been synthesized and characterized structurally and magnetically. The crystal structures reveal that all three structures contain syn-anti carboxylatebridged dimeric [Mn-2(salen)(2)(OOCR)](+) cations (OOCR = bridging carboxylate) that are joined together by weak Mn center dot center dot center dot O(phenoxo) interactions to form infinite alternating chain structures in 1 and 2, but the relatively long Mn center dot center dot center dot O(phenoxo) distance [3.621(2)angstrom] in 3 restricts this structure to tetranuclear units. Magnetic studies showed that 1 and 2 exhibited magnetic long-range order at T-N = 4.0 and 4.6 K (T-N = Neel transition temperature), respectively, to give spin-canted antiferromagnetic structures. Antiferromagnetic coupling was also observed in 3 but no peaks were recorded in the field-cooled magnetization (FCM) or zero-field-cooled magnetization (ZFCM) data, indicating that 3 remained paramagnetic down to 2 K. This dominant antiferromagnetic coupling is attributed to the carboxylate bridges. The ferromagnetic coupling expected due to the Mn-O(phenoxo)center dot center dot center dot Mn bridge plays an auxiliary role in the magnetic chain, but is an essential component of the bulk magnetic properties of the material.

SANS/WANS time-resolving neutron scattering studies of polymer phase transitions using NIMROD

We use new neutron scattering instrumentation to follow in a single quantitative time-resolving experiment, the three key scales of structural development which accompany the crystallisation of synthetic polymers. These length scales span 3 orders of magnitude of the scattering vector. The study of polymer crystallisation dates back to the pioneering experiments of Keller and others who discovered the chain-folded nature of the thin lamellae crystals which are normally found in synthetic polymers. The inherent connectivity of polymers makes their crystallisation a multiscale transformation. Much understanding has developed over the intervening fifty years but the process has remained something of a mystery. There are three key length scales. The chain folded lamellar thickness is ~ 10nm, the crystal unit cell is ~ 1nm and the detail of the chain conformation is ~ 0.1nm. In previous work these length scales have been addressed using different instrumention or were coupled using compromised geometries. More recently researchers have attempted to exploit coupled time-resolved small-angle and wide-angle x-ray experiments. These turned out to be challenging experiments much related to the challenge of placing the scattering intensity on an absolute scale. However, they did stimulate the possibility of new phenomena in the very early stages of crystallisation. Although there is now considerable doubt on such experiments, they drew attention to the basic question as to the process of crystallisation in long chain molecules. We have used NIMROD on the second target station at ISIS to follow all three length scales in a time-resolving manner for poly(e-caprolactone). The technique can provide a single set of data from 0.01 to 100Å-1 on the same vertical scale. We present the results using a multiple scale model of the crystallisation process in polymers to analyse the results.

Response of the northern stratospheric polar vortex to the seasonal alignment of QBO phase transitions

This study considers the strength of the Northern Hemisphere Holton-Tan effect (HTE) in terms of the phase alignment of the quasi-biennial oscillation (QBO) with respect to the annual cycle. Using the ERA-40 Reanalysis, it is found that the early winter (Nov–Dec) and late winter (Feb–Mar) relation between QBO phase and the strength of the stratospheric polar vortex is optimized for subsets of the 44-year record that are chosen on the basis of the seasonality of QBO phase transitions at the 30 hPa level. The timing of phase transitions serves as a proxy for changes in the vertical structure of the QBO over the whole depth of the tropical stratosphere. The statistical significance of the Nov–Dec (Feb–Mar) HTE is greatest when 30 hPa QBO phase transitions occur 9–14 (4–9) months prior to the January of the NH winter in question. This suggests that there exists for both early and late winter a vertical structure of tropical stratospheric winds that is most effective at influencing the interannual variability of the polar vortex, and that an early (late) winter HTE is associated with an early (late) progression of QBO phase towards that structure. It is also shown that the seasonality of QBO phase transitions at 30 hPa varies on a decadal timescale, with transitions during the first half of the calendar year being relatively more common during the first half of the tropical radiosonde wind record. Combining these two results suggests that decadal changes in HTE strength could result from the changing seasonality of QBO phase transitions. Citation: Anstey, J. A., and T. G. Shepherd (2008), Response of the northern stratospheric polar vortex to the seasonal alignment of QBO phase transitions, Geophys. Res. Lett., 35, L22810, doi:10.1029/2008GL035721.

Phase transitions of frozen camu-camu (Myrciaria dubia (HBK) McVaugh) pulp: Effect of cryostabilizer addition

Camu-camu is a tropical fruit with very high vitamin C content and commercialized as frozen pulp. Enthalpies of freezing, temperatures of the onset of ice melting, and glass transition temperatures of the maximally freeze-concentrated phase (T`(g)) of camu-camu pulp and of samples containing maltodextrin (DE20) and sucrose were measured by differential scanning calorimetry. Maltodextrin exhibited the largest freeze stabilization potential, increasing T`(g) from -58.2 degrees C (natural pulp) to -39.6 degrees C when 30% (w/w) maltodextrin DE 20 was added. Sucrose showed negligible effect on T`(g) but enhanced considerably the freezing point depression and less amount of ice was formed.

Phase transitions of cassava starch dispersions prepared with glycerol solutions

The aim of this work was to study the glass transition, the glass transition of the maximally freeze-concentrated fractions, the ice melting and the gelatinization phenomenon in dispersions of starch prepared using glycerol- water solutions. The starch concentration was maintained constant at 50 g cassava starch/100 g starch dispersions, but the concentration of the glycerol solutions was variable (C-g= 20, 40, 60, 80 and 100 mass/mass%). The phase transitions of these dispersions were studied by calorimetric methods, using a conventional differential scanning calorimeter (DSC) and a more sensitive equipment (micro-DSC). Apparently, in the glycerol diluted solutions (20 and 40%), the glycerol molecules interacted strongly with the glucose molecules of starch. While in the more concentrated glycerol domains (C-g> 40%), the behaviour was controlled by migration of water molecules from the starch granules, due to a hypertonic character of glycerol, which affected all phase transitions.

Bootstrap percolation on homogeneous trees has 2 phase transitions

We study the threshold theta bootstrap percolation model on the homogeneous tree with degree b + 1, 2 <= theta <= b, and initial density p. It is known that there exists a nontrivial critical value for p, which we call p(f), such that a) for p > p(f), the final bootstrapped configuration is fully occupied for almost every initial configuration, and b) if p < p(f) , then for almost every initial configuration, the final bootstrapped configuration has density of occupied vertices less than 1. In this paper, we establish the existence of a distinct critical value for p, p(c), such that 0 < p(c) < p(f), with the following properties: 1) if p <= p(c), then for almost every initial configuration there is no infinite cluster of occupied vertices in the final bootstrapped configuration; 2) if p > p(c), then for almost every initial configuration there are infinite clusters of occupied vertices in the final bootstrapped configuration. Moreover, we show that 3) for p < p(c), the distribution of the occupied cluster size in the final bootstrapped configuration has an exponential tail; 4) at p = p(c), the expected occupied cluster size in the final bootstrapped configuration is infinite; 5) the probability of percolation of occupied vertices in the final bootstrapped configuration is continuous on [0, p(f)] and analytic on (p(c), p(f) ), admitting an analytic continuation from the right at p (c) and, only in the case theta = b, also from the left at p(f).

Spectroscopic study of the effects of pressure media on high-pressure phase transitions in natrolite

Structural phase transitions in natrolite have been investigated as a function of pressure and different hydrostatic media using micro-Raman scattering and synchrotron infrared (IR) spectroscopy. Natrolite undergoes two reversible phase transitions at 0.86 and 1.53 GPa under pure water pressure medium. These phase transitions are characterized by the changes in the vibrational frequencies of four- and eight-membered rings related to the variations in the bridging T−O−T angles and the geometry of the elliptical eight-ring channels under pressure. Concomitant to the changes in the framework vibrational modes, the number of the O−H stretching vibrational modes of natrolite changes as a result of the rearrangements of the hydrogen bonds in the channels caused by a successive increase in the hydration level under hydrostatic pressure. Similar phase transitions were also observed at relatively higher pressures (1.13 and 1.59 GPa) under alcohol−water pressure medium. Furthermore, no phase transition was found up to 2.52 GPa if a lower volume ratio of the alcohol−water to natrolite was employed. This indicates that the water content in the pressure media plays a crucial role in triggering the pressure-induced phase transitions in natrolite. In addition, the average of the mode Grüneisen parameters is calculated to be about 0.6, while the thermodynamic Grüneisen parameter is found to be 1.33. This might be attributed to the contrast in the rigidity between the TO4 tetrahedral primary building units and other flexible secondary building units in the natrolite framework upon compression and subsequent water insertion.

New insights into the thermal behaviour of organic ionic plastic crystals: magnetic resonance imaging of polycrystalline morphology alterations induced by solid-solid phase transitions

Organic ionic plastic crystals (OIPCs) show strong potential as solid-state electrolytes for lithium battery applications, demonstrating promising electrochemical performance and eliminating the need for a volatile and flammable liquid electrolyte. The ionic conductivity (σ) in these systems has recently been shown to depend strongly on polycrystalline morphology, which is largely determined by the sample's thermal history. [K. Romanenko et al., J. Am. Chem. Soc., 2014, 136, 15638]. Tailoring this morphology could lead to conductivities sufficiently high for battery applications, so a more complete understanding of how phenomena such as solid-solid phase transitions can affect the sample morphology is of significant interest. Anisotropic relaxation of nuclear spin magnetisation provides a new MRI based approach for studies of polycrystalline materials at both a macroscopic and molecular level. In this contribution, morphology alterations induced by solid-solid phase transitions in triisobutyl(methyl)phosphonium bis(fluorosulfonyl)imide (P1444FSI) and diethyl(methyl)(isobutyl)phosphonium hexafluorophosphate (P1224PF6) are examined using magnetic resonance imaging (MRI), alongside nuclear magnetic resonance (NMR) spectroscopy, diffusion measurements and conductivity data. These observations are linked to molecular dynamics and structural behaviour crucial for the conductive properties of OIPCs. A distinct correlation is established between the conductivity at a given temperature, σ(T), and the intensity of the narrow NMR signal that is attributed to a mobile fraction, fm(T), of ions in the OIPC. To explain these findings we propose an analogy with the well-studied relationship between permeability (k) and void fraction (θ) in porous media, with k(θ) commonly quantified by a power-law dependence that can also be employed to describe σ(fm).

Structural phase transitions of PbZr0.52Ti0.48O3 ceramic: An infrared spectroscopy study

The recent discovery of a ferroelectric monoclinic phase in the PbZr1-xTixO3 (PZT) system attained the attention of several researchers due to the possibility of understanding the relationships between structural features and piezoelectric properties. The nature of the monoclinic phase in some PZT compositions remains controversial and unclear. In this work, structural phase transitions of PbZr0.52Ti0.48O3 ceramic were investigated by infrared spectroscopy as a function of temperature. Studies were centered on nu(1)-stretching modes and corresponding half width Wi as a function of temperature. The occurrence of the anomalies in the infrared spectra as a function of temperature suggests the following monoclinic ( LT) -> monoclinic ( HT) -> tetragonal phase transition were observed at 183 K and at 263 K.