An official website of the United States government
We present a methodology for simulating multidimensional electronic spectra of molecular aggregates with coupling of electronic excitation to a structured environment using the stochastic non-Markovian quantum state diffusion (NMQSD) method in combination with perturbation theory for the response functions. A crucial aspect of our approach is that we propagate the NMQSD equation in a doubled system Hilbert space but with the same noise. We demonstrate that our approach shows fast convergence with respect to the number of stochastic trajectories, providing a promising technique for numerical calculation of two-dimensional electronic spectra of large molecular aggregates.
more »
« less
Bethe–Salpeter equation spectra for very large systems
We present a highly efficient method for the extraction of optical properties of very large molecules via the Bethe–Salpeter equation. The crutch of this approach is the calculation of the action of the effective Coulombic interaction, W, through a stochastic time-dependent Hartree propagation, which uses only ten stochastic orbitals rather than propagating the full sea of occupied states. This leads to a scaling that is at most cubic in system size with trivial parallelization of the calculation. We apply this new method to calculate the spectra and electronic density of the dominant excitons of a carbon-nanohoop bound fullerene system with 520 electrons using less than 4000 core hours.
more »
« less
- Award ID(s):
- 1763176
- PAR ID:
- 10425677
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 157
- Issue:
- 3
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 031104
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Corminboeuf, Clémence; Gagliardi, Laura (Ed.)We develop an efficient approach to evaluate range-separated exact exchange for grid- or plane-wave-based representations within the generalized Kohn–Sham–density functional theory (GKS–DFT) framework. The Coulomb kernel is fragmented in reciprocal space, and we employ a mixed deterministic-stochastic representation, retaining long-wavelength (low-k) contributions deterministically and using a sparse (“fragmented”) stochastic basis for the high-k part. Coupled with a projection of the Hamiltonian onto a subspace of valence and conduction states from a prior local-DFT calculation, this method allows for the calculation of the long-range exchange of large molecular systems with hundreds and potentially thousands of coupled valence states delocalized over millions of grid points. We find that even a small number of valence and conduction states is sufficient for converging the HOMO and LUMO energies of the GKS–DFT. Excellent tuning of long-range separated hybrids (RSH) is easily obtained in the method for very large systems, as exemplified here for the chlorophyll hexamer of Photosystem II with 1320 electrons.more » « less
-
The seismic fragility of a system is the probability that the system enters a damage state under seismic ground motions with specified characteristics. Plots of the seismic fragilities with respect to scalar ground motion intensity measures are called fragility curves. Recent studies show that fragility curves may not be satisfactory measures for structural seismic performance, since scalar intensity measures cannot comprehensively characterize site seismicity. The limitations of traditional seismic intensity measures, e.g., peak ground acceleration or pseudo-spectral acceleration, are shown and discussed in detail. A bivariate vector with coordinates moment magnitude m and source-to-site distance r is proposed as an alternative seismic intensity measure. Implicitly, fragility surfaces in the (m, r)-space could be used as graphical representations of seismic fragility. Unlike fragility curves, which are functions of scalar intensity measures, fragility surfaces are characterized by two earthquake-hazard parameters, (m, r). The calculation of fragility surfaces may be computationally expensive for complex systems. Thus, as solutions to this issue, a bi-variate log-normal parametric model and an efficient calculation method, based on stochastic-reduced-order models, for fragility surfaces are proposed.more » « less
-
We propose a free energy calculation method for receptor–ligand binding, which have multiple binding poses that avoids exhaustive enumeration of the poses. For systems with multiple binding poses, the standard procedure is to enumerate orientations of the binding poses, restrain the ligand to each orientation, and then, calculate the binding free energies for each binding pose. In this study, we modify a part of the thermodynamic cycle in order to sample a broader conformational space of the ligand in the binding site. This modification leads to more accurate free energy calculation without performing separate free energy simulations for each binding pose. We applied our modification to simple model host–guest systems as a test, which have only two binding poses, by using a single decoupling method (SDM) in implicit solvent. The results showed that the binding free energies obtained from our method without knowing the two binding poses were in good agreement with the benchmark results obtained by explicit enumeration of the binding poses. Our method is applicable to other alchemical binding free energy calculation methods such as the double decoupling method (DDM) in explicit solvent. We performed a calculation for a protein–ligand system with explicit solvent using our modified thermodynamic path. The results of the free energy simulation along our modified path were in good agreement with the results of conventional DDM, which requires a separate binding free energy calculation for each of the binding poses of the example of phenol binding to T4 lysozyme in explicit solvent. © 2019 Wiley Periodicals, Inc.more » « less
-
We present a scalable methodology to verify stochastic hybrid systems for inequality linear temporal logic (iLTL) or inequality metric interval temporal logic (iMITL). Using the Mori–Zwanzig reduction method, we construct a finite-state Markov chain reduction of a given stochastic hybrid system and prove that this reduced Markov chain is approximately equivalent to the original system in a distributional sense. Approximate equivalence of the stochastic hybrid system and its Markov chain reduction means that analyzing the Markov chain with respect to a suitably strengthened property allows us to conclude whether the original stochastic hybrid system meets its temporal logic specifications. Based on this, we propose the first statistical model checking algorithms to verify stochastic hybrid systems against correctness properties, expressed in iLTL or iMITL. The scalability of the proposed algorithms is demonstrated by a case study.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0100213
