Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials.

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes.

The photolytic stability of explosives and energetic functional groups is of importance for those who regularly handle or are exposed to explosives in typical environmental conditions. This study models the photolytic degradation of dodecane substituted with various energetic functional groups: azide, nitro, nitrate ester, and nitramine. For the studied molecules, it was found that excitons localize on the energetic functional group, no matter where they were initially formed, and thus, the predominant degradation pathway involves the degradation of the energetic functional group. The relative trends for both 4 and 8 eV excitation energies followed with what is expected from the relative stability of the energetic functional groups to thermal and sub-shock degradation. The one notable exception was the azide functional group; more work should be done to further understand the photolytic effects on the azide functional group.

Sexually Transmitted Infections in Association With Area-Level Prostitution and Drug-Related Arrests.

OBJECTIVES: This study aimed to examine the mediators and moderators of area-level prostitution arrests and sexually transmitted infections (STIs) using population-level data. METHODS: Using justice and public health STI/HIV data in Marion County (Indianapolis), Indiana, during an 18-year period, we assessed the overall association of area-level prostitution and drug-related arrests and STI/HIV, and mediators and moderators of the relationship. Point-level arrests were geocoded and aggregated by a census block group. RESULTS: Results indicate a positive relationship between numbers of prostitution arrests and area-level STI rates. There was a dose-response relationship between prostitution arrests and STI rates when accounting for drug-related arrests. The highest quintile block groups had significantly higher rates of reported chlamydia (incident rate ratio [IRR], 3.29; 95% confidence interval [CI], 2.82-3.84), gonorrhea (IRR, 4.73; 95% CI, 3.90-5.57), syphilis (IRR, 4.28; 95% CI, 3.47-5.29), and HIV (IRR, 2.76; 95% CI, 2.24-3.39) compared with the lowest quintile. When including drug arrests, the second (IRR, 1.19; 95% CI, 1.03-1.38) and the third (IRR, 1.20; 95% CI, 1.02-1.41) highest quintile block groups had lower IRR for reported rates of chlamydia, indicating that drug arrests mediated the prostitution arrest effect. CONCLUSIONS: These findings inform public health agencies and community-based organizations that conduct outreach in these areas to expand their efforts to include harm reduction and HIV/STI testing for both sex workers and individuals experiencing substance use disorder. Another implication of these data is the importance of greater collaboration in public health and policing efforts to address overlapping epidemics that engage both health and legal interventions.

Excitation Energy Transfer between bodipy Dyes in a Symmetric Molecular Excitonic Seesaw.

We examine the redistribution of energy between electronic and vibrational degrees of freedom that takes place between a π-conjugated oligomer, a phenylene-butadiynylene, and two identical boron-dipyrromethene (bodipy) end-caps using femtosecond transient absorption spectroscopy, single-molecule spectroscopy, and nonadiabatic excited-state molecular dynamics (NEXMD) modeling techniques. The molecular structure represents an excitonic seesaw in that the excitation energy on the oligomer backbone can migrate to either one end-cap or the other, but not to both. The NEXMD simulations closely reproduce the characteristic time scale for redistribution of electronic and vibrational energy of 2.2 ps and uncover the vibrational modes contributing to the intramolecular relaxation. The calculations indicate that the dihedral angle between the bodipy dye and the oligomer change upon excitation of the oligomer. Single-molecule experiments reveal a difference in photoluminescence lifetime of the bodipy dyes depending on whether they are excited by direct absorption or by redistribution of energy from the backbone. This difference in lifetime may be attributed to the difference in dihedral angle. The simulations also suggest that a strong coupling can occur between the two end-caps, giving rise to a reversible shuttling of excitation energy between them. Strong coupling should lead to a pronounced loss in polarization memory of the fluorescence since the oligomer backbone tends to be slightly distorted and the two bodipy transition dipoles have different orientations. A sensitive single-molecule technique is presented to test for such coupling. However, although redistribution of electronic and vibrational energy between the end-caps can occur, it appears to be unidirectional and irreversible, suggesting that an additional localization mechanism is at play which is, as yet, not fully accounted for in the simulations.

Correlating Mechanical Sensitivity with Spin Transition in the Explosive Spin Crossover Complex [Fe(Htrz)3]n[ClO4]2n.

Spin crossover complexes are known to undergo bond length, volume, and enthalpy changes during spin transition. In an explosive spin crossover complex, these changes could affect the mechanical and initiation sensitivity of the explosive and lead to the development of a new class of sensitivity switchable materials. To explore this relationship, the well-known spin crossover compound [Fe(Htrz)3]n[ClO4]2n (1) was re-evaluated for its explosive properties, and its mechanical impact sensitivity was correlated to spin transition. A variable temperature impact test was developed and used to evaluate the impact sensitivity of 1 in the low spin (LS, S = 0), thermally accessed high spin (HS, S = 2), and mixed LS and HS states. For comparison, the structurally similar Ni compound, [Ni(Htrz)3]n[ClO4]2n (2), which does not undergo a spin transition at accessible temperatures, was synthesized and characterized, and its explosive properties and variable temperature impact sensitivity measured. These results reveal a correlation between impact sensitivity and spin transition, where 1 exhibits lower impact sensitivity in the LS state and increases in sensitivity upon transition to the HS state. Density functional theory was used to predict structural changes that occur upon spin transition that correlate to the change in sensitivity. This demonstrates, for the first time, an explosive spin crossover compound (ExSCO) that exhibits switchable impact sensitivity with a fully reversible internal switching mechanism.

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born-Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor-acceptor-donor (D-A-D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are "sufficiently good" to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

Non-adiabatic molecular dynamics of molecules in the presence of strong light-matter interactions.

When the interaction between a molecular system and confined light modes in an optical or plasmonic cavity is strong enough to overcome the dissipative process, hybrid light-matter states (polaritons) become the fundamental excitations in the system. The mixing between the light and matter characters modifies the photophysical and photochemical properties. Notably, it was reported that these polaritons can be employed to control photochemical reactions, charge and energy transfer, and other processes. In addition, according to recent studies, vibrational strong coupling can be employed to resonantly enhance the thermally-activated chemical reactions. In this work, a theoretical model and an efficient numerical method for studying the dynamics of molecules strongly interacting with quantum light are developed based on nonadiabatic excited-state molecular dynamics. The methodology was employed to study the cis-trans photoisomerization of a realistic molecule in a cavity. Numerical simulations demonstrate that the photochemical reactions can be controlled by tuning the properties of the cavity. In the calculated example, the isomerization is suppressed when polaritonic states develop a local minimum on the lower polaritonic state. Moreover, the observed reduction of isomerization is tunable via the photon energy and light-molecule coupling strength. However, the fluctuation in the transition dipole screens the effect of light-matter, which makes it harder to tune the photochemical properties via the coupling strength. These insights suggest quantum control of photochemical reactions is possible by specially designed photonic or plasmonic cavities.

Site-Specific Photodecomposition in Conjugated Energetic Materials.

Nonadiabatic excited-state molecular dynamics (NEXMD) has been used to study photodecomposition in a class of recently synthesized bicyclic conjugated energetic materials (CEMs) composed of fused tetrazole and tetrazine derivatives with increasing oxygen substitutions. Modification by oxygen functionalization has already been demonstrated to increase the two-photon absorption intensity in the target CEMs while simultaneously improving oxygen balance. Photodecomposition mechanisms in materials that undergo nonlinear absorption could be used to achieve controlled, direct optical initiation. Here, we use NEXMD simulations to model the nonradiative relaxation and photodecomposition in CEMs following photoexcitation by a simulated Nd:YAG laser pulse. Excess electronic energy is quickly converted into vibrational energy on a sub-100 fs time scale resulting in bond dissociation. We find that, for the studied tetrazine derivatives, the bicyclic framework is an important structural feature that enhances the photochemical quantum yield and the high atomic oxygen content increases the relaxation lifetime and opens additional photodissociation pathways targeting the oxygen-substituted sites. The presented analysis scheme based on bond orders in the swarm of NEXMD trajectories is a useful tool for determining photochemical reactions.

NEXMD Modeling of Photoisomerization Dynamics of 4-Styrylquinoline.

Isomerization of molecular systems is ubiquitous in chemistry and biology, and is also important for many applications. Atomistic simulations can help determine the tunable parameters influencing this process. In this paper, we use the Nonadiabatic EXcited state Molecular Dynamics (NEXMD) software to study the photoisomerization of a representative molecule, 4-styrylquinoline (SQ). trans-SQ transforms into dihydrobenzophenanthridine (DHBP) upon irradiation with laser light, with the cis conformer acting as an intermediate. We study how varying three different external stimuli (i.e., apolar versus polar solvent, low versus high photoexcitation energy, and vacuum versus a constant temperature thermostat) affects the trans-to- cis photoisomerization of SQ. Our results show that polarization effects due to implicit solvation and the thermostat play a crucial role in the isomerization process, whereas photoexcitation energy plays a lesser role on the outcome and efficiency. We also show that NEXMD captures the correct energy profile between the ground and first singlet excited state, showing that there are two distinct reaction pathways to the final stable product that vary by the number of photons absorbed, in agreement with experiment. Ultimately, NEXMD proves to be an effective tool for investigating excited state single molecule dynamics subject to various environments and initial conditions.

Phonon bottleneck and long-lived excited states in π-conjugated pyrene hoop.

In the last decade, recent synthetic advances have launched carbon-based π-conjugated hoops to the forefront of theoretical and experimental investigation not only for their potential use as bottom-up templates for carbon nanotube (CNT) growth, but also for the interesting excitonic effects arising from the cyclic geometry, unique π-system orientation, and unusual electronic interactions and couplings. In particular, cyclic materials based on pyrene, a common component in organic electronics, are popular candidates for the future design of π-conjugated nanorings for optoelectronic applications. Understanding the photophysical response in cyclic oligopyrenes can be achieved using non-adiabatic excited state molecular dynamics (NA-ESMD). Through NA-ESMD modeling, we reveal details of the nonradiative relaxation processes in the circular pyrene tetramer [4]cyclo-2,7-pyrenylene ([4]CPY) where we find that the strong non-adiabatic coupling combined with the dense manifold of excited states creates an internal conversion mechanism dominated by ultrafast sequential quantum transitions. However, we observe two long-lived electronic excited states that introduce a phonon bottleneck in the electronic relaxation process. In fact, the timescale for the electronic relaxation is almost exclusively dominated by the lifetimes of the long-lived states. We find that the states associated with the phonon bottleneck are separated from lower energy states by large energy gaps and are characterized by localization on a single pyrene unit resulting in a spatial mismatch with strongly delocalized neighboring states.

Ultrafast Photodissociation Dynamics of Nitromethane.

Nitromethane (NM), a high explosive (HE) with low sensitivity, is known to undergo photolysis upon ultraviolet (UV) irradiation. The optical transparency, homogeneity, and extensive study of NM make it an ideal system for studying photodissociation mechanisms in conventional HE materials. The photochemical processes involved in the decomposition of NM could be applied to the future design of controllable photoactive HE materials. In this study, the photodecomposition of NM from the nπ* state excited at 266 nm is being investigated on the femtosecond time scale. UV femtosecond transient absorption (TA) spectroscopy and excited state femtosecond stimulated Raman spectroscopy (FSRS) are combined with nonadiabatic excited state molecular dynamics (NA-ESMD) simulations to provide a unified picture of NM photodecomposition. The FSRS spectrum of the photoproduct exhibits peaks in the NO2 region and slightly shifted C-N vibrational peaks pointing to methyl nitrite formation as the dominant photoproduct. A total photolysis quantum yield of 0.27 and an nπ* state lifetime of ∼20 fs were predicted from NA-ESMD simulations. Predicted time scales revealed that NO2 dissociation occurs in 81 ± 4 fs and methyl nitrite formation is much slower having a time scale of 452 ± 9 fs corresponding to the excited state absorption feature with a decay of 480 ± 17 fs observed in the TA spectrum. Although simulations predict C-N bond cleavage as the primary photochemical process, the relative time scales are consistent with isomerization occurring via NO2 dissociation and subsequent rebinding of the methyl radical and nitrogen dioxide.

Two-Photon Absorption in Conjugated Energetic Molecules.

Time-dependent density functional theory (TD-DFT) was used to investigate the relationship between molecular structure and the one- and two-photon absorption (OPA and TPA, respectively) properties of novel and recently synthesized conjugated energetic molecules (CEMs). The molecular structures of CEMs can be strategically altered to influence the heat of formation and oxygen balance, two factors that can contribute to the sensitivity and strength of an explosive material. OPA and TPA are sensitive to changes in molecular structure as well, influencing the optical range of excitation. We found calculated vertical excitation energies to be in good agreement with experiment for most molecules. Peak TPA intensities were found to be significant and on the order of 10(2) GM. Natural transition orbitals for essential electronic states defining TPA peaks of relatively large intensity were used to examine the character of relevant transitions. Modification of molecular substituents, such as additional oxygen or other functional groups, produces significant changes in electronic structure, OPA, and TPA and improves oxygen balance. The results show that certain molecules are apt to undergo nonlinear absorption, opening the possibility for controlled, direct optical initiation of CEMs through photochemical pathways.

Nonadiabatic excited-state molecular dynamics: modeling photophysics in organic conjugated materials.

To design functional photoactive materials for a variety of technological applications, researchers need to understand their electronic properties in detail and have ways to control their photoinduced pathways. When excited by photons of light, organic conjugated materials (OCMs) show dynamics that are often characterized by large nonadiabatic (NA) couplings between multiple excited states through a breakdown of the Born-Oppenheimer (BO) approximation. Following photoexcitation, various nonradiative intraband relaxation pathways can lead to a number of complex processes. Therefore, computational simulation of nonadiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons. Over the years, we have developed a nonadiabatic excited-state molecular dynamics (NA-ESMD) framework that efficiently and accurately describes photoinduced phenomena in extended conjugated molecular systems. We use the fewest-switches surface hopping (FSSH) algorithm to treat quantum transitions among multiple adiabatic excited state potential energy surfaces (PESs). Extended molecular systems often contain hundreds of atoms and involve large densities of excited states that participate in the photoinduced dynamics. We can achieve an accurate description of the multiple excited states using the configuration interaction single (CIS) formalism with a semiempirical model Hamiltonian. Analytical techniques allow the trajectory to be propagated "on the fly" using the complete set of NA coupling terms and remove computational bottlenecks in the evaluation of excited-state gradients and NA couplings. Furthermore, the use of state-specific gradients for propagation of nuclei on the native excited-state PES eliminates the need for simplifications such as the classical path approximation (CPA), which only uses ground-state gradients. Thus, the NA-ESMD methodology offers a computationally tractable route for simulating hundreds of atoms on ~10 ps time scales where multiple coupled excited states are involved. In this Account, we review recent developments in the NA-ESMD modeling of photoinduced dynamics in extended conjugated molecules involving multiple coupled electronic states. We have successfully applied the outlined NA-ESMD framework to study ultrafast conformational planarization in polyfluorenes where the rate of torsional relaxation can be controlled based on the initial excitation. With the addition of the state reassignment algorithm to identify instances of unavoided crossings between noninteracting PESs, NA-ESMD can now be used to study systems in which these so-called trivial unavoided crossings are expected to predominate. We employ this technique to analyze the energy transfer between poly(phenylene vinylene) (PPV) segments where conformational fluctuations give rise to numerous instances of unavoided crossings leading to multiple pathways and complex energy transfer dynamics that cannot be described using a simple Förster model. In addition, we have investigated the mechanism of ultrafast unidirectional energy transfer in dendrimers composed of poly(phenylene ethynylene) (PPE) chromophores and have demonstrated that differential nuclear motion favors downhill energy transfer in dendrimers. The use of native excited-state gradients allows us to observe this feature.

Photoactive high explosives: linear and nonlinear photochemistry of petrin tetrazine chloride.

Pentaerythritol tetranitrate (PETN), a high explosive, initiates with traditional shock and thermal mechanisms. In this study, the tetrazine-substituted derivative of PETN, pentaerythritol trinitrate chlorotetrazine (PetrinTzCl), is being investigated for a photochemical initiation mechanism that could allow control over the chemistry contributing to decomposition leading to initiation. PetrinTzCl exhibits a photochemical quantum yield (QYPC) at 532 nm not evident with PETN. Using static spectroscopic methods, we observe energy absorption on the tetrazine (Tz) ring that results in photodissociation yielding N2, Cl-CN, and Petrin-CN as the major photoproducts. The QYPC was enhanced with increasing irradiation intensity. Experiment and theoretical calculations imply this excitation mechanism follows sequential photon absorption. Dynamic simulations demonstrate that the relaxation mechanism leading to the observed photochemistry in PetrinTzCl is due to vibrational excitation during internal conversion. PetrinTzCl's single photon stability and intensity dependence suggest this material could be stable in ambient lighting, yet possible to initiate with short-pulsed lasers.

Signature of nonadiabatic coupling in excited-state vibrational modes.

Using analytical excited-state gradients, vibrational normal modes have been calculated at the minimum of the electronic excited-state potential energy surfaces for a set of extended conjugated molecules with different coupling between them. Molecular model systems composed of units of polyphenylene ethynylene (PPE), polyphenylenevinylene (PPV), and naphthacene/pentacene (NP) have been considered. In all cases except the NP model, the influence of the nonadiabatic coupling on the excited-state equilibrium normal modes is revealed as a unique highest frequency adiabatic vibrational mode that overlaps with the coupling vector. This feature is removed by using a locally diabatic representation in which the effect of NA interaction is removed. Comparison of the original adiabatic modes with a set of vibrational modes computed in the locally diabatic representation demonstrates that the effect of nonadiabaticity is confined to only a few modes. This suggests that the nonadiabatic character of a molecular system may be detected spectroscopically by identifying these unique state-specific high frequency vibrational modes.

Obesity Prevention in Early Life (OPEL) study: linking longitudinal data to capture obesity risk in the first 1000 days.

To develop robust prediction models for infant obesity risk, we need data spanning multiple levels of influence, including child clinical health outcomes (eg, height and weight), information about maternal pregnancy history, detailed sociodemographic information of parents and community-level factors. Few data sources contain all of this information. This manuscript describes the creation of the Obesity Prevention in Early Life (OPEL) database, a longitudinal, population-based database that links clinical data with birth certificates and geocoded area-level indicators for 19 437 children born in Marion County, Indiana between 2004 and 2019. This brief describes the methodology of linking administrative data, the establishment of the OPEL database, and the clinical and public health implications facilitated by these data. The OPEL database provides a strong basis for further longitudinal child health outcomes studies and supports the continued development of intergenerational linked clinical-public health databases.

Extremely long nonradiative relaxation of photoexcited graphane is greatly accelerated by oxidation: time-domain ab initio study.

Graphane and its derivatives are stable and extremely thin, wide band gap semiconductors that promise to replace conventional semiconductors in electronics, catalysis, and energy applications, greatly reducing device size and power consumption. In order to be useful, band-gap excitations in these materials should be long lived and nonradiative energy losses to heat should be slow. We use state-of-the-art nonadiabatic molecular dynamics combined with time-dependent density functional theory in order to determine the nonradiative lifetime and radiative line width of the lowest energy singlet excitations in pure and oxidized graphanes. We predict that pure graphane has a very long nonradiative decay time, on the order of 100 ns, while epoxy- and hydroxy-graphanes lose electronic excitation energy to heat 10-20 times faster. The luminescence line width is 1.5 times larger in pristine graphane compared to its oxidized forms, and at room temperature, it is on the order of 50 meV. Hydroxylation lowers graphane's band gap, while epoxidation increases the gap. The nonradiative decay and luminescence line width of pure graphane are governed by electron coupling to the 1200 cm(-1) vibrational mode. In the oxidized forms of graphane, the electronic excitations couple to a broad range of vibrational modes, rationalizing the more rapid nonradiative decay in these systems. The slow electron-phonon energy losses in graphane compared to other graphene derivatives, such as carbon nanotubes and nanoribbons, indicate that graphanes are excellent candidates for semiconductor applications.

Conformational disorder in energy transfer: beyond Förster theory.

Energy transfer in donor-acceptor chromophore pairs, where the absorption of each species is well separated while donor emission and acceptor absorption overlap, can be understood through a Förster resonance energy transfer model. The picture is more complex for organic conjugated polymers, where the total absorption spectrum can be described as a sum of the individual contributions from each subunit (chromophore), whose absorption is not well separated. Although excitations in these systems tend to be well localized, traditional donors and acceptors cannot be defined and energy transfer can occur through various pathways where each subunit (chromophore) is capable of playing either role. In addition, fast torsional motions between individual monomers can break conjugation and lead to reordering of excited state energy levels. Fast torsional fluctuations occur on the same timescale as electronic transitions leading to multiple trivial unavoided crossings between excited states during dynamics. We use the non-adiabatic excited state molecular dynamics (NA-ESMD) approach to simulate energy transfer between two poly-phenylene vinylene (PPV) oligomers composed of 3-rings and 4-rings, respectively, separated by varying distances. The change in the spatial localization of the transient electronic transition density, initially localized on the donors, is used to determine the transfer rate. Our analysis shows that evolution of the intramolecular transition density can be decomposed into contributions from multiple transfer pathways. Here we present a detailed analysis of ensemble dynamics as well as a few representative trajectories which demonstrate the intertwined role of electronic and conformational processes. Our study reveals the complex nature of energy transfer in organic conjugated polymer systems and emphasizes the caution that must be taken in performing such an analysis when a single simple unidirectional pathway is unlikely.

Nonadiabatic excited-state molecular dynamics: treatment of electronic decoherence.

Within the fewest switches surface hopping (FSSH) formulation, a swarm of independent trajectories is propagated and the equations of motion for the quantum coefficients are evolved coherently along each independent nuclear trajectory. That is, the phase factors, or quantum amplitudes, are retained. At a region of strong coupling, a trajectory can branch into multiple wavepackets. Directly following a hop, the two wavepackets remain in a region of nonadiabatic coupling and continue exchanging population. After these wavepackets have sufficiently separated in phase space, they should begin to evolve independently from one another, the process known as decoherence. Decoherence is not accounted for in the standard surface hopping algorithm and leads to internal inconsistency. FSSH is designed to ensure that at any time, the fraction of classical trajectories evolving on each quantum state is equal to the average quantum probability for that state. However, in many systems this internal consistency requirement is violated. Treating decoherence is an inherent problem that can be addressed by implementing some form of decoherence correction to the standard FSSH algorithm. In this study, we have implemented two forms of the instantaneous decoherence procedure where coefficients are reinitialized following hops. We also test the energy-based decoherence correction (EDC) scheme proposed by Granucci et al. and a related version where the form of the decoherence time is taken from Truhlar's Coherent Switching with Decay of Mixing method. The sensitivity of the EDC results to changes in parameters is also evaluated. The application of these computationally inexpensive ad hoc methods is demonstrated in the simulation of nonradiative relaxation in two conjugated oligomer systems, specifically poly-phenylene vinylene and poly-phenylene ethynylene. We find that methods that have been used successfully for treating small systems do not necessarily translate to large polyatomic systems and their success depends on the particular system under study.

HIV Care Continuum Among People Living With HIV and History of Arrest and Mental Health Diagnosis.

OBJECTIVE: Justice involvement and psychiatric comorbidities contribute to excess HIV morbidity, yet their interaction is poorly understood. We examined associations of this overlap with HIV outcomes among people living with HIV (PLWH). METHODS: We conducted a retrospective cohort study of PLWH aged 13 years and older residing in Marion County (Indianapolis), IN, during 2018 (n = 5730) using linked HIV surveillance, arrest, and clinical data. We used univariable and multivariable regression to evaluate main and interaction effects of 2010-2017 arrest and mental health diagnosis on 2018 linkage to care (LTC), retention in care (RIC), and undetectable viral load (UVL). RESULTS: LTC decreased among those with, versus without, an arrest (P = 0.02), although mental health diagnoses had no significant effect on LTC. When controlling for demographics and substance use disorder, analyses indicated a protective effect of arrest history on odds of RIC (adjusted odds ratio [aOR] = 1.54) and UVL (aOR = 1.26). Mental health diagnosis also increased odds of RIC (aOR = 2.02) and UVL (aOR = 1.95). Post hoc tests demonstrated that these results were mediated by outpatient care utilization, although an arrest or mental health diagnosis did increase odds of RIC among PLWH and a history of low outpatient utilization. CONCLUSIONS: Outpatient care utilization improves HIV outcomes, even among those with justice involvement and psychiatric comorbidities. Holistic approaches to care can increase utilization. Implementation of "no wrong door" approaches, such as integration of mental health care in the primary care setting, simplifies health care navigation and improves access. Among those arrested, access to a Behavioral Court program can improve, rather than disrupt, HIV care.