Showing 1–30 of 30 results for author: Martínez, T J

Attosecond Probing of Coherent Vibrational Dynamics in CBr$_4$

Authors: Jen-Hao Ou, Diptarka Hait, Patrick Rupprecht, John E. Beetar, Todd J. Martínez, Stephen R. Leone

Abstract: A coherent vibrational wavepacket is launched and manipulated in the symmetric stretch (a$_1$) mode of CBr$_4$, by impulsive stimulated Raman scattering from non-resonant 400 nm laser pump pulses with various peak intensities on the order of tens of 10$^{12}$ W/cm$^2$. Extreme ultraviolet (XUV) attosecond transient absorption spectroscopy (ATAS) records the wavepacket dynamics as temporal oscillat… ▽ More A coherent vibrational wavepacket is launched and manipulated in the symmetric stretch (a$_1$) mode of CBr$_4$, by impulsive stimulated Raman scattering from non-resonant 400 nm laser pump pulses with various peak intensities on the order of tens of 10$^{12}$ W/cm$^2$. Extreme ultraviolet (XUV) attosecond transient absorption spectroscopy (ATAS) records the wavepacket dynamics as temporal oscillations in XUV absorption energy at the bromine M$_{4,5}$ 3d$_{3/2,5/2}$ edges around 70 eV. The results are augmented by nuclear time-dependent Schrödinger equation simulations. Slopes of the (Br-3d$_{3/2,5/2}$)$^{-1}$10a$_1^*$ core-excited state potential energy surface (PES) along the a$_1$ mode are calculated to be -9.4 eV/Å from restricted open-shell Kohn-Sham calculations. Using analytical relations derived for the small-displacement limit with the calculated slopes of the core-excited state PES, a deeper insight into the vibrational dynamics is obtained by retrieving the experimental excursion amplitude of the vibrational wavepacket and the amount of population transferred to the vibrational first-excited state, as a function of pump-pulse peak intensity. Experimentally, the results show that XUV ATAS is capable of easily resolving oscillations in the XUV absorption energy on the order of few to tens of meV and tens of femtosecond time precision, limited only by the averaging times in the experimental scans. This corresponds to oscillations of C-Br bond length on the order of 10$^{-4}$ to 10$^{-3}$ Å. The results and the analytic relationships offer a clear physical picture, on multiple levels of understanding, for how the pump-pulse intensity controls the vibrational dynamics launched by non-resonant ISRS in the small-displacement limit. △ Less

Submitted 1 August, 2024; originally announced August 2024.

Extending GPU-Accelerated Gaussian Integrals in the TeraChem Software Package to f Type Orbitals: Implementation and Applications

Authors: Yuanheng Wang, Diptarka Hait, K. Grace Johnson, O. Jonathan Fajen, Rubén D. Guerrero, Todd J. Martínez

Abstract: The increasing availability of GPUs for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. The complexity of integral kernels for high angular momentum basis functions however often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report implementation of $f$ function support i… ▽ More The increasing availability of GPUs for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. The complexity of integral kernels for high angular momentum basis functions however often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report implementation of $f$ function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles (CCSD) calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH$_3$) with complete active space self-consistent field (CASSCF) are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules. △ Less

Submitted 15 August, 2024; v1 submitted 21 June, 2024; originally announced June 2024.

Comments: Updated benchmarks from previous version

Unexpected hydrogen dissociation in thymine: predictions from a novel coupled cluster theory

Authors: Eirik F. Kjønstad, O. Jonathan Fajen, Alexander C. Paul, Sara Angelico, Dennis Mayer, Markus Gühr, Thomas J. A. Wolf, Todd J. Martínez, Henrik Koch

Abstract: The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright $ππ^*$ state to a dark $nπ^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we sim… ▽ More The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright $ππ^*$ state to a dark $nπ^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the ultrafast dynamics in thymine at the highest level of theory to date, performing wavepacket dynamics with a new coupled cluster method. Our simulation confirms an ultrafast $ππ^*$ to $nπ^*$ transition ($τ = 41 \pm 14$ fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with the experimental results. Our simulation also predicts an as-yet uncharacterized photochemical pathway: a $πσ^*$ channel that leads to hydrogen dissociation at one of the two N-H bonds in thymine. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel. △ Less

Submitted 7 March, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 42 pages, 23 figures

Prediction of Photodynamics of 200 nm Excited Cyclobutanone with Linear Response Electronic Structure and Ab Initio Multiple Spawning

Authors: Diptarka Hait, Dean Lahana, O. Jonathan Fajen, Amiel S. P. Paz, Pablo A. Unzueta, Bhaskar Rana, Lixin Lu, Yuanheng Wang, Todd J. Martinez

Abstract: Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S$_2$). A picosecond of nonadiabatic dynam… ▽ More Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S$_2$). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. We use both time dependent density functional theory and equation-of-motion coupled cluster for the underlying electronic structure theory. We find that the lifetime of the S$_2$ state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted UED spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Journal ref: J. Chem. Phys. 160, 244101 (2024)

Femtosecond electronic and hydrogen structural dynamics in ammonia imaged with ultrafast electron diffraction

Authors: Elio G. Champenois, Nanna H. List, Matthew Ware, Mathew Britton, Philip H. Bucksbaum, Xinxin Cheng, Martin Centurion, James P. Cryan, Ruaridh Forbes, Ian Gabalski, Kareem Hegazy, Matthias C. Hoffmann, Andrew J. Howard, Fuhao Ji, Ming-Fu Lin, J. Pedro Nunes, Xiaozhe Shen, Jie Yang, Xijie Wang, Todd J. Martinez, Thomas J. A. Wolf

Abstract: Directly imaging structural dynamics involving hydrogen atoms by ultrafast diffraction methods is complicated by their low scattering cross-sections. Here we demonstrate that megaelectronvolt ultrafast electron diffraction is sufficiently sensitive to follow hydrogen dynamics in isolated molecules. In a study of the photodissociation of gas phase ammonia, we simultaneously observe signatures of th… ▽ More Directly imaging structural dynamics involving hydrogen atoms by ultrafast diffraction methods is complicated by their low scattering cross-sections. Here we demonstrate that megaelectronvolt ultrafast electron diffraction is sufficiently sensitive to follow hydrogen dynamics in isolated molecules. In a study of the photodissociation of gas phase ammonia, we simultaneously observe signatures of the nuclear and corresponding electronic structure changes resulting from the dissociation dynamics in the time-dependent diffraction. Both assignments are confirmed by ab initio simulations of the photochemical dynamics and the resulting diffraction observable. While the temporal resolution of the experiment is insufficient to resolve the dissociation in time, our results represent an important step towards the observation of proton dynamics in real space and time. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Rehybridization dynamics into the pericyclic minimum of an electrcyclic reaction imaged in real-time

Authors: Yusong Liu, David M. Sanchez, Matthew R. Ware, Elio G. Champenois, Jie Yang, J. Pedro F. Nunes, Andrew Attar, Martin Centurion, James P. Cryan, Ruaridh G. Forbes, Kareem Hegazy, Matthias C. Hoffmann, Fuhao Ji, Ming-Fu Lin, Duan Luo, Sajib K. Saha, Xiaozhe Shen, Xijie Wang, Todd J. Martínez, Thomas J. A. Wolf

Abstract: Electrocyclic reactions are characterized by the concerted formation and cleavage of both σ and π bonds through a cyclic structure. This structure is known as a pericyclic transition state for thermal reactions and a pericyclic minimum in the excited state for photochemical reactions. However, the structure of the pericyclic geometry has yet to be observed experimentally. We use a combination of u… ▽ More Electrocyclic reactions are characterized by the concerted formation and cleavage of both σ and π bonds through a cyclic structure. This structure is known as a pericyclic transition state for thermal reactions and a pericyclic minimum in the excited state for photochemical reactions. However, the structure of the pericyclic geometry has yet to be observed experimentally. We use a combination of ultrafast electron diffraction and excited state wavepacket simulations to image structural dynamics through the pericyclic minimum of a photochemical electrocyclic ring-opening reaction in the molecule α-terpinene. The structural motion into the pericyclic minimum is dominated by rehybridization of two carbon atoms, which is required for the transformation from two to three conjugated π bonds. The σ bond dissociation largely happens after internal conversion from the pericyclic minimum to the electronic ground state. These findings may be transferrable to electrocyclic reactions in general. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Comments: Combined manuscript and supplementary information

Rank-reduced coupled-cluster III. Tensor hypercontraction of the doubles amplitudes

Authors: Edward G. Hohenstein, B. Scott Fales, Robert M. Parrish, Todd J. Martinez

Abstract: We develop a quartic-scaling implementation of coupled-cluster singles and doubles based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the g… ▽ More We develop a quartic-scaling implementation of coupled-cluster singles and doubles based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the gain in computational efficiency as it is sufficient, in conjunction with a Cholesky decomposition of the ERIs, to reduce the computational complexity of most contributions to the CCSD amplitude equations. Further THC factorization of the ERIs reduces the complexity of certain terms arising from nested commutators between the doubles excitation operator and the two-electron operator. We implement this new algorithm using graphical processing units (GPUs) and demonstrate that it enables CCSD calculations for molecules with 250 atoms and 2500 basis functions using a single computer node. Further, we show that the new method computes correlation energies with comparable accuracy to the underlying RR-CCSD method. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Conformer-specific Chemistry Imaged in Real Space and Time

Authors: E. G. Champenois, D. M. Sanchez, J. Yang, J. P. F. Nunes, A. Attar, M. Centurion, R. Forbes, M. Gühr, K. Hegazy, F. Ji, S. K. Saha, Y. Liu, M. -F. Lin, D. Luo, B. Moore, X. Shen, M. R. Ware, X. J. Wang, T. J. Martínez, T. J. A. Wolf

Abstract: Conformational isomers or conformers of molecules play a decisive role in chemistry and biology. However, experimental methods to investigate chemical reaction dynamics are typically not conformer-sensitive. Here, we report on a gas-phase megaelectronvolt ultrafast electron diffraction investigation of α-phellandrene undergoing an electrocyclic ring-opening reaction. We directly image the evolutio… ▽ More Conformational isomers or conformers of molecules play a decisive role in chemistry and biology. However, experimental methods to investigate chemical reaction dynamics are typically not conformer-sensitive. Here, we report on a gas-phase megaelectronvolt ultrafast electron diffraction investigation of α-phellandrene undergoing an electrocyclic ring-opening reaction. We directly image the evolution of a specific set of α-phellandrene conformers into the product isomer predicted by the Woodward-Hoffmann rules in real space and time. Our experimental results are in quantitative agreement with nonadiabatic quantum molecular dynamics simulations, which provide unprecedented detail of how conformation influences time scale and quantum efficiency of photoinduced ring-opening reactions. Due to the prevalence of large numbers of conformers in organic chemistry, our findings impact our general understanding of reaction dynamics in chemistry and biology. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Transient Resonant Auger-Meitner Spectra of Photoexcited Thymine

Authors: Thomas J. A. Wolf, Alexander C. Paul, Sarai D. Folkestad, Rolf H. Myhre, James P. Cryan, Nora Berrah, Phil H. Bucksbaum, Sonia Coriani, Giacomo Coslovich, Raimund Feifel, Todd J. Martinez, Stefan P. Moeller, Melanie Mucke, Razib Obaid, Oksana Plekan, Richard J. Squibb, Henrik Koch, Markus Gühr

Abstract: We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a ππ* excited state, thymine is known to undergo internal conversion to an… ▽ More We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a ππ* excited state, thymine is known to undergo internal conversion to an nπ* excited state with a strong resonance at the oxygen K-edge, red-shifted from the ground state π* resonances of thymine (see our previous study Wolf et al., Nat. Commun., 2017, 8, 29). We resolve and compare the Auger-Meitner electron spectra associated both with the excited state and ground state resonances, and distinguish participator and spectator decay contributions. Furthermore, we observe simultaneously with the decay of the nπ* state signatures the appearance of additional resonant Auger-Meitner contributions at photon energies between the nπ* state and the ground state resonances. We assign these contributions to population transfer from the nπ* state to a ππ* triplet state via intersystem crossing on the picosecond timescale based on simulations of the X-ray absorption spectra in the vibrationally hot triplet state. Moreover, we identify signatures from the initially excited ππ* singlet state which we have not observed in our previous study. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Matter-wave interference of a native polypeptide

Authors: Armin Shayeghi, Philipp Rieser, Georg Richter, Ugur Sezer, Jonas H. Rodewald, Philipp Geyer, Todd. J. Martinez, Markus Arndt

Abstract: The de Broglie wave nature of matter is a paradigmatic example of fundamental quantum physics and enables precise measurements of forces, fundamental constants and even material properties. However, even though matter-wave interferometry is nowadays routinely realized in many laboratories, this feat has remained an outstanding challenge for the vast class of native polypeptides, the building block… ▽ More The de Broglie wave nature of matter is a paradigmatic example of fundamental quantum physics and enables precise measurements of forces, fundamental constants and even material properties. However, even though matter-wave interferometry is nowadays routinely realized in many laboratories, this feat has remained an outstanding challenge for the vast class of native polypeptides, the building blocks of life, which are ubiquitous in biology but fragile and difficult to handle. Here, we demonstrate the quantum wave nature of gramicidin, a natural antibiotic composed of 15 amino acids. Femtosecond laser desorption of a thin biomolecular film with intensities up to 1~TW/cm$^2$ transfers these molecules into a cold noble gas jet. Even though the peptide's de Broglie wavelength is as tiny as 350~fm, the molecular coherence is delocalized over more than 20 times the molecular size in our all-optical time-domain Talbot-Lau interferometer. We compare the observed interference fringes for two different interference orders with a model that includes both a rigorous treatment of the peptide's quantum wave nature as well as a quantum chemical assessment of its optical properties to distinguish our result from classical predictions. The successful realization of quantum optics with this polypeptide as a prototypical biomolecule paves the way for quantum-assisted molecule metrology and in particular the optical spectroscopy of a large class of biologically relevant molecules. △ Less

Submitted 31 October, 2019; originally announced October 2019.

Intermolecular Vibrations Drive Ultrafast Singlet Fission

Authors: Hong-Guang Duan, Ajay Jha 1, Xin Li, Vandana Tiwari, Hanyang Ye, Pabitra K. Nayak, Xiao-Lei Zhu, Zheng Li, Todd J. Martinez, Michael Thorwart, R. J. Dwayne Miller

Abstract: Singlet fission is a spin-allowed exciton multiplication process in organic semiconductors that converts one spin-singlet exciton to two triplet excitons. It offers the potential to enhance solar energy conversion by circumventing the Shockley-Queisser limit on efficiency. Recently, the mechanism of the primary singlet fission process in pentacene and its derivatives have been extensively investig… ▽ More Singlet fission is a spin-allowed exciton multiplication process in organic semiconductors that converts one spin-singlet exciton to two triplet excitons. It offers the potential to enhance solar energy conversion by circumventing the Shockley-Queisser limit on efficiency. Recently, the mechanism of the primary singlet fission process in pentacene and its derivatives have been extensively investigated, however, the nature of the primary ultrafast process in singlet fission is still a matter of debate. Here, we study the singlet fission process in a pentacene film by employing a combination of transient-grating (TG) and two-dimensional (2D) electronic spectroscopy complemented by quantum chemical and nonadiabatic dynamics calculations. The high sensitivity of heterodyne detected TG spectroscopy enabled us to capture the vibrational coherence and to show that it mediates the transition from the singlet excited electronic state to the triplet-pair state. This coherent process is further examined by 2D electronic spectroscopy. Detailed analysis of the experimental data reveals that significant vibronic couplings of a few key modes in the low- and high-frequency region connect the excited singlet and triplet-pair states. Based on quantum chemical calculations, we identify these key intermolecular rocking modes along the longitudinal molecular axis between the pentacene molecules. They play the essential role of an electronic bridge between the singlet and triplet-pair states. Along with high-frequency local vibrations acting as tuning modes, these rocking motions drive the ultrafast dynamics at the multidimensional conical intersection in the singlet fission process. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Diffractive imaging of dissociation and ground state dynamics in a complex molecule

Authors: Kyle Wilkin, Robert Parrish, Jie Yang, Thomas J. A. Wolf, Pedro Nunes, Markus Guehr, Renkai Li, Xiaozhe Shen, Qiang Zheng, Xijie Wang, Todd J. Martinez, Martin Centurion

Abstract: We have investigated the structural dynamics in photoexcited 1,2-diiodotetrafluoroethane molecules (C2F4I2) in the gas phase experimentally using ultrafast electron diffraction and theoretically using FOMO-CASCI excited state dynamics simulations. The molecules are excited by an ultra-violet femtosecond laser pulse to a state characterized by a transition from the iodine 5p orbital to a mixed 5p||… ▽ More We have investigated the structural dynamics in photoexcited 1,2-diiodotetrafluoroethane molecules (C2F4I2) in the gas phase experimentally using ultrafast electron diffraction and theoretically using FOMO-CASCI excited state dynamics simulations. The molecules are excited by an ultra-violet femtosecond laser pulse to a state characterized by a transition from the iodine 5p orbital to a mixed 5p|| hole and CF2 antibonding orbital, which results in the cleavage of one of the carbon-iodine bonds. We have observed, with sub-Angstrom resolution, the motion of the nuclear wavepacket of the dissociating iodine atom followed by coherent vibrations in the electronic ground state of the C2F4I radical. The radical reaches a stable classical (non-bridged) structure in less than 200 fs. △ Less

Submitted 26 July, 2019; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: 13 pages, 11 figures

Journal ref: Phys. Rev. A 100, 023402 (2019)

Imaging the Photochemical Ring-Opening of 1,3-Cyclohexadiene by Ultrafast Electron Diffraction

Authors: T. J. A. Wolf, D. M. Sanchez, J. Yang, R. M. Parrish, J. P. F. Nunes, M. Centurion, R. Coffee, J. P. Cryan, M. Gühr, K. Hegazy, A. Kirrander, R. K. Li, J. Ruddock, X. Shen, T. Veccione, S. P. Weathersby, P. M. Weber, K. Wilkin, H. Yong, Q. Zheng, X. J. Wang, M. P. Minitti, T. J. Martínez

Abstract: The ultrafast photoinduced ring-opening of 1,3-cyclohexadiene constitutes a textbook example of electrocyclic reactions in organic chemistry and a model for photobiological reactions in vitamin D synthesis. Here, we present direct and unambiguous observation of the ring-opening reaction path on the femtosecond timescale and sub-Ångström length scale by megaelectronvolt ultrafast electron diffracti… ▽ More The ultrafast photoinduced ring-opening of 1,3-cyclohexadiene constitutes a textbook example of electrocyclic reactions in organic chemistry and a model for photobiological reactions in vitamin D synthesis. Here, we present direct and unambiguous observation of the ring-opening reaction path on the femtosecond timescale and sub-Ångström length scale by megaelectronvolt ultrafast electron diffraction. We follow the carbon-carbon bond dissociation and the structural opening of the 1,3-cyclohexadiene ring by direct measurement of time-dependent changes in the distribution of interatomic distances. We observe a substantial acceleration of the ring-opening motion after internal conversion to the ground state due to steepening of the electronic potential gradient towards the product minima. The ring-opening motion transforms into rotation of the terminal ethylene groups in the photoproduct 1,3,5-hexatriene on the sub-picosecond timescale. Our work demonstrates the potential of megaelectronvolt ultrafast electron diffraction to elucidate photochemical reaction paths in organic chemistry. △ Less

Submitted 5 October, 2018; originally announced October 2018.

Reaction dynamics of cyanohydrins with hydrosulfide in water

Authors: Stéphanie Valleau, Todd J. Martínez

Abstract: 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 th… ▽ More 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. △ Less

Submitted 22 June, 2018; originally announced June 2018.

Crossing conditions in coupled cluster theory

Authors: Eirik F. Kjønstad, Rolf H. Myhre, Todd J. Martinez, Henrik Koch

Abstract: We derive the crossing conditions at conical intersections between electronic states in coupled cluster theory, and show that if the coupled cluster Jacobian matrix is nondefective, two (three) independent conditions are correctly placed on the nuclear degrees of freedom for an inherently real (complex) Hamiltonian. Calculations using coupled cluster theory on an $2 {^{1}}A' / 3 {^{1}}A'$ conical… ▽ More We derive the crossing conditions at conical intersections between electronic states in coupled cluster theory, and show that if the coupled cluster Jacobian matrix is nondefective, two (three) independent conditions are correctly placed on the nuclear degrees of freedom for an inherently real (complex) Hamiltonian. Calculations using coupled cluster theory on an $2 {^{1}}A' / 3 {^{1}}A'$ conical intersection in hypofluorous acid illustrate the nonphysical artifacts associated with defects at accidental same-symmetry intersections. In particular, the observed intersection seam is folded about a space of the correct dimensionality, indicating that minor modifications to the theory are required for it to provide a correct description of conical intersections in general. We find that an accidental symmetry allowed $1 {^{1}}A" / 2 {^{1}}A"$ intersection in hydrogen sulfide is properly described, showing no artifacts as well as linearity of the energy gap to first order in the branching plane. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Comments: 9 pages and 4 figures

Probing ultrafast ππ*/nπ* internal conversion in organic chromophores via K-edge resonant absorption

Authors: T. J. A. Wolf, R. H. Myhre, J. P. Cryan, S. Coriani, R. J. Squibb, A. Battistoni, N. Berrah, C. Bostedt, P. Bucksbaum, G. Coslovich, R. Feifel, K. J. Gaffney, J. Grilj, T. J. Martinez, S. Miyabe, S. P. Moeller, M. Mucke, A. Natan, R. Obaid, T. Osipov, O. Plekan, S. Wang, H. Koch, M. Gühr

Abstract: Organic chromophores with heteroatoms possess an important excited state relaxation channel from an optically allowed ππ* to a dark nπ*state. We exploit the element and site specificity of soft x-ray absorption spectroscopy to selectively follow the electronic change during the ππ*/nπ* internal conversion. As a hole forms in the n orbital during ππ*/nπ* internal conversion, the near edge x-ray abs… ▽ More Organic chromophores with heteroatoms possess an important excited state relaxation channel from an optically allowed ππ* to a dark nπ*state. We exploit the element and site specificity of soft x-ray absorption spectroscopy to selectively follow the electronic change during the ππ*/nπ* internal conversion. As a hole forms in the n orbital during ππ*/nπ* internal conversion, the near edge x-ray absorption fine structure (NEXAFS) spectrum at the heteroatom K-edge exhibits an additional resonance. We demonstrate the concept with the nucleobase thymine, a prototypical heteroatomic chromophore. With the help of time resolved NEXAFS spectroscopy at the oxygen K-edge, we unambiguously show that ππ*/nπ* internal conversion takes place within (60 \pm 30) fs. High-level coupled cluster calculations on the isolated molecules used in the experiment confirm the superb electronic structure sensitivity of this new method for excited state investigations. △ Less

Submitted 25 October, 2016; originally announced October 2016.

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

Authors: Fang Liu, Nathan Luehr, Heather J. Kulik, Todd J. Martínez

Abstract: 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 graphic… ▽ More 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. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: 36 pages, 12 figures

How large should the QM region be in QM/MM calculations? The case of catechol O-methyltransferase

Authors: Heather J. Kulik, Jianyu Zhang, Judith P. Klinman, Todd J. Martinez

Abstract: Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties i… ▽ More Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a representative example of a methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500-600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (ca. 11-16 residues or 200-300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes. △ Less

Submitted 3 August, 2016; v1 submitted 21 May, 2015; originally announced May 2015.

Comments: 27 pages, 10 figures

Multiple Time Step Integrators in Ab Initio Molecular Dynamics

Authors: Nathan Luehr, Thomas E. Markland, Todd J. Martinez

Abstract: Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties… ▽ More Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we show that such a timescale separation is possible using two different schemes: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy. △ Less

Submitted 10 November, 2013; originally announced December 2013.

Comments: 28 pages, 7 Figures

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

Authors: Miroslav Šulc, Henar Hernández, Todd J. Martínez, Jiří Vaníček

Abstract: 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 qu… ▽ More 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. △ Less

Submitted 3 June, 2013; originally announced June 2013.

Comments: 15 pages, 7 figures

Journal ref: J. Chem. Phys.139, 034112 (2013)

Tensor hypercontraction: A universal technique for the resolution of matrix elements of local, finite-range $N$-body potentials in many-body quantum problems

Authors: Robert M. Parrish, Edward G. Hohenstein, Nicolas F. Schunck, C. David Sherrill, Todd J. Martinez

Abstract: Configuration-space matrix elements of N-body potentials arise naturally and ubiquitously in the Ritz-Galerkin solution of many-body quantum problems. For the common specialization of local, finite-range potentials, we develop the eXact Tensor HyperContraction (X-THC) method, which provides a quantized renormalization of the coordinate-space form of the N-body potential, allowing for a highly sepa… ▽ More Configuration-space matrix elements of N-body potentials arise naturally and ubiquitously in the Ritz-Galerkin solution of many-body quantum problems. For the common specialization of local, finite-range potentials, we develop the eXact Tensor HyperContraction (X-THC) method, which provides a quantized renormalization of the coordinate-space form of the N-body potential, allowing for a highly separable tensor factorization of the configuration-space matrix elements. This representation allows for substantial computational savings in chemical, atomic, and nuclear physics simulations, particularly with respect to difficult "exchange-like" contractions. △ Less

Submitted 4 September, 2013; v1 submitted 21 January, 2013; originally announced January 2013.

Comments: Third version of the manuscript after referee's comments. In press in PRL. Main text: 4 pages, 2 figures, 1 table; Supplemental material (also included): 14 pages, 2 figures, 2 tables

Journal ref: Phys. Rev. Lett. 111, 132505 (2013)

Ultrafast Internal Conversion in Ethylene. II. Mechanisms and Pathways for Quenching and Hydrogen Elimination

Authors: T. K. Allison, H. Tao, W. J. Glover, T. W. Wright, A. M. Stooke, C. Khurmi, J. van Tilborg, Y. Liu, R. W. Falcone, T. J. Martinez, A. Belkacem

Abstract: Through a combined experimental and theoretical approach, we study the nonadiabatic dynamics of the prototypical ethylene (C$_2$H$_4$) molecule upon $π\rightarrow π^*$ excitation with 161 nm light. Using a novel experimental apparatus, we combine femtosecond pulses of vacuum ultraviolet (VUV) and extreme ultraviolet (XUV) radiation with variable delay to perform time resolved photo-ion fragment sp… ▽ More Through a combined experimental and theoretical approach, we study the nonadiabatic dynamics of the prototypical ethylene (C$_2$H$_4$) molecule upon $π\rightarrow π^*$ excitation with 161 nm light. Using a novel experimental apparatus, we combine femtosecond pulses of vacuum ultraviolet (VUV) and extreme ultraviolet (XUV) radiation with variable delay to perform time resolved photo-ion fragment spectroscopy. In this second part of a two part series, the extreme ultraviolet (17 eV$ < h ν< 23$ eV) probe pulses are sufficiently energetic to break the C-C bond in photoionization, or photoionize the dissociation products of the vibrationally hot ground state. The experimental data is directly compared to ab initio molecular dynamics simulations accounting for both the pump and probe steps. Enhancements of the CH$_2^+$ and CH$_3^+$ photoion fragment yields, corresponding to molecules photoionized in ethylene (CH$_2$CH$_2$) and ethylidene (CH$_3$CH) like geometries are observed within 100 fs after $π\rightarrow π^*$ excitation. Quantitative agreement between theory and experiment on the relative CH$_2^+$ and CH$_3^+$ yields provides experimental confirmation of the theoretical prediction of two distinct transition states and their branching ratio (Tao, et al. J. Phys. Chem. A. 113, 13656 (2009)). Fast, non-statistical, elimination of H$_2$ molecules and H atoms is observed in the time resolved H$_2^+$ and H$^+$ signals. △ Less

Submitted 1 October, 2011; originally announced October 2011.

Control of 1,3-Cyclohexadiene Photoisomerization Using Light-Induced Conical Intersections

Authors: Jaehee Kim, Hongli Tao, James L. White, Vladimir S. Petrovic, Todd J. Martinez, Philip H. Bucksbaum

Abstract: We have studied the photo-induced isomerization from 1,3-cyclohexadiene to 1,3,5-hexatriene in the presence of an intense ultrafast laser pulse. We find that the laser field maximally suppresses isomerization if it is both polarized parallel to the excitation dipole and present 50 fs after the initial photoabsorption, at the time when the system is expected to be in the vicinity of a conical inter… ▽ More We have studied the photo-induced isomerization from 1,3-cyclohexadiene to 1,3,5-hexatriene in the presence of an intense ultrafast laser pulse. We find that the laser field maximally suppresses isomerization if it is both polarized parallel to the excitation dipole and present 50 fs after the initial photoabsorption, at the time when the system is expected to be in the vicinity of a conical intersection that mediates this structural transition. A modified ab initio multiple spawning (AIMS) method shows that the laser induces a resonant coupling between the excited state and the ground state, i.e., a light-induced conical intersection. The theory accounts for the timing and direction of the effect. △ Less

Submitted 8 November, 2011; v1 submitted 26 September, 2011; originally announced September 2011.

Ultrafast Internal Conversion in Ethylene I. The Excited State Lifetime

Authors: H. Tao, T. K. Allison, T. W. Wright, A. M. Stooke, C. Khurmi, J. van Tilborg, Y. Liu, R. W. Falcone, A. Belkacem, T. J. Martinez

Abstract: Using a combined theoretical and experimental approach, we investigate the non-adiabatic dynamics of the prototypical ethylene (C2H4) molecule upon π \to π* excitation. In this first part of a two part series, we focus on the lifetime of the excited electronic state. The femtosecond Time-Resolved Photoelectron Spectrum (TRPES) of ethylene is simulated based on our recent molecular dynamics simulat… ▽ More Using a combined theoretical and experimental approach, we investigate the non-adiabatic dynamics of the prototypical ethylene (C2H4) molecule upon π \to π* excitation. In this first part of a two part series, we focus on the lifetime of the excited electronic state. The femtosecond Time-Resolved Photoelectron Spectrum (TRPES) of ethylene is simulated based on our recent molecular dynamics simulation using the ab initio multiple spawning method (AIMS) with Multi-State Second Order Perturbation Theory (Tao, et al. J. Phys. Chem. A 113 13656 2009). We find excellent agreement between the TRPES calculation and the photoion signal observed in a pump-probe experiment using femtosecond vacuum ultraviolet (hν = 7.7 eV) pulses for both pump and probe. These results explain the apparent discrepancy over the excited state lifetime between theory and experiment that has existed for ten years, with experiments (e.g., Farmanara, et al. Chem. Phys. Lett. 288 518 1998 and Kosma, et al. J. Phys. Chem. A 112 7514 2008) reporting much shorter lifetimes than predicted by theory. Investigation of the TRPES indicates that the fast decay of the photoion yield originates from both energetic and electronic factors, with the energetic factor playing a larger role in shaping the signal. △ Less

Submitted 3 October, 2011; v1 submitted 5 June, 2011; originally announced June 2011.

Journal ref: J. Chem. Phys, volume 134, page 244306 (2011)

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

Authors: Jiahao Chen, Todd J. Martínez

Abstract: 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… ▽ More 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. △ Less

Submitted 6 March, 2009; originally announced March 2009.

Comments: 9 pages, submitted to J. Chem. Phys

Journal ref: J. Chem. Phys. 131 044114 (2009)

Size-extensive polarizabilities with intermolecular charge transfer in a fluctuating-charge model

Authors: Jiahao Chen, Todd J Martinez

Abstract: Fluctuating-charge models have been used to model polarization effects in molecular mechanics methods. However, they overestimate polarizabilities in large systems. Previous attempts to remedy this have been at the expense of forbidding intermolecular charge-transfer. Here, we investigate this lack of size-extensivity and show that the neglect of terms arising from charge conservation is partly… ▽ More Fluctuating-charge models have been used to model polarization effects in molecular mechanics methods. However, they overestimate polarizabilities in large systems. Previous attempts to remedy this have been at the expense of forbidding intermolecular charge-transfer. Here, we investigate this lack of size-extensivity and show that the neglect of terms arising from charge conservation is partly responsible; these terms are also vital for maintaining the correct translational symmetries of the dipole moment and polarizability that classical electrostatic theory requires. Also, QTPIE demonstrates linear-scaling polarizabilities when coupling the external electric field in a manner that treats its potential as a perturbation of the atomic electronegativities. Thus for the first time, we have a fluctuating-charge model that predicts size-extensive dipole polarizabilities, yet allows intermolecular charge-transfer. △ Less

Submitted 8 December, 2008; originally announced December 2008.

Comments: 16 pages, 4 figures

The dissociation catastrophe in fluctuating-charge models and its implications for the concept of atomic electronegativity

Authors: Jiahao Chen, Todd J Martinez

Abstract: We have recently developed the QTPIE (charge transfer with polarization current equilibration) fluctuating-charge model, a new model with correct dissociation behavior for nonequilibrium geometries. The correct asymptotics originally came at the price of representing the solution in terms of charge-transfer variables instead of atomic charges. However, we have found an exact reformulation of flu… ▽ More We have recently developed the QTPIE (charge transfer with polarization current equilibration) fluctuating-charge model, a new model with correct dissociation behavior for nonequilibrium geometries. The correct asymptotics originally came at the price of representing the solution in terms of charge-transfer variables instead of atomic charges. However, we have found an exact reformulation of fluctuating-charge models in terms of atomic charges again, which is made possible by the symmetries of classical electrostatics. We show how this leads to the distinguishing between two types of atomic electronegativities in our model. While one is a intrinsic property of individual atoms, the other takes into account the local electrical surroundings. This suggests that this distinction could resolve some confusion surrounding the concept of electronegativity as to whether it is an intrinsic property of elements, or otherwise. △ Less

Submitted 8 December, 2008; originally announced December 2008.

Comments: 17 pages, prepared for "Proceedings of QSCP-XIII" in Prog. Theor. Chem. Phys

Journal ref: Prog. Theor. Chem. Phys., 19 (2009), 397-416

Multiple Spawning with Optimal Basis Set Expansion

Authors: Sandy Yang, Joshua D. Coe, Benjamin Kaduk, Todd J. Martinez

Abstract: The Full Multiple Spawning (FMS) method is designed to simulate quantum dynamics in the multi-state electronic problem. The FMS nuclear wavefunction is represented in a basis of coupled, frozen Gaussians, and the spawning procedure prescribes a means of adaptively increasing the size of the basis in order to capture population transfer between electronic states. Parent trajectories create childr… ▽ More The Full Multiple Spawning (FMS) method is designed to simulate quantum dynamics in the multi-state electronic problem. The FMS nuclear wavefunction is represented in a basis of coupled, frozen Gaussians, and the spawning procedure prescribes a means of adaptively increasing the size of the basis in order to capture population transfer between electronic states. Parent trajectories create children when passing through regions of significant nonadiabatic coupling. In order to converge branching ratios without allowing the basis to reach an impractical size, population transfer at individual spawning events should be made as effective as possible. Herein we detail a new algorithm for specifying the initial conditions of freshly spawned basis functions, one that minimizes the number of spawns needed for convergence by maximizing the efficiency of individual spawning events. Optimization is achieved by maximizing the coupling between parent and child trajectories, as a function of child position and momentum, at the point of spawning. The method is tested with a two-state, one-mode avoided crossing model and a two-state, two-mode conical intersection model. △ Less

Submitted 9 September, 2008; originally announced September 2008.

Comments: 35 pages, 8 figures, for submission to JCP

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

Authors: Jiahao Chen, Todd J. Martínez

Abstract: 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 elect… ▽ More 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. △ Less

Submitted 13 July, 2008; originally announced July 2008.

Comments: 36 pages, 7 figures; submitted to J. Chem. Phys

Journal ref: J. Chem. Phys. 129, 214113 (2008)

QTPIE: Charge Transfer with Polarization Current Equalization. A fluctuating charge model with correct asymptotics

Authors: Jiahao Chen, Todd J. Martínez

Abstract: Polarization and charge transfer are important effects which are difficult to describe using conventional force fields. Charge equilibration models can include both of these effects in large-scale molecular simulations. However, these models behave incorrectly when bonds are broken, making it difficult to use them in the context of reactive force fields. We develop a new method for describing ch… ▽ More Polarization and charge transfer are important effects which are difficult to describe using conventional force fields. Charge equilibration models can include both of these effects in large-scale molecular simulations. However, these models behave incorrectly when bonds are broken, making it difficult to use them in the context of reactive force fields. We develop a new method for describing charge flow in molecules -- QTPIE. The QTPIE method is based on charge transfer variables (as opposed to atomic charges) and correctly treats asymptotic behavior near dissociation. It is also able to provide a realistic description of in-plane polarizabilities. △ Less

Submitted 13 July, 2008; originally announced July 2008.

Comments: 22 pages, 5 figures

Journal ref: Chem. Phys. Lett. 438 (2007), 315-320