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.

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.

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.

Out-of-equilibrium thermodynamics of quantum optomechanical systems

We address the out-of-equilibrium thermodynamics of an isolated quantum system consisting of a cavity optomechanical device. We explore the dynamical response of the system when driven out of equilibrium by a sudden quench of the coupling parameter and compute analytically the full distribution of the work generated by the process. We consider linear and quadratic optomechanical coupling, where the cavity field is parametrically coupled to either the position or the square of the position of a mechanical oscillator, respectively. In the former case we find that the average work generated by the quench is zero, whilst the latter leads to a non-zero average value. Through fluctuations theorems we access the most relevant thermodynamical figures of merit, such as the free energy difference and the amount of irreversible work generated. We thus provide a full charac- terization of the out-of-equilibrium thermodynamics in the quantum regime for nonlinearly coupled bosonic modes. Our study is the first due step towards the construction and full quantum analysis of an optomechanical machine working fully out of equilibrium.

Some current problems in perovskite nano-ferroelectrics and multiferroics: kinetically-limited systems of finite lateral size

We describe some unsolved problems of current interest; these involve quantum critical points in
ferroelectrics and problems which are not amenable to the usual density functional theory, nor to
classical Landau free energy approaches (they are kinetically limited), nor even to the Landau–
Kittel relationship for domain size (they do not satisfy the assumption of infinite lateral diameter)
because they are dominated by finite aperiodic boundary conditions.

Nonequilibrium Quantum Landauer Principle

Using the operational framework of completely positive, trace preserving operations and thermodynamic fluctuation relations, we derive a lower bound for the heat exchange in a Landauer erasure process on a quantum system. Our bound comes from a nonphenomenological derivation of the Landauer principle which holds for generic nonequilibrium dynamics. Furthermore, the bound depends on the nonunitality of dynamics, giving it a physical significance that differs from other derivations. We apply our framework to the model of a spin-1/2 system coupled to an interacting spin chain at finite temperature.