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.

Interfacial Properties of LiTFSI and LiPF6-Based Electrolytes in Binary and Ternary Mixtures of Alkylcarbonates on Graphite Electrodes and Celgard Separator

For a better understanding of the adsorption behavior of alkylcarbonate-based electrolytes on graphite electrodes and Celgard separator for Li-ion batteries applications, the interface parameters are determined by contact angle and surface tension measurements. The correlation between these parameters and chemical compositions made of alkyl carbonate with a varying nature of lithium salts (LiPF6 and LiTFSI) and volume fractions of binary and ternary mixtures containing propylene carbonate (PC), ethylene carbonate (EC), and dimethyl carbonate (DMC) is investigated. From the obtained contact angle and surface tension (?L) values for each liquid, the dispersive and polar components of the surface tension (?Ld and ?Lp) of the electrolyte and interfacial free energy between the solid and liquid (?SL) were then calculated using the Young’s equation. The variation of contact angle (?) and the surface tension, as well as the work of adhesion (WA) of binary PC/DMC mixtures on PP, PE, and PET model surfaces were also measured and commented as function of volume fraction of PC in DMC. Finally, the Zisman’s critical surface tension (?C) for studied surfaces was then obtained showing positives slopes of cos ? versus ?L. This behavior is explained by a relative higher adsorption of alkylcarbonates to the hydrogenated supports or graphite. These results are decisive to understand the performance of electrolyte/electrode material/separator interfaces in lithium-ion battery devices.

Comparison of Dynamics of Extracellular Accesses to the β(1) and β(2) Adrenoceptors Binding Sites Uncovers the Potential of Kinetic Basis of Antagonist Selectivity.

From the molecular mechanism of antagonist unbinding in the ß(1) and ß(2) adrenoceptors investigated by steered molecular dynamics, we attempt to provide further possibilities of ligand subtype and subspecies selectivity. We have simulated unbinding of ß(1) -selective Esmolol and ß(2) -selective ICI-118551 from both receptors to the extracellular environment and found distinct molecular features of unbinding. By calculating work profiles, we show different preference in antagonist unbinding pathways between the receptors, in particular, perpendicular to the membrane pathway is favourable in the ß(1) adrenoceptor, whereas the lateral pathway involving helices 5, 6 and 7 is preferable in the ß(2) adrenoceptor. The estimated free energy change of unbinding based on the preferable pathway correlates with the experimental ligand selectivity. We then show that the non-conserved K347 (6.58) appears to facilitate in guiding Esmolol to the extracellular surface via hydrogen bonds in the ß(1) adrenoceptor. In contrast, hydrophobic and aromatic interactions dominate in driving ICI-118551 through the easiest pathway in the ß(2) adrenoceptor. We show how our study can stimulate design of selective antagonists and discuss other possible molecular reasons of ligand selectivity, involving sequential binding of agonists and glycosylation of the receptor extracellular surface. © 2012 John Wiley & Sons A/S.

Excess electron interactions with solvated DNA nucleotides:Strand breaks possible at room temperature

When biological matter is subjected to ionizing radiation, a wealth of secondary low-energy (<20 eV) electrons are produced. These electrons propagate inelastically, losing energy to the medium until they reach energies low enough to localize in regions of high electron affinity. We have recently shown that in fully solvated DNA fragments, nucleobases are particularly attractive for such excess electrons. The next question is what is their longer-term effect on DNA. It has been advocated that they can lead to strand breaks by cleavage of the phosphodiester C-3'-O-3' bond. Here we present a first-principles study of free energy barriers for the cleavage of this bond in fully solvated nucleotides. We have found that except for dAMP, the barriers are on the order of 6 kcal/mol, suggesting that bond cleavage is a regular feature at 300 K. Such low barriers are possible only as a result of solvent and thermal fluctuations. These findings support the notion that low-energy electrons can indeed lead to strand breaks in DNA.

Construction of Drug–Polymer Thermodynamic Phase Diagrams Using Flory–Huggins Interaction Theory: Identifying the Relevance of Temperature and Drug Weight Fraction to Phase Separation within Solid Dispersions

Amorphous drug-polymer solid dispersions have the potential to enhance the dissolution performance and thus bioavailability of BCS class II drug compounds. The principle drawback of this approach is the limited physical stability of amorphous drug within the dispersion. Accurate determination of the solubility and miscibility of drug in the polymer matrix is the key to the successful design and development of such systems. In this paper, we propose a novel method, based on Flory-Huggins theory, to predict and compare the solubility and miscibility of drug in polymeric systems. The systems chosen for this study are (1) hydroxypropyl methylcellulose acetate succinate HF grade (HPMCAS-HF)-felodipine (FD) and (2) Soluplus (a graft copolymer of polyvinyl caprolactam-polyvinyl acetate-polyethylene glycol)-FD. Samples containing different drug compositions were mixed, ball milled, and then analyzed by differential scanning calorimetry (DSC). The value of the drug-polymer interaction parameter ? was calculated from the crystalline drug melting depression data and extrapolated to lower temperatures. The interaction parameter ? was also calculated at 25 °C for both systems using the van Krevelen solubility parameter method. The rank order of interaction parameters of the two systems obtained at this temperature was comparable. Diagrams of drug-polymer temperature-composition and free energy of mixing (?G mix) were constructed for both systems. The maximum crystalline drug solubility and amorphous drug miscibility may be predicted based on the phase diagrams. Hyper-DSC was used to assess the validity of constructed phase diagrams by annealing solid dispersions at specific drug loadings. Three different samples for each polymer were selected to represent different regions within the phase diagram

Using sketch-map coordinates to analyze and bias molecular dynamics simulations

When examining complex problems, such as the folding of proteins, coarse grained descriptions of the system drive our investigation and help us to rationalize the results. Oftentimes collective variables (CVs), derived through some chemical intuition about the process of interest, serve this purpose. Because finding these CVs is the most difficult part of any investigation, we recently developed a dimensionality reduction algorithm, sketch-map, that can be used to build a low-dimensional map of a phase space of high-dimensionality. In this paper we discuss how these machine-generated CVs can be used to accelerate the exploration of phase space and to reconstruct free-energy landscapes. To do so, we develop a formalism in which high-dimensional configurations are no longer represented by low-dimensional position vectors. Instead, for each configuration we calculate a probability distribution, which has a domain that encompasses the entirety of the low-dimensional space. To construct a biasing potential, we exploit an analogy with metadynamics and use the trajectory to adaptively construct a repulsive, history-dependent bias from the distributions that correspond to the previously visited configurations. This potential forces the system to explore more of phase space by making it desirable to adopt configurations whose distributions do not overlap with the bias. We apply this algorithm to a small model protein and succeed in reproducing the free-energy surface that we obtain from a parallel tempering calculation.

Locating binding poses in protein-ligand systems using reconnaissance metadynamics

A molecular dynamics-based protocol is proposed for finding and scoring protein-ligand binding poses. This protocol uses the recently developed reconnaissance metadynamics method, which employs a self-learning algorithm to construct a bias that pushes the system away from the kinetic traps where it would otherwise remain. The exploration of phase space with this algorithm is shown to be roughly six to eight times faster than unbiased molecular dynamics and is only limited by the time taken to diffuse about the surface of the protein. We apply this method to the well-studied trypsin-benzamidine system and show that we are able to refind all the poses obtained from a reference EADock blind docking calculation. These poses can be scored based on the length of time the system remains trapped in the pose. Alternatively, one can perform dimensionality reduction on the output trajectory and obtain a map of phase space that can be used in more expensive free-energy calculations.

The phase diagram of water at negative pressures: Virtual ices

The phase diagram of water at negative pressures as obtained from computer simulations for two models of water, TIP4P/2005 and TIP5P is presented. Several solid structures with lower densities than ice Ih, so-called virtual ices, were considered as possible candidates to occupy the negative pressure region of the phase diagram of water. In particular the empty hydrate structures sI, sII, and sH and another, recently proposed, low-density ice structure. The relative stabilities of these structures at 0 K was determined using empirical water potentials and density functional theory calculations. By performing free energy calculations and Gibbs-Duhem integration the phase diagram of TIP4P/2005 was determined at negative pressures. The empty hydrates sII and sH appear to be the stable solid phases of water at negative pressures. The phase boundary between ice Ih and sII clathrate occurs at moderate negative pressures, while at large negative pressures sH becomes the most stable phase. This behavior is in reasonable agreement with what is observed in density functional theory calculations.

Binding of Calcium and Carbonate to Polyacrylates

Polyacrylate molecules can be used to slow the growth of calcium carbonate. However, little is known about the mechanism by which the molecules impede the growth rate. A recent computational study (Bulo et al. Macromolecules 2007, 40, 3437) used metadynamics to investigate the binding of calcium to polyacrylate chains and has thrown some light on the coiling and precipitation of these polymers. We extend these simulations to examine the binding of calcium and carbonate to polyacrylate chains. We show that calcium complexed with both carbonate and polyacrylate is a very stable species. The free energies of calcium-carbonate-polyacrylate complexes, with different polymer configurations, are calculated, and differences in the free energy of the binding of carbonate are shown to be due to differences in the amount of steric hindrance about the calcium, which prevents the approach of the carbonate ion.

Measuring the Characteristic Function of the Work Distribution

We propose an interferometric setting for the ancilla-assisted measurement of the characteristic function of the work distribution following a time-dependent process experienced by a quantum system. We identify how the configuration of the effective interferometer is linked to the symmetries enjoyed by the Hamiltonian ruling the process and provide the explicit form of the operations to implement in order to accomplish our task. We finally discuss two physical settings, based on hybrid optomechanical-electromechanical devices, where the theoretical proposals discussed in our work could find an experimental demonstration.

Molecular dynamics simulation of nanoindentation of Fe3C and tetrahedral Fe4C

Study of nanomechanical response of iron carbides is important because presence of iron carbides greatly influences the performance and longevity of steel components. This work contributes to the literature by exploring nanoindentation of nanocrystalline Fe3C and tetrahedral-Fe4C using molecular dynamics simulation. The chemical interactions of iron and carbon were described through an analytical bond order inter-atomic potential (ABOP) energy function. The indentations were performed at an indentation speed of 50 m/sec and a repeat trial was performed at 5 m/sec. Load-displacement (P-h) curve for both these carbides showed residual indentation depth and maximum indentation depth (hf/hmax) ratio to be higher than 0.7 i.e. a circumstance where Oliver and Pharr method was not appropriate to be applied to evaluate the material properties. Alternate evaluation revealed Fe3C to be much harder than Fe4C. Gibbs free energy of formation and radial distribution function, coupled with state of the average local temperature and von Mises stresses indicate the formation of a new phase of iron-carbide. Formation of this newer phase was found to be due to deviatoric strain rather than the high temperature induced in the substrate during nanoindentation

PLUMED 2: New feathers for an old bird

Enhancing sampling and analyzing simulations are central issues in molecular simulation. Recently, we introduced PLUMED, an open-source plug-in that provides some of the most popular molecular dynamics (MD) codes with implementations of a variety of different enhanced sampling algorithms and collective variables (CVs). The rapid changes in this field, in particular new directions in enhanced sampling and dimensionality reduction together with new hardware, require a code that is more flexible and more efficient. We therefore present PLUMED 2 here a,complete rewrite of the code in an object-oriented programming language (C++). This new version introduces greater flexibility and greater modularity, which both extends its core capabilities and makes it far easier to add new methods and CVs. It also has a simpler interface with the MD engines and provides a single software library containing both tools and core facilities. Ultimately, the new code better serves the ever-growing community of users and contributors in coping with the new challenges arising in the field.

Program summary

Program title: PLUMED 2

Catalogue identifier: AEEE_v2_0

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEEE_v2_0.html

Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland

Licensing provisions: Yes

No. of lines in distributed program, including test data, etc.: 700646

No. of bytes in distributed program, including test data, etc.: 6618136

Distribution format: tar.gz

Programming language: ANSI-C++.

Computer: Any computer capable of running an executable produced by a C++ compiler.

Operating system: Linux operating system, Unix OSs.

Has the code been vectorized or parallelized?: Yes, parallelized using MPI.

RAM: Depends on the number of atoms, the method chosen and the collective variables used.

Classification: 3, 7.7, 23. Catalogue identifier of previous version: AEEE_v1_0.

Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 1961.

External routines: GNU libmatheval, Lapack, Bias, MPI. (C) 2013 Elsevier B.V. All rights reserved.

Detecting the work statistics through Ramsey-like interferometry

Out-of-equilibrium statistical mechanics is attracting considerable interest due to the recent advances in the control and manipulations of systems at the quantum level. Recently, an interferometric scheme for the detection of the characteristic function of the work distribution following a time-dependent process has been proposed [L. Mazzola et al., Phys. Rev. Lett. 110 (2013) 230602]. There, it was demonstrated that the work statistics of a quantum system undergoing a process can be reconstructed by effectively mapping the characteristic function of work on the state of an ancillary qubit. Here, we expand that work in two important directions. We first apply the protocol to an interesting specific physical example consisting of a superconducting qubit dispersively coupled to the field of a microwave resonator, thus enlarging the class of situations for which our scheme would be key in the task highlighted above. We then account for the interaction of the system with an additional one (which might embody an environment), and generalize the protocol accordingly.

Low pressure methane solubility in lithium-ion batteries based solvents and electrolytes as a function of temperature. Measurement and prediction

The methane solubility in five pure electrolyte solvents and one binary solvent mixture for lithium ion batteries – such as ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), ethyl methyl carbonate (EMC), diethyl carbonate (DEC) and the (50:50 wt%) mixture of EC:DMC was studied experimentally at pressures close to atmospheric and as a function of temperature between (280 and 343) K by using an isochoric saturation technique. The effect of the selected anions of a lithium salt LiX (X = hexafluorophosphate,