Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Concertedness and solvent effects in multiple proton transfer reactions: The formic acid dimer in solution

The issue of multiple proton transfer (PT) reactions in solution is addressed by performing molecular dynamics simulations for a formic acid dimer embedded in a water cluster. The reactant species is treated quantum mechanically, within a density functional approach, while the solvent is represented by a classical model. By constraining different distances within the dimer we analyze the PT process in a variety of situations representative of more complex environments. Free energy profiles are presented, and analyzed in terms of typical solvated configurations extracted from the simulations. A decrease in the PT barrier height upon solvation is rationalized in terms of a transition state which is more polarized than the stable states. The dynamics of the double PT process is studied in a low-barrier case and correlated with solvent polarization fluctuations. Cooperative effects in the motion of the two protons are observed in two different situations: when the solvent polarization does not favor the transfer of one of the two protons and when the motion of the two protons is not synchronized. This body of observations is correlated with local structural and dynamical properties of the solvent in the vicinity of the reactant. (C) 2000 American Institute of Physics. [S0021-9606(00)51121-0].

Solvation Structure and Transport of Acidic Protons in Ionic Liquids: A First-principles Simulation Study

Ab initio simulations of a single molecule of HCl in liquid dimethyl imidazolium chloride [dmim][Cl] show that the acidic proton exists as a symmetric, linear ClHCl- species. Details of the solvation structure around this molecule are given. The proton-transfer process was investigated by applying a force along the antisymmetric stretch coordinate until the molecule broke. Changes in the free energy and local solvation structure during this process were investigated. In the reaction mechanism identified, a free chloride approaches the proton from the side. As the original ClHCl- distorts and the incoming chloride forms a new bond to the proton, one of the original chlorine atoms is expelled and a new linear molecule is formed.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Thermodynamic Studies of Ionic Interactions in Aqueous Solutions of Imidazolium-Based Ionic Liquids [Emim][Br] and [Bmim][Cl]

Experimental measurements of density at different temperatures ranging from 293.15 to 313.15 K, the speed of sound and osmotic coefficients at 298.15 K for aqueous solution of 1-ethyl-3-methylimidazolium bromide ([Emim][Br]), and osmotic coefficients at 298.15 K for aqueous solutions of 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) in the dilute concentration region are taken. The data are used to obtain compressibilities, expansivity, apparent and limiting molar properties, internal pressure, activity, and activity coefficients for [Emim][Br] in aqueous solutions. Experimental activity coefficient data are compared with that obtained from Debye-Hückel and Pitzer models. The activity data are further used to obtain the hydration number and the osmotic second virial coefficients of ionic liquids. Partial molar entropies of [Bmim][Cl] are also obtained using the free-energy and enthalpy data. The distance of the closest approach of ions is estimated using the activity data for ILs in aqueous solutions and is compared with that of X-ray data analysis in the solid phase. The measured data show that the concentration dependence for aqueous solutions of [Emim][Br] can be accounted for in terms of the hydrophobic hydration of ions and that this IL exhibits Coulombic interactions as well as hydrophobic hydration for both the cations and anions. The small hydration numbers for the studied ILs indicate that the low charge density of cations and their hydrophobic nature is responsible for the formation of the water-structure-enforced ion pairs.

Effect of solution pH, ionic strength, and temperature on adsorption behavior of reactive dyes on activated carbon

The adsorption behavior of C.I. Reactive Blue 2, C.I. Reactive Red 4, and C.I. Reactive Yellow 2 from aqueous solution onto activated carbon was investigated under various experimental conditions. The adsorption capacity of activated carbon for reactive dyes was found to be relatively high. At pH 7.0 and 298 K, the maximum adsorption capacity for C.I. Reactive Blue 2, C.I. Reactive Yellow 2 and C.I. Reactive Red 4 dyes was found to be 0.27, 0.24, and 0.11 mmol/g, respectively. The shape of the adsorption isotherms indicated an L2-type isotherm according to the Giles and Smith classification. The experimental adsorption data showed good correlation with the Langmuir and Ferundlich isotherm models. Further analysis indicated that the formation of a complete monolayer was not achieved, with the fraction of surface coverage found to be 0.45, 0.42, and 0.22 for C.I. Reactive Blue 2, C.I. Reactive Yellow 2 and C.I. Reactive Red 4 dyes, respectively. Experimental data indicated that the adsorption capacity of activated carbon for the dyes was higher in acidic rather than in basic solutions, and further indicated that the removal of dye increased with increase in the ionic strength of solution, this was attributed to aggregation of reactive dyes in solution. Thermodynamic studies indicated that the adsorption of reactive dyes onto activated carbon was an endothermic process. The adsorption enthalpy (?H) for C.I. Reactive Blue 2 and C.I. Reactive Yellow 2 dyes were calculated at 42.2 and 36.2 kJ/mol, respectively. The negative values of free energy (?G) determined for these systems indicated that adsorption of reactive dyes was spontaneous at the temperatures under investigation (298-328 K). © 2007 Elsevier Ltd. All rights reserved.

Scaling of Superdomain Bands in Ferroelectric Dots

Bundles of 90° stripe domains have been observed to form into distinct groups, or bands, in mesoscale BaTiO3 single crystal dots. Vector piezoresponse force microscopy (PFM) shows that each band region, when considered as a single entity, possesses a resolved polarization that lies approximately along the pseudocubic direction; antiparallel alignment of this resultant polarization in adjacent bands means that these regions can be considered as 180° “superdomains.” For dots with sidewall dimensions below ~2 microns, Landau–Kittel like scaling in the width of these superdomains was observed, strongly suggesting that they form in response to lateral depolarizing fields. In larger dot structures, scaling laws break down. We have rationalized these observations by considering changes in the driving force for the adoption of equilibrium superdomain periodicities implied by Landau–Kittel-free energy models; we conclude that the formation of ordered bands of superdomains is a uniquely meso/nanoscale phenomenon. We also note that the superdomain bands found by PFM imaging in air contrast with the quadrant arrangements seen previously by Schilling et al. (Nano Lett., 9, 3359 (2009)) through transmission electron microscopy imaging in vacuum. The importance of the exact nature of the boundary conditions in determining the domain patterns that spontaneously form in nanostructures is therefore clearly implied.

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers.

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150–200 nucleotides at the 3' end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.

Volumetric properties and enthalpies of solution of alcohols CkH2k+1OH (k = 1, 2, 6) in 1-methyl-3-alkylimidazolium bis(trifluoromethylsulfonyl)imide {[C1CnIm][NTf2] n = 2, 4, 6, 8, 10} ionic liquids

We present a study on the effect of the alkyl chain length of the imidazolium ring in 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquids, [C1CnIm][NTf2] (n = 2 to 10), on the mixing properties of (ionic liquid + alcohol) mixtures (enthalpy and volume). We have measured small excess molar volumes with highly asymmetric curves as a function of mole fraction composition (S-shape) with more negative values in the alcohol-rich regions. The excess molar volumes increase with the increase of the alkyl-chain length of the imidazolium cation of the ionic liquid. The values of the partial molar excess enthalpy and the enthalpy of mixing are positive and, for the case of methanol, do not vary monotonously with the length of the alkyl side-chain of the cation on the ionic liquid – increasing from n = 2 to 6 and then decreasing from n = 8. This non-monotonous variation is explained by a more favourable interaction of methanol with the cation head group of the ionic liquid for alkyl chains longer than eight carbon atoms. It is also observed that the mixing is less favourable for the smaller alcohols, the enthalpy of mixing decreasing to less positive values as the alkyl chain of the alcohol increases. Based on the data from this work and on the knowledge of the vapour pressure of {[C1CnIm][NTf2] + alcohol} binary mixtures at T = 298 K reported in the literature, the excess Gibbs free energy, excess enthalpy and excess entropy could be then calculated and it was observed that these mixtures behave like the ones constituted by a non-associating and a non-polar component, with its solution behaviour being determined by the enthalpy.

Understanding the Performance of Melt-Extruded Poly(ethylene oxide)–Bicalutamide Solid Dispersions: Characterisation of Microstructural Properties Using Thermal, Spectroscopic and Drug Release Methods

In this article, we have prepared hot-melt-extruded solid dispersions of bicalutamide (BL) using poly(ethylene oxide) (PEO) as a matrix platform. Prior to preparation, miscibility of PEO and BL was assessed using differential scanning calorimetry (DSC). The onset of BL melting was signi?cantly depressed in the presence of PEO, and using Flory– Huggins (FH) theory, we identi?ed a negative value of -3.4, con?rming miscibility. Additionally, using FH lattice theory, we estimated the Gibbs free energy of mixing which was shown to be negative, passing through a minimum at a polymer fraction of 0.55. Using these data, solid dispersions at drug-to-polymer ratios of 1:10, 2:10 and 3:10 were prepared via hot-melt extrusion. Using a combination of DSC, powder X-ray diffractometry and scanning electron
microscopy, amorphous dispersions of BL were con?rmed at the lower two drug loadings. At the 3:10 BL to PEO ratio, crystalline BL was detected. The percent crystallinity of PEO was reduced by approximately 10% in all formulations following extrusion. The increased amorphous content within PEO following extrusion accommodated amorphous BL at drug to polymer loadings up to 2:10; however, the increased amorphous domains with PEO following extrusion were not suf?cient to fully accommodate BL at drug-to-polymer ratios of 3:10.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).

Domain annihilation due to temperature and thickness gradients in single-crystal BaTiO3

The manner in which 90? ferroelectric-ferroelastic domains respond to changes in temperature has been mapped in BaTiO3 single crystals using atomic force microscopy. Domain periodicity remains unaltered until approximately 2 ? C below the Curie temperature (TC ), whereupon domains coarsened dramatically. This behavior was successfully rationalized by considering the temperature dependence of the parameters associated with standard models of ferroelastic domain formation. However, while successful in describing the expected radical increase in equilibrium period with temperature, the model did not predict the unusual mechanism by which domain coarsening occurred; this was not continuous at a local level but instead involved discrete domain annihilation events. Subsequent insights from a combination of free energy analysis for the system and further experimental data from an analogous situation, in which domain period increases with increasing crystal thickness, suggested that domain annihilation is inevitable whenever a component of the relevant gradient that affects domain period is orientated parallel to the domain walls. Consistent with this thesis, we note that, for the observations presented herein, the thermal gradient possessed a significant component parallel to the domain walls. We suggest that domain annihilation is a general feature of domain structures in gradient fields.