All Publications


  • Initial Conditions for Excited-State Dynamics in Solvated Systems: A Case Study. The journal of physical chemistry. B Curtis, E. R., Jones, C. M., Martínez, T. J. 2025

    Abstract

    Simulating excited-state dynamics or computing spectra for molecules in condensed phases requires sampling the ground state to generate initial conditions. Initial conditions (or snapshots for spectra) are typically produced by QM/MM Boltzmann sampling following MM equilibration or optimization. Given the switch from a MM to a QM/MM potential energy surface, one should discard a set period of time (which we call the "healing time") from the beginning of the QM/MM trajectory. Ideally, the healing time is as short as possible (to avoid unnecessary computational effort), but long enough to equilibrate to the QM/MM ground state distribution. Healing times in previous studies range from tens of femtoseconds to tens of picoseconds, suggesting the need for guidelines to choose a healing time. We examine the effect of healing time on the nonadiabatic dynamics and spectrum of a first-generation Donor-Acceptor Stenhouse Adduct in chloroform. Insufficient healing times skew the branching ratio of ground state products and alter the relaxation time for one pathway. The influence of the healing time on the absorption spectrum is less pronounced, warning that the spectrum is not a sensitive indicator for the quality of a set of initial conditions for dynamics. We demonstrate that a reasonable estimate for the healing time can be obtained by monitoring the solute temperature during the healing trajectory. We suggest that this procedure should become standard practice for determining healing times to generate initial conditions for nonadiabatic QM/MM simulations in large molecules and condensed phases.

    View details for DOI 10.1021/acs.jpcb.4c06536

    View details for PubMedID 39931914

  • Performance of Coupled-Cluster Singles and Doubles on Modern Stream Processing Architectures. Journal of chemical theory and computation Fales, B. S., Curtis, E. R., Johnson, K. G., Lahana, D., Seritan, S., Wang, Y., Weir, H., Martinez, T. J., Hohenstein, E. G. 2020

    Abstract

    We develop a new implementation of coupled-cluster singles and doubles (CCSD) optimized for the most recent graphical processing unit (GPU) hardware. We find that a single node with 8 NVIDIA V100 GPUs is capable of performing CCSD computations on roughly 100 atoms and 1300 basis functions in less than 1 day. Comparisons against massively parallel implementations of CCSD suggest that more than 64 CPU-based nodes (each with 16 cores) are required to match this performance.

    View details for DOI 10.1021/acs.jctc.0c00336

    View details for PubMedID 32567305