Researcher profile

Todd J. Martínez

Todd J. Martínez contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.chem-ph physics.comp-ph quant-ph

Trust snapshot

Quick read

Trust 19 - Baseline

5works

0followers

3topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2018arXiv

Reaction dynamics of cyanohydrins with hydrosulfide in water

We studied the reaction dynamics of a proposed prebiotic reaction theoretically. The chemical process involves acetone cyanohydrin or formalcyanohydrin reacting with hydrosulfide in an aqueous environment. Rate constants and populations of reactant and product bimolecular geometric orientations for the reactions were obtained by using density functional theory for the energies, transition state theory for the rates and matrix exponentiation as well as the hybrid tau-leaping algorithm for the population dynamics. The role of including the solvent explicitly versus implicitly was also investigated. We found that adding water or hydrogen sulfide molecules explicitly lowers the activation energy barrier and leads to a more efficient reaction pathway. In particular, hydrogen sulfide was a better proton donor. Further, when adding the solvent explicitly a synchronous reaction mechanism was observed while a concerted reaction mechanism was seen when using an implicit solvent model. Finally, we studied the role of including more than one reactant and product bimolecular orientation geometry in the dynamics. Including all pathways, reactant to reactant, product to product, reactant to product and product to reactant reduced the reaction half-time by two orders of magnitude and was therefore determined to be fundamental for an accurate description of the dynamics. Overall, we found that most reactions which involve formalcyanohydrin occur more rapidly or at the same speed as reactions which involve acetone cyanohydrin at room temperature.

preprint2015arXiv

Quantum Chemistry for Solvated Molecules on Graphical Processing Units (GPUs)using Polarizable Continuum Models

The conductor-like polarization model (C-PCM) with switching/Gaussian smooth discretization is a widely used implicit solvation model in chemical simulations. However, its application in quantum mechanical calculations of large-scale biomolecular systems can be limited by computational expense of both the gas phase electronic structure and the solvation interaction. We have previously used graphical processing units (GPUs) to accelerate the first of these steps. Here, we extend the use of GPUs to accelerate electronic structure calculations including C-PCM solvation. Implementation on the GPU leads to significant acceleration of the generation of the required integrals for C-PCM. We further propose two strategies to improve the solution of the required linear equations: a dynamic convergence threshold and a randomized block-Jacobi preconditioner. These strategies are not specific to GPUs and are expected to be beneficial for both CPU and GPU implementations. We benchmark the performance of the new implementation using over 20 small proteins in solvent environment. Using a single GPU, our method evaluates the C-PCM related integrals and their derivatives more than 10X faster than a conventional CPU based implementation. Our improvements to the linear solver provide a further 3X acceleration. The overall calculations including C-PCM solvation require typically 20-40% more effort than their gas phase counterparts for moderate basis set and molecule surface discretization level. The relative cost of the C-PCM solvation correction decreases as the basis sets and/or cavity radii increase. Therefore description of solvation with this model should be routine. We also discuss applications to the study of the conformational landscape of an amyloid fibril.

preprint2013arXiv

Relation of exact Gaussian basis methods to the dephasing representation: Theory and application to time-resolved electronic spectra

We recently showed that the Dephasing Representation (DR) provides an efficient tool for computing ultrafast electronic spectra and that further acceleration is possible with cellularization [M. Šulc and J. Vaníček, Mol. Phys. 110, 945 (2012)]. Here we focus on increasing the accuracy of this approximation by first implementing an exact Gaussian basis method, which benefits from the accuracy of quantum dynamics and efficiency of classical dynamics. Starting from this exact method, the DR is derived together with ten other methods for computing time-resolved spectra with intermediate accuracy and efficiency. These methods include the Gaussian DR, an exact generalization of the DR, in which trajectories are replaced by communicating frozen Gaussian basis functions evolving classically with an average Hamiltonian. The newly obtained methods are tested numerically on time correlation functions and time-resolved stimulated emission spectra in the harmonic potential, pyrazine S0/S1 model, and quartic oscillator. Numerical results confirm that both the Gaussian basis method and the Gaussian DR increase the accuracy of the DR. Surprisingly, in chaotic systems the Gaussian DR can outperform the presumably more accurate Gaussian basis method, in which the two bases are evolved separately.

preprint2009arXiv

A note on the role of charge conservation in electronegativity equalization and its implications for the translational symmetries of electrostatic properties in fluctuating-charge models

An analytical solution of fluctuating-charge models using Gaussian elimination allows us to isolate the contribution of charge conservation effects in determining the charge distribution. We use this analytical solution to calculate dipole moments and polarizabilities and show that charge conservation plays a critical role in maintaining the correct translational properties of the electrostatic properties predicted by these models.

preprint2008arXiv

A unified theoretical framework for fluctuating-charge models in atom-space and in bond-space

Our previously introduced QTPIE (charge transfer with polarization current equilibration) model (J. Chen and T. J. Martinez, Chem. Phys. Lett. 438, 315 (2007)) is a fluctuating-charge model with correct asymptotic behavior. Unlike most other fluctuating-charge models, QTPIE is formulated in terms of charge-transfer variables and pairwise electronegativities, not atomic charge variables and electronegativities. The pairwise character of the electronegativities in QTPIE avoids spurious charge transfer when bonds are broken. However, the increased number of variables leads to considerable computational expense and a rank-deficient set of working equations, which is numerically inconvenient. Here, we show that QTPIE can be exactly reformulated in terms of atomic charge variables, leading to a considerable reduction in computational complexity. The transformation between atomic and bond variables is generally applicable to arbitrary fluctuating charge models, and uncovers an underlying topological framework that can be used to understand the relation between fluctuating-charge models and the classical theory of electrical circuits.