On the numerical solution of the exact factorization equations.

The exact factorization (EF) approach to coupled electron-ion dynamics recasts the time-dependent molecular Schrödinger equation as two coupled equations, one for the nuclear wavefunction and one for the conditional electronic wavefunction. The potentials appearing in these equations have provided insight into non-adiabatic processes, and new practical non-adiabatic dynamics methods have been formulated starting from these equations. Here, we provide a first demonstration of a self-consistent solution of the exact equations, with a preliminary analysis of their stability and convergence properties. The equations have an unprecedented mathematical form, involving a Hamiltonian outside the class of Hermitian Hamiltonians usually encountered in time-propagation, and so the usual numerical methods for time-dependent Schrödinger fail when applied in a straightforward way to the EF equations. We find an approach that enables stable propagation long enough to witness non-adiabatic behavior in a model system before non-trivial instabilities take over. Implications for the development and analysis of EF-based methods are discussed.

Suppression of Th1 Priming by TLR2 Agonists during Cutaneous Immunization Is Mediated by Recruited CCR2+ Monocytes.

Effective subunit vaccines require the incorporation of adjuvants that stimulate cells of the innate immune system to generate protective adaptive immune responses. Pattern recognition receptor agonists are a growing class of potential adjuvants that can shape the character of the immune response to subunit vaccines by directing the polarization of CD4 T cell differentiation to various functional subsets. In the current study, we applied a high-throughput in vitro screen to assess murine CD4 T cell polarization by a panel of pattern recognition receptor agonists. This identified lipopeptides with TLR2 agonist activity as exceptional Th1-polarizing adjuvants. In vivo, we demonstrated that i.v. administration of TLR2 agonists with Ag in mice replicated the findings from in vitro screening by promoting strong Th1 polarization. In contrast, TLR2 agonists inhibited priming of Th1 responses when administered cutaneously in mice. This route-specific suppression was associated with infiltrating CCR2+ cells in the skin-draining lymph nodes and was not uniquely dependent on any of the well characterized subsets of dendritic cells known to reside in the skin. We further demonstrated that priming of CD4 T cells to generate Th1 effectors following immunization with the Mycobacterium bovis bacillus Calmette-Guérin (BCG) strain, a lipoprotein-rich bacterium recognized by TLR2, was dependent on the immunization route, with significantly greater Th1 responses with i.v. compared with intradermal administration of BCG. A more complete understanding of route-dependent TLR2 responses may be critical for informed design of novel subunit vaccines and for improvement of BCG and other vaccines based on live-attenuated organisms.

Coupled-Trajectory Mixed Quantum-Classical Algorithm: A Deconstruction.

We analyze a mixed quantum-classical algorithm recently derived from the exact factorization equations [Min, Agostini, Gross, PRL 115, 073001 (2015)] to show the role of the different terms in the algorithm in bringing about decoherence and wavepacket branching. The algorithm has the structure of Ehrenfest equations plus a "coupled-trajectory" term for both the electronic and nuclear equations, and we analyze the relative roles played by the different nonadiabatic terms in these equations, including how they are computed in practice. In particular, we show that while the coupled-trajectory term in the electronic equation is essential in yielding accurate dynamics, that in the nuclear equation has a much smaller effect. A decoherence time is extracted from the electronic equations and compared with that of augmented fewest-switches surface-hopping. We revisit a series of nonadiabatic Tully model systems to illustrate our analysis.