1
|
Walsh MR. Comparing brute force to transition path sampling for gas hydrate nucleation with a flat interface: comments on time reversal symmetry. Phys Chem Chem Phys 2024; 26:5762-5772. [PMID: 38214888 DOI: 10.1039/d3cp05059a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Fluid to solid nucleation is often investigated with the rare event method transition path sampling (TPS). I claim that the inherent irreversibility of solid nucleation, even at stationary conditions, calls into question TPS's applicability for determining solid nucleation mechanisms, especially for pre-critical behavior. Even when applied to a phenomenon which displays time reversal asymmetry like solid nucleation, TPS is a good means of exploring phase space and giving trends in post-critical structure, and its ability to facilitate nucleation rate and free energy calculations remains outstanding. Forward-only splitting and ratcheting methods such as forward flux sampling are more attractive for understanding nucleation mechanisms as they do not require time reversal symmetry, but at low driving forces may suffer from the same limitations as brute force: they may never make it to the first ratchet. Here I briefly summarize the TPS method and gas hydrate nucleation simulation literature, focusing on topics within both to facilitate a comparison of brute force hydrate nucleation to transition path sampling of hydrate nucleation. Perhaps anecdotally, the brute force technique results in more crystalline trajectories despite having higher driving forces than TPS. I maintain this difference is because of the inherent irreversibility of hydrate nucleation, meaning its pre-critical behavior cannot accurately be determined by the melting trajectories that comprise approximately half of the configurations in TPS's path ensemble.
Collapse
|
2
|
Lin N, Mani T. Anti-Arrhenius behavior of electron transfer reactions in molecular dimers. Chem Sci 2023; 14:13095-13107. [PMID: 38023507 PMCID: PMC10664467 DOI: 10.1039/d3sc03609j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Accepted: 10/29/2023] [Indexed: 12/01/2023] Open
Abstract
Rates of chemical reactions typically accelerate as the temperature rises, following the Arrhenius law. However, electron transfer reactions may exhibit weak temperature dependence or counterintuitive behavior, known as anti-Arrhenius behavior, wherein reaction rates decrease as temperature increases. Solvent reorganization energy and torsion-induced changes in electronic couplings could contribute to this unusual behavior, but how each contributes to the overall temperature dependence is unclear. One can decelerate the charge recombination process in photogenerated radical pairs or charge-separated states by harnessing this often-overlooked phenomenon. This means that we could achieve long-lived radical pairs without relying on conventional cooling. Using a series of homo molecular dimers, we showed that the degree of torsional hindrance dictates temperature-dependent torsion-induced changes in electronic coupling and, therefore, charge recombination rates. The overall temperature dependence is controlled by how changes in electronic coupling and the temperature-dependent solvent reorganization energy contribute to the rates of charge recombination. Our findings pave the way for rationally designing molecules that exhibit anti-Arrhenius behavior to slow down charge recombination, opening possibilities for applications in energy-related and quantum information technologies.
Collapse
Affiliation(s)
- Neo Lin
- Department of Chemistry, University of Connecticut Storrs CT 06269 USA
| | - Tomoyasu Mani
- Department of Chemistry, University of Connecticut Storrs CT 06269 USA
- Chemistry Division, Brookhaven National Laboratory Upton NY 11973 USA
| |
Collapse
|
3
|
Naleem N, Abreu CRA, Warmuz K, Tong M, Kirmizialtin S, Tuckerman ME. An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions. J Chem Phys 2023; 159:034102. [PMID: 37458344 DOI: 10.1063/5.0147597] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 06/23/2023] [Indexed: 07/20/2023] Open
Abstract
Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.
Collapse
Affiliation(s)
- Nawavi Naleem
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
| | - Charlles R A Abreu
- Chemical Engineering Department, Escola de Química, Universidade Federal do Rio de Janeiro, 21941-909 Rio de Janeiro, RJ, Brazil
| | - Krzysztof Warmuz
- Computer Science Program, Science Division, New York University, Abu Dhabi, UAE
| | - Muchen Tong
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
| | - Serdal Kirmizialtin
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Center for Smart Engineering Materials, New York University, Abu Dhabi, UAE
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Rd. North, Shanghai 200062, China
- Simons Center for Computational Physical Chemistry at New York University, New York, New York 10003, USA
| |
Collapse
|
4
|
Shayesteh Zadeh A, Khan SA, Vandervelden C, Peters B. Site-Averaged Ab Initio Kinetics: Importance Learning for Multistep Reactions on Amorphous Supports. J Chem Theory Comput 2023; 19:2873-2886. [PMID: 37093705 DOI: 10.1021/acs.jctc.3c00160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Single-atom centers on amorphous supports include catalysts for polymerization, partial oxidation, metathesis, hydrogenolysis, and more. The disordered environment makes each site different, and the kinetics exponentially magnifies these differences to make ab initio site-averaged kinetics calculations extremely difficult. This work extends the importance learning algorithm for efficient and precise site-averaged kinetics estimates to ab initio calculations and multistep reaction mechanisms. Specifically, we calculate site-averaged proton transfer relaxation rates on an ensemble of cluster models representing Brønsted acid sites on silica-alumina. We include direct and water-assisted proton transfer pathways and simultaneously estimate the water adsorption and activation enthalpies for forward and backward proton transfers. We use density functional theory (DFT) to obtain a site-averaged rate, somewhat like a turnover frequency, for the proton transfer relaxation rate. Finally, we show that importance learning can provide orders-of-magnitude acceleration over standard sampling methods for site-averaged rate calculations in cases where the rate is dominated by a few highly active sites.
Collapse
Affiliation(s)
- Armin Shayesteh Zadeh
- Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Salman A Khan
- Delaware Energy Institute (DEI), University of Delaware, Newark, Delaware 19711, United States
| | | | - Baron Peters
- Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
5
|
Janjic P, Solev D, Kocarev L. Non-trivial dynamics in a model of glial membrane voltage driven by open potassium pores. Biophys J 2023; 122:1470-1490. [PMID: 36919241 PMCID: PMC10147837 DOI: 10.1016/j.bpj.2023.03.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2022] [Revised: 02/01/2023] [Accepted: 03/07/2023] [Indexed: 03/16/2023] Open
Abstract
Despite the molecular evidence that a nearly linear steady-state current-voltage relationship in mammalian astrocytes reflects a total current resulting from more than one differentially regulated K+ conductance, detailed ordinary differential equation (ODE) models of membrane voltage Vm are still lacking. Various experimental results reporting altered rectification of the major Kir currents in glia, dominated by Kir4.1, have motivated us to develop a detailed model of Vm dynamics incorporating the weaker potassium K2P-TREK1 current in addition to Kir4.1, and study the stability of the resting state Vr. The main question is whether, with the loss of monotonicity in glial I-V curve resulting from altered Kir rectification, the nominal resting state Vr remains stable, and the cell retains the trivial, potassium electrode behavior with Vm after EK. The minimal two-dimensional model of Vm near Vr showed that an N-shape deformed Kir I-V curve induces multistability of Vm in a model that incorporates K2P activation kinetics, and nonspecific K+ leak currents. More specifically, an asymmetrical, nonlinear decrease of outward Kir4.1 conductance, turning the channels into inward rectifiers, introduces instability of Vr. That happens through a robust bifurcation giving birth to a second, more depolarized stable resting state Vdr > -10 mV. Realistic recordings from electrographic seizures were used to perturb the model. Simulations of the model perturbed by constant current through gap junctions and seizure-like discharges as local field potentials led to depolarization and switching of Vm between the two stable states, in a downstate-upstate manner. In the event of prolonged depolarizations near Vdr, such catastrophic instability would affect all aspects of the glial function, from metabolic support to membrane transport, and practically all neuromodulatory roles assigned to glia.
Collapse
Affiliation(s)
- Predrag Janjic
- Laboratory for Complex Systems and Networks, Research Centre for Computer Science and Information Technologies, Macedonian Academy of Sciences and Arts, Skopje, North Macedonia.
| | - Dimitar Solev
- Laboratory for Complex Systems and Networks, Research Centre for Computer Science and Information Technologies, Macedonian Academy of Sciences and Arts, Skopje, North Macedonia
| | - Ljupco Kocarev
- Laboratory for Complex Systems and Networks, Research Centre for Computer Science and Information Technologies, Macedonian Academy of Sciences and Arts, Skopje, North Macedonia
| |
Collapse
|
6
|
Kamat K, Naullage PM, Molinero V, Peters B. Oriented attachment kinetics for rod-like particles at a flat surface: Buffon's needle at the nanoscale. J Chem Phys 2022; 157:214113. [PMID: 36511557 DOI: 10.1063/5.0124531] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
The adsorption of large rod-like molecules or crystallites on a flat crystal face, similar to Buffon's needle, requires the rods to "land," with their binding sites in precise orientational alignment with matching sites on the surface. An example is provided by long, helical antifreeze proteins (AFPs), which bind at specific facets and orientations on the ice surface. The alignment constraint for adsorption, in combination with the loss in orientational freedom as the molecule diffuses toward the surface, results in an entropic barrier that hinders the adsorption. Prior kinetic models do not factor in the complete geometry of the molecule, nor explicitly enforce orientational constraints for adsorption. Here, we develop a diffusion-controlled adsorption theory for AFP molecules binding at specific orientations to flat ice surfaces. We formulate the diffusion equation with relevant boundary conditions and present analytical solutions to the attachment rate constant. The resulting rate constant is a function of the length and aspect ratio of the AFP, the distance threshold associated with binding, and solvent conditions such as temperature and viscosity. These results and methods of calculation may also be useful for predicting the kinetics of crystal growth through oriented attachment.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Pavithra M Naullage
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Baron Peters
- Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
7
|
Cabanas AM, Flores Araya JC, Jessop IA, Humire F. Anomalous (Exergonic) Behavior in the Transfer of Electrons between Donors and Acceptors: Mobility, Energy, Caloric Capacity, and Entropy. ACS OMEGA 2022; 7:35153-35158. [PMID: 36211079 PMCID: PMC9535709 DOI: 10.1021/acsomega.2c04094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 09/08/2022] [Indexed: 06/16/2023]
Abstract
Understanding the kinetics of electron transfer reactions involves active research in physics, chemistry, biology, and nano-tech. Here, we propose a model to apply in a broader framework by establishing a connection between thermodynamics and kinetics. From a purely thermodynamic point of view, electronic transfer Marcus' theory is revisited; consequently, calculations of thermodynamic variables such as mobility, energy, and entropy are provided. More significantly, two different regimes are explicitly established. In the anomalous region, an exergonic process associated with negative heat capacity appears. Further, in the same region, mobility, energy, and entropy decrease when the temperature increases.
Collapse
Affiliation(s)
- Ana M. Cabanas
- Departamento
de Física, FACI, Universidad de Tarapacá, Arica 1000965, Chile
| | | | - Ignacio A. Jessop
- Departamento
de Química, FACI, Universidad de
Tarapacá, Arica 1000007, Chile
| | - Fernando Humire
- Departamento
de Física, FACI, Universidad de Tarapacá, Arica 1000965, Chile
| |
Collapse
|
8
|
Marinova V, Freeman CL, Harding JH. Significance of atomic-scale defects in flexible surfaces on local solvent and ion behaviour. Faraday Discuss 2022; 235:289-306. [PMID: 35380136 DOI: 10.1039/d1fd00082a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Many factors can affect the course of heterogeneous nucleation, such as surface chemistry, flexibility and topology, substrate concentration and solubility. Atomic-scale defects are rarely investigated in detail and are often considered to be unimportant surface features. In this work, we set out to investigate the significance of atomic-scale defects in a flexible self-assembled monolayer surface for the behaviour of clusters of Ca2+ and CO32- ions in water. To this end, we use molecular dynamics simulations to estimate the diffusion coefficients of ion clusters at different topological surface features and obtain ionic radial distribution functions around features of interest. Well-tempered metadynamics is used to gain insight into the free energy of ions around selected surface defects. We find that certain defects, which we refer to as active defects, can impair ionic surface diffusion, as well as affect the diffusion of ions in close proximity to the surface feature in question. Our findings suggest that this effect can result in an ability of such topological features to promote ion clustering and increase local ionic concentration at specific surface sites. The work reported here shows how the presence of small atomic-scale defects can affect the role of a surface in the process of heterogeneous nucleation and contributes towards a rational definition of surfaces as effective nucleating agents.
Collapse
Affiliation(s)
- Veselina Marinova
- Department of Materials Science and Engineering, The University of Sheffield, Sheffield, S1 3JD, UK.
| | - Colin L Freeman
- Department of Materials Science and Engineering, The University of Sheffield, Sheffield, S1 3JD, UK.
| | - John H Harding
- Department of Materials Science and Engineering, The University of Sheffield, Sheffield, S1 3JD, UK.
| |
Collapse
|
9
|
Bal KM. Nucleation rates from small scale atomistic simulations and transition state theory. J Chem Phys 2021; 155:144111. [PMID: 34654300 DOI: 10.1063/5.0063398] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.
Collapse
Affiliation(s)
- Kristof M Bal
- Department of Chemistry and NANOlab Center of Excellence, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
| |
Collapse
|
10
|
Abstract
We describe a nonparametric approach for accurate determination of the slowest relaxation eigenvectors of molecular dynamics. The approach is blind as it uses no system specific information. In particular, it does not require a functional form with many parameters to closely approximate eigenvectors, e.g., linear combinations of molecular descriptors or a deep neural network, and thus no extensive expertise with the system. We suggest a rigorous and sensitive validation/optimality criterion for an eigenvector. The criterion uses only eigenvector time series and can be used to validate eigenvectors computed by other approaches. The power of the approach is illustrated on long atomistic protein folding trajectories. The determined eigenvectors pass the validation test at a time scale of 0.2 ns, much shorter than alternative approaches.
Collapse
Affiliation(s)
- Sergei V Krivov
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
11
|
Zhang C, Sibert EL, Gruebele M. A phase diagram for energy flow-limited reactivity. J Chem Phys 2021; 154:104301. [PMID: 33722023 DOI: 10.1063/5.0043665] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Intramolecular energy flow (also known as intramolecular vibrational redistribution or IVR) is often assumed in Rice-Ramsperger-Kassel-Marcus, transition state, collisional energy transfer, and other rate calculations not to be an impediment to reaction. In contrast, experimental spectroscopy, computational results, and models based on Anderson localization have shown that ergodicity is achieved rather slowly during molecular energy flow. The statistical assumption in rate theories might easily fail due to quantum localization. Here, we develop a simple model for the interplay of IVR and energy transfer and simulate the model with near-exact quantum dynamics for a 10-degree of freedom system composed of two five-mode molecular fragments. The calculations are facilitated by applying the van Vleck transformation to local random matrix models of the vibrational Hamiltonian. We find that there is a rather sharp "phase transition" as a function of molecular anharmonicity "a" between a region of facile energy transfer and a region limited by IVR and incomplete accessibility of the state space (classically, the phase space). The very narrow transition range of the order parameter a happens to lie right in the middle of the range expected for molecular torsion, bending, and stretching vibrations, thus demonstrating that reactive energy transfer dynamics several kBT above the thermal energy occurs not far from the localization boundary, with implications for controllability of reactions.
Collapse
Affiliation(s)
- Chenghao Zhang
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Edwin L Sibert
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Martin Gruebele
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
12
|
Khan SA, Dickson BM, Peters B. How fluxional reactants limit the accuracy/efficiency of infrequent metadynamics. J Chem Phys 2020; 153:054125. [DOI: 10.1063/5.0006980] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Salman A. Khan
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106-5080, USA
| | | | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Chemistry and Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
13
|
Numerical Analysis of Solute–Solvent Coupling Magnitude in the Thermally Backward Ring Closing Reaction of Spirooxazines. J SOLUTION CHEM 2020. [DOI: 10.1007/s10953-020-00986-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
14
|
Khan SA, Vandervelden CA, Scott SL, Peters B. Grafting metal complexes onto amorphous supports: from elementary steps to catalyst site populationsviakernel regression. REACT CHEM ENG 2020. [DOI: 10.1039/c9re00357f] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We use transition state theory, kernel regression, and population balance modeling techniques to model the grafting of metal complexes onto amorphous catalyst supports.
Collapse
Affiliation(s)
- Salman A. Khan
- Department of Chemical Engineering
- University of California
- Santa Barbara
- USA
| | | | - Susannah L. Scott
- Department of Chemical Engineering
- University of California
- Santa Barbara
- USA
- Department of Chemistry & Biochemistry
| | - Baron Peters
- Department of Chemical & Biomolecular Engineering
- University of Illinois at Urbana-Champaign
- Urbana
- USA
- Department of Chemistry
| |
Collapse
|
15
|
Yappert R, Kamat K, Peters B. The overdamped transmission coefficient: Recovering the true mean first passage time from an inaccurate reaction coordinate. J Chem Phys 2019; 151:184108. [DOI: 10.1063/1.5117237] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Ryan Yappert
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Kartik Kamat
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
16
|
Zaklika J, Komorowski L, Ordon P. Bond Fragility Spectra for the Double Proton-Transfer Reaction in the Formic Acid-Type Dimers. J Phys Chem A 2019; 123:4274-4283. [PMID: 31008601 DOI: 10.1021/acs.jpca.9b00595] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The newly developed method of fragility spectra for observation of bond breaking and formation upon a reaction has been applied to the canonical reaction series of the double proton transfer (DPT). Formic acid and its thio-analogues HCXYH (X, Y = O, S) have been chosen for the analysis. Very accurate linear correlations have been determined between the nondiagonal elements of the connectivity matrix, essential for the method, and the Wiberg bond orders for the corresponding bonds. Relation of the slope of this correlation to the global softness and to the atomic numbers of the bonded atoms has been proved, thus corroborating the c-DFT formula describing the fragility spectra. The electron density changes in bonds, as observed by the fragility spectra, are in harmony with the curvature diagrams reported by other authors.
Collapse
Affiliation(s)
- Jarosław Zaklika
- Department of Physical and Quantum Chemistry , Wrocław University of Science and Technology , Wyb. Wyspiańskiego 27 , 50-370 Wrocław , Poland
| | - Ludwik Komorowski
- Department of Physical and Quantum Chemistry , Wrocław University of Science and Technology , Wyb. Wyspiańskiego 27 , 50-370 Wrocław , Poland
| | - Piotr Ordon
- Department of Physics and Biophysics , Wrocław University of Environmental and Life Sciences , ul. Norwida 25 , 50-373 Wrocław , Poland
| |
Collapse
|
17
|
Craven GT, Chen R, Nitzan A. Upside/Downside statistical mechanics of nonequilibrium Brownian motion. II. Heat transfer and energy partitioning of a free particle. J Chem Phys 2018; 149:104103. [PMID: 30219017 DOI: 10.1063/1.5045361] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The energy partitioning during activation and relaxation events under steady-state conditions for a Brownian particle driven by multiple thermal reservoirs of different local temperatures is investigated. Specifically, we apply the formalism derived in Paper I [G. T. Craven and A. Nitzan, J. Chem. Phys. 148, 044101 (2018)] to examine the thermal transport properties of two sub-ensembles of Brownian processes, distinguished at any given time by the specification that all the trajectories in each group have, at that time, energy either above (upside) or below (downside) a preselected energy threshold. Dynamical properties describing energy accumulation and release during activation/relaxation events and relations for upside/downside energy partitioning between thermal reservoirs are derived. The implications for heat transport induced by upside and downside events are discussed.
Collapse
Affiliation(s)
- Galen T Craven
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Renai Chen
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
18
|
Zimmermann NER, Vorselaars B, Espinosa JR, Quigley D, Smith WR, Sanz E, Vega C, Peters B. NaCl nucleation from brine in seeded simulations: Sources of uncertainty in rate estimates. J Chem Phys 2018; 148:222838. [DOI: 10.1063/1.5024009] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Nils. E. R. Zimmermann
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Bart Vorselaars
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln LN6 7TS, United Kingdom
| | - Jorge R. Espinosa
- Departmento de Quimica-Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - David Quigley
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - William R. Smith
- Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario N1G2W1, Canada
- Department of Chemistry, University of Guelph, Guelph, Ontario N1G2W1, Canada
| | - Eduardo Sanz
- Departmento de Quimica-Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departmento de Quimica-Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Baron Peters
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
19
|
Craven GT, Nitzan A. Upside/Downside statistical mechanics of nonequilibrium Brownian motion. I. Distributions, moments, and correlation functions of a free particle. J Chem Phys 2018; 148:044101. [PMID: 29390855 DOI: 10.1063/1.5007854] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Statistical properties of Brownian motion that arise by analyzing, separately, trajectories over which the system energy increases (upside) or decreases (downside) with respect to a threshold energy level are derived. This selective analysis is applied to examine transport properties of a nonequilibrium Brownian process that is coupled to multiple thermal sources characterized by different temperatures. Distributions, moments, and correlation functions of a free particle that occur during upside and downside events are investigated for energy activation and energy relaxation processes and also for positive and negative energy fluctuations from the average energy. The presented results are sufficiently general and can be applied without modification to the standard Brownian motion. This article focuses on the mathematical basis of this selective analysis. In subsequent articles in this series, we apply this general formalism to processes in which heat transfer between thermal reservoirs is mediated by activated rate processes that take place in a system bridging them.
Collapse
Affiliation(s)
- Galen T Craven
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
20
|
Feng J, Shukla D. Characterizing Conformational Dynamics of Proteins Using Evolutionary Couplings. J Phys Chem B 2018; 122:1017-1025. [DOI: 10.1021/acs.jpcb.7b07529] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiangyan Feng
- Department
of Chemical and Biomolecular Engineering, ‡Center for Biophysics and Quantitative
Biology, §Department of Plant Biology, and ∥National Center for Supercomputing Applications, University of Illinois, Urbana, Illinois 61801, United States
| | - Diwakar Shukla
- Department
of Chemical and Biomolecular Engineering, ‡Center for Biophysics and Quantitative
Biology, §Department of Plant Biology, and ∥National Center for Supercomputing Applications, University of Illinois, Urbana, Illinois 61801, United States
| |
Collapse
|
21
|
Computational Analysis of Solute–Solvent Coupling Magnitude in the Z/E Isomerization Reaction of Nitroazobenzene and Benzylideneanilines. J SOLUTION CHEM 2018. [DOI: 10.1007/s10953-018-0711-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
Paesani F. Making Ice from Stacking-Disordered Crystallites. Chem 2017. [DOI: 10.1016/j.chempr.2017.12.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
23
|
Goldsmith BR, Peters B, Johnson JK, Gates BC, Scott SL. Beyond Ordered Materials: Understanding Catalytic Sites on Amorphous Solids. ACS Catal 2017. [DOI: 10.1021/acscatal.7b01767] [Citation(s) in RCA: 99] [Impact Index Per Article: 14.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bryan R. Goldsmith
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
- Department
of Chemical Engineering, University of Michigan, Ann Arbor, Michigan 48105, United States
| | - Baron Peters
- Department
of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - J. Karl Johnson
- Department
of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Bruce C. Gates
- Department
of Chemical Engineering, University of California, Davis, California 95616, United States
| | - Susannah L. Scott
- Department
of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
24
|
Chen R, Craven GT, Nitzan A. Electron-transfer-induced and phononic heat transport in molecular environments. J Chem Phys 2017; 147:124101. [DOI: 10.1063/1.4990410] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Renai Chen
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104,
USA
| | - Galen T. Craven
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104,
USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104,
USA
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
25
|
Zinovjev K, Tuñón I. Reaction coordinates and transition states in enzymatic catalysis. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1329] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Kirill Zinovjev
- Departament de Química FísicaUniversitat de València Valencia Spain
| | - Iñaki Tuñón
- Departament de Química FísicaUniversitat de València Valencia Spain
| |
Collapse
|
26
|
Craven GT, Nitzan A. Electrothermal Transistor Effect and Cyclic Electronic Currents in Multithermal Charge Transfer Networks. PHYSICAL REVIEW LETTERS 2017; 118:207201. [PMID: 28581803 DOI: 10.1103/physrevlett.118.207201] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2016] [Indexed: 06/07/2023]
Abstract
A theory is developed to describe the coupled transport of energy and charge in networks of electron donor-acceptor sites which are seated in a thermally heterogeneous environment, where the transfer kinetics are dominated by Marcus-type hopping rates. It is found that the coupling of heat and charge transfer in such systems gives rise to exotic transport phenomena which are absent in thermally homogeneous systems and cannot be described by standard thermoelectric relations. Specifically, the directionality and extent of thermal transistor amplification and cyclical electronic currents in a given network can be controlled by tuning the underlying temperature gradient in the system. The application of these findings toward the optimal control of multithermal currents is illustrated on a paradigmatic nanostructure.
Collapse
Affiliation(s)
- Galen T Craven
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
27
|
Craven GT, Nitzan A. Electron transfer at thermally heterogeneous molecule-metal interfaces. J Chem Phys 2017. [DOI: 10.1063/1.4971293] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Galen T. Craven
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
28
|
McGibbon RT, Husic BE, Pande VS. Identification of simple reaction coordinates from complex dynamics. J Chem Phys 2017; 146:044109. [PMID: 28147508 PMCID: PMC5272828 DOI: 10.1063/1.4974306] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2016] [Accepted: 01/05/2017] [Indexed: 11/14/2022] Open
Abstract
Reaction coordinates are widely used throughout chemical physics to model and understand complex chemical transformations. We introduce a definition of the natural reaction coordinate, suitable for condensed phase and biomolecular systems, as a maximally predictive one-dimensional projection. We then show that this criterion is uniquely satisfied by a dominant eigenfunction of an integral operator associated with the ensemble dynamics. We present a new sparse estimator for these eigenfunctions which can search through a large candidate pool of structural order parameters and build simple, interpretable approximations that employ only a small number of these order parameters. Example applications with a small molecule's rotational dynamics and simulations of protein conformational change and folding show that this approach can filter through statistical noise to identify simple reaction coordinates from complex dynamics.
Collapse
Affiliation(s)
- Robert T McGibbon
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Brooke E Husic
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
29
|
Lupi L, Peters B, Molinero V. Pre-ordering of interfacial water in the pathway of heterogeneous ice nucleation does not lead to a two-step crystallization mechanism. J Chem Phys 2016; 145:211910. [DOI: 10.1063/1.4961652] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Laura Lupi
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| | - Baron Peters
- Department of Chemical Engineering and Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
30
|
Feng X, Krylov AI. On couplings and excimers: lessons from studies of singlet fission in covalently linked tetracene dimers. Phys Chem Chem Phys 2016; 18:7751-61. [PMID: 26910414 DOI: 10.1039/c6cp00177g] [Citation(s) in RCA: 84] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Electronic factors controlling singlet fission (SF) rates are investigated in covalently linked dimers of tetracene. Using covalent linkers, relative orientation of the individual chromophores can be controlled, maximizing the rates of SF. Structures with coplanar and staggered arrangements of tetracene moieties are considered. The electronic structure calculations and three-state kinetic model for SF rates provide explanations for experimentally observed low SF yields in coplanar dimers and efficient SF in staggered dimers. The calculations illuminate the role of the excimer formation in SF process. The structural relaxation in the S1 state leads to the increased rate of the multi-exciton (ME) state formation, but impedes the second step, separation of the ME state into independent triplets. The slower second step reduces SF yield by allowing other processes, such as radiationless relaxation, to compete with triplet generation. The calculations of electronic couplings also suggest an increased rate of radiationless relaxation at the excimer geometries. Thus, the excimer serves as a trap of the ME state. The effect of covalent linkers on the electronic factors and SF rates is investigated. In all considered structures, the presence of the linker leads to larger couplings, however, the effect on the overall rate is less straightforward, since the linkers generally result in less favorable energetics. This complex behavior once again illustrates the importance of integrative approaches that evaluate the overall rate, rather than focusing on specific electronic factors such as energies or couplings.
Collapse
Affiliation(s)
- Xintian Feng
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA.
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA.
| |
Collapse
|
31
|
Abstract
Charge transfer is a fundamental process that underlies a multitude of phenomena in chemistry and biology. Recent advances in observing and manipulating charge and heat transport at the nanoscale, and recently developed techniques for monitoring temperature at high temporal and spatial resolution, imply the need for considering electron transfer across thermal gradients. Here, a theory is developed for the rate of electron transfer and the associated heat transport between donor-acceptor pairs located at sites of different temperatures. To this end, through application of a generalized multidimensional transition state theory, the traditional Arrhenius picture of activation energy as a single point on a free energy surface is replaced with a bithermal property that is derived from statistical weighting over all configurations where the reactant and product states are equienergetic. The flow of energy associated with the electron transfer process is also examined, leading to relations between the rate of heat exchange among the donor and acceptor sites as functions of the temperature difference and the electronic driving bias. In particular, we find that an open electron transfer channel contributes to enhanced heat transport between sites even when they are in electronic equilibrium. The presented results provide a unified theory for charge transport and the associated heat conduction between sites at different temperatures.
Collapse
Affiliation(s)
- Galen T Craven
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104; School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
32
|
Banushkina PV, Krivov SV. Optimal reaction coordinates. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1276] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Polina V. Banushkina
- Astbury Center for Structural Molecular Biology; Faculty of Biological Sciences, University of Leeds; Leeds UK
| | - Sergei V. Krivov
- Astbury Center for Structural Molecular Biology; Faculty of Biological Sciences, University of Leeds; Leeds UK
| |
Collapse
|
33
|
Affiliation(s)
- Baron Peters
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106;
| |
Collapse
|
34
|
Komorowski L, Ordon P, Jędrzejewski M. The reaction fragility spectrum. Phys Chem Chem Phys 2016; 18:32658-32663. [DOI: 10.1039/c6cp06519h] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report an original method that provides a new insight into the reaction mechanism by direct observation of bond breaking and formation. Example: HONS → ONSH.
Collapse
Affiliation(s)
- Ludwik Komorowski
- Department of Physical and Quantum Chemistry
- Wrocław University of Science and Technology
- 50-370 Wrocław
- Poland
| | - Piotr Ordon
- Department of Physics and Biophysics
- Wrocław University of Environmental and Life Sciences
- 50-375 Wrocław
- Poland
| | - Mateusz Jędrzejewski
- Department of Physical and Quantum Chemistry
- Wrocław University of Science and Technology
- 50-370 Wrocław
- Poland
| |
Collapse
|
35
|
Matyushov DV. Protein electron transfer: is biology (thermo)dynamic? JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:473001. [PMID: 26558324 DOI: 10.1088/0953-8984/27/47/473001] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Simple physical mechanisms are behind the flow of energy in all forms of life. Energy comes to living systems through electrons occupying high-energy states, either from food (respiratory chains) or from light (photosynthesis). This energy is transformed into the cross-membrane proton-motive force that eventually drives all biochemistry of the cell. Life's ability to transfer electrons over large distances with nearly zero loss of free energy is puzzling and has not been accomplished in synthetic systems. The focus of this review is on how this energetic efficiency is realized. General physical mechanisms and interactions that allow proteins to fold into compact water-soluble structures are also responsible for a rugged landscape of energy states and a broad distribution of relaxation times. Specific to a protein as a fluctuating thermal bath is the protein-water interface, which is heterogeneous both dynamically and structurally. The spectrum of interfacial fluctuations is a consequence of protein's elastic flexibility combined with a high density of surface charges polarizing water dipoles into surface nanodomains. Electrostatics is critical to the protein function and the relevant questions are: (i) What is the spectrum of interfacial electrostatic fluctuations? (ii) Does the interfacial biological water produce electrostatic signatures specific to proteins? (iii) How is protein-mediated chemistry affected by electrostatics? These questions connect the fluctuation spectrum to the dynamical control of chemical reactivity, i.e. the dependence of the activation free energy of the reaction on the dynamics of the bath. Ergodicity is often broken in protein-driven reactions and thermodynamic free energies become irrelevant. Continuous ergodicity breaking in a dense spectrum of relaxation times requires using dynamically restricted ensembles to calculate statistical averages. When applied to the calculation of the rates, this formalism leads to the nonergodic activated kinetics, which extends the transition-state theory to dynamically dispersive media. Releasing the grip of thermodynamics in kinetic calculations through nonergodicity provides the mechanism for an efficient optimization between reaction rates and the spectrum of relaxation times of the protein-water thermal bath. Bath dynamics, it appears, play as important role as the free energy in optimizing biology's performance.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, USA
| |
Collapse
|
36
|
Banushkina PV, Krivov SV. Nonparametric variational optimization of reaction coordinates. J Chem Phys 2015; 143:184108. [DOI: 10.1063/1.4935180] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Polina V. Banushkina
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Sergei V. Krivov
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
37
|
Zinovjev K, Tuñón I. Transition state ensemble optimization for reactions of arbitrary complexity. J Chem Phys 2015; 143:134111. [DOI: 10.1063/1.4931596] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|