Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment.

ABSTRACT: The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.

SUBMITTER: Gelman D

PROVIDER: S-EPMC3041154 | biostudies-other | 2011 Jan

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

1.14*CM1A-LBCC: Localized Bond-Charge Corrected CM1A Charges for Condensed-Phase Simulations.

Project description:The quality of the 1.14*CM1A and 1.20*CM5 charge models was evaluated for calculations of free energies of hydration. For a set of 426 neutral molecules, 1.14*CM1A and 1.20*CM5 yield MADs of 1.26 and 1.21 kcal/mol, respectively. The 1.14*CM1A charges, which can be readily obtained for large systems, exhibit large deviations only for a subset of functional groups. The results for these cases were systematically improved using localized bond-charge corrections (LBCC) by which offsetting adjustments are made to the partial charges for atoms in specified bond types. Only 19 LBCCs were needed to yield 1.14*CM1A-LBCC charges that reduce the errors for the 426 ?Ghyd values to only 0.61 kcal/mol. The modified charge method was also tested in computation of heats of vaporization and densities for pure organic liquids, yielding average errors of 1.40 kcal/mol and 0.024 g/cm3, similar to those for 1.14*CM1A.

| S-EPMC5813481 | biostudies-literature

Phase-space methods for the spin dynamics in condensed matter systems.

Project description:Using the phase-space formulation of quantum mechanics, we derive a four-component Wigner equation for a system composed of spin-[Formula: see text] fermions (typically, electrons) including the Zeeman effect and the spin-orbit coupling. This Wigner equation is coupled to the appropriate Maxwell equations to form a self-consistent mean-field model. A set of semiclassical Vlasov equations with spin effects is obtained by expanding the full quantum model to first order in the Planck constant. The corresponding hydrodynamic equations are derived by taking velocity moments of the phase-space distribution function. A simple closure relation is proposed to obtain a closed set of hydrodynamic equations.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

| S-EPMC5360899 | biostudies-literature

Finite-Size-Corrected Rotational Diffusion Coefficients of Membrane Proteins and Carbon Nanotubes from Molecular Dynamics Simulations.

Project description:We investigate system-size effects on the rotational diffusion of membrane proteins and other membrane-embedded molecules in molecular dynamics simulations. We find that the rotational diffusion coefficient slows down relative to the infinite-system value by a factor of one minus the ratio of protein and box areas. This correction factor follows from the hydrodynamics of rotational flows under periodic boundary conditions and is rationalized in terms of Taylor-Couette flow. For membrane proteins like transporters, channels, or receptors in typical simulation setups, the protein-covered area tends to be relatively large, requiring a significant finite-size correction. Molecular dynamics simulations of the protein adenine nucleotide translocase (ANT1) and of a carbon nanotube porin in lipid membranes show that the hydrodynamic finite-size correction for rotational diffusion is accurate in standard-use cases. The dependence of the rotational diffusion on box size can be used to determine the membrane viscosity.

| S-EPMC6750896 | biostudies-literature

Condensed-Phase Molecular Representation to Link Structure and Thermodynamics in Molecular Dynamics.

Project description:Molecular design requires systematic and broadly applicable methods to extract structure-property relationships. The focus of this study is on learning thermodynamic properties from molecular-liquid simulations. The methodology relies on an atomic representation originally developed for electronic properties: the Spectrum of London and Axilrod-Teller-Muto representation (SLATM). SLATM's expansion in one-, two-, and three-body interactions makes it amenable to probing structural ordering in molecular liquids. We show that such representation encodes enough critical information to permit the learning of thermodynamic properties via linear methods. We demonstrate our approach on the preferential insertion of small solute molecules toward cardiolipin membranes and monitor selectivity against a similar lipid. Our analysis reveals simple, interpretable relationships between two- and three-body interactions and selectivity, identifies key interactions to build optimal prototypical solutes, and charts a two-dimensional projection that displays clearly separated basins. The methodology is generally applicable to a variety of thermodynamic properties.

| S-EPMC10373487 | biostudies-literature

(Meta-)stability and Core-Shell Dynamics of Gold Nanoclusters at Finite Temperature.

Project description:Gold nanoclusters have been the focus of numerous computational studies, but an atomistic understanding of their structural and dynamical properties at finite temperature is far from satisfactory. To address this deficiency, we investigate gold nanoclusters via ab initio molecular dynamics, in a range of sizes where a core-shell morphology is observed. We analyze their structure and dynamics using state-of-the-art techniques, including unsupervised machine-learning nonlinear dimensionality reduction (sketch-map) for describing the similarities and differences among the range of sampled configurations. In the examined temperature range between 300 and 600 K, we find that whereas the gold nanoclusters exhibit continuous structural rearrangement, they are not amorphous. Instead, they clearly show persistent motifs: a cationic core of 1-5 atoms is loosely bound to a shell which typically displays a substructure resulting from the competition between locally spherical versus planar fragments. Besides illuminating the properties of core-shell gold nanoclusters, the present study proposes a set of useful tools for understanding their nature in operando.

| S-EPMC6372262 | biostudies-literature

Double proton transfer dynamics of model DNA base pairs in the condensed phase.

Project description:The dynamics of excited-state double proton transfer of model DNA base pairs, 7-azaindole dimers, is reported using femtosecond fluorescence spectroscopy. To elucidate the nature of the transfer in the condensed phase, here we examine variation of solvent polarity and viscosity, solute concentration, and isotopic fractionation. The rate of proton transfer is found to be significantly dependent on polarity and on the isotopic composition in the pair. Consistent with a stepwise mechanism, the results support the presence of an ionic intermediate species which forms on the femtosecond time scale and decays to the final tautomeric form on the picosecond time scale. We discuss the results in relation to the molecular motions involved and comment on recent claims of concerted transfer in the condensed phase. The nonconcerted mechanism is in agreement with previous isolated-molecule femtosecond dynamics and is also consistent with the most-recent high-level theoretical study on the same pair.

| S-EPMC1885566 | biostudies-literature

Hybrid finite element and Brownian dynamics method for diffusion-controlled reactions.

Project description:Diffusion is often the rate determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. This paper proposes a new hybrid diffusion method that couples the strengths of each of these two methods. The method is derived for a general multidimensional system, and is presented using a basic test case for 1D linear and radially symmetric diffusion systems.

| S-EPMC3350536 | biostudies-other

Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method.

Project description:To relate the subcellular molecular events to organ level physiology in heart, we have developed a three-dimensional finite-element-based simulation program incorporating the cellular mechanisms of excitation-contraction coupling and its propagation, and simulated the fluid-structure interaction involved in the contraction and relaxation of the human left ventricle. The FitzHugh-Nagumo model and four-state model representing the cross-bridge kinetics were adopted for cellular model. Both ventricular wall and blood in the cavity were modeled by finite element mesh. An arbitrary Lagrangian Eulerian finite element method with automatic mesh updating has been formulated for large domain changes, and a strong coupling strategy has been taken. Using electrical analog of pulmonary circulation and left atrium as a preload and the windkessel model as an afterload, dynamics of ventricular filling as well as ejection was simulated. We successfully reproduced the biphasic filling flow consisting of early rapid filling and atrial contraction similar to that reported in clinical observation. Furthermore, fluid-structure analysis enabled us to analyze the wave propagation velocity of filling flow. This simulator can be a powerful tool for establishing a link between molecular abnormality and the clinical disorder at the macroscopic level.

| S-EPMC1304609 | biostudies-literature

Ab-Initio Molecular Dynamics Simulation of Condensed-Phase Reactivity: The Electrolysis of Amino Acids and Peptides.

Project description:Electrolysis is a potential candidate for a quick method of wastewater cleansing. However, it is necessary to know what compounds might be formed from bioorganic matter. We want to know if there are toxic intermediates and if it is possible to influence the product formation by the variation in initial conditions. In the present study, we use Car-Parrinello molecular dynamics to simulate the fastest reaction steps under such circumstances. We investigate the behavior of amino acids and peptides under anodic conditions. Such highly reactive situations lead to chemical reactions within picoseconds, and we can model the reaction mechanisms in full detail. The role of the electric current is to discharge charged species and, hence, to produce radicals from ions. This leads to ultra-fast radical reactions in a bulk environment, which can also be seen as redox reactions as the oxidation states change. In the case of amino acids, the educts can be zwitterionic, so we also observe complex acid-base chemistry. Hence, we obtain the full spectrum of condensed-phase chemistry.

| S-EPMC7699423 | biostudies-literature

Evidence for a Finite-Temperature Insulator.

Project description:In superconductors the zero-resistance current-flow is protected from dissipation at finite temperatures (T) by virtue of the short-circuit condition maintained by the electrons that remain in the condensed state. The recently suggested finite-T insulator and the "superinsulating" phase are different because any residual mechanism of conduction will eventually become dominant as the finite-T insulator sets-in. If the residual conduction is small it may be possible to observe the transition to these intriguing states. We show that the conductivity of the high magnetic-field insulator terminating superconductivity in amorphous indium-oxide exhibits an abrupt drop, and seem to approach a zero conductance at T < 0.04 K. We discuss our results in the light of theories that lead to a finite-T insulator.

| S-EPMC4550897 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data