Proton affinity and gas phase basicity of diamandoid molecules: diamantane to C131H116.

Calculated proton affinities (PAs) and gas phase basicities (GPBs) are reported for diamantane (C14H20), triamantane (C18H24), 'globular and planar' isomers of tetramantane (C22H28) and pentamantane (C26H32), and for one 'globular' isomer of each of the larger diamondoid molecules: C51H58, C78H72, C102H90, and C131H116. Assuming CxHy as the parent diamondoid molecule, we calculated PA and GPB values for a variety of CxHy+1+ isomers, as well as for the reaction CxHy + H+ yielding CxHy-1+ + H2(g); the latter is slightly favored based on GPB values for diamantane through pentamantane, but less favored compared to certain CxHy+1+ isomers of C51H58, C102H90, and C131H116. Indeed, the GPB values of C51H58, C102H90, and C131H116 classifiy them as 'superbases'. Calculations that had the initial location of the proton in an interstitial site inside the diamondoid molecule always showed the H having moved to the outside of the diamondoid molecule; for this reason, we focused on testing a variety of initial configurations with the proton placed in an initial position on the surface. Additional protons were added to determine the limiting number that could be, per these calculations, taken up by the diamondoid molecules and the maximum number of protons are shown in parentheses: C14H20(2), C18H24(3), C22H28(3), C26H32(3), C51H58(4). Bader charge distributions obtained for CxHy+1+ isomers (for diamantane through pentamantane) suggest that the positive charge is essentially completely delocalized over all the H atoms. NMR spectra were calculated for different isomers of C14H19+, and compared to the published NMR spectrum for when diamantane was mixed with magic acid and H2(g) was produced.

Independent trajectory mixed quantum-classical approaches based on the exact factorization.

Mixed quantum-classical dynamics based on the exact factorization exploits the "derived" electron-nuclear correlation (ENC) term, aiming for the description of quantum coherences. The ENC contains interactions between the phase of electronic states and nuclear quantum momenta, which depend on the spatial shape of the nuclear density. The original surface hopping based on the exact factorization (SHXF) [Ha et al., J. Phys. Chem. Lett. 9, 1097 (2018)] exploits frozen Gaussian functions to construct the nuclear density in the ENC term, while the phase of electronic states is approximated as a fictitious nuclear momentum change. However, in reality, the width of nuclear wave packets varies in time depending on the shape of potential energy surfaces. In this work, we present a modified SHXF approach and a newly developed Ehrenfest dynamics based on the exact factorization (EhXF) with time-dependent Gaussian functions and phases by enforcing total energy conservation. We perform numerical tests for various one-dimensional two-state model Hamiltonians. Overall, the time-dependent width of Gaussian functions and the energy conserving phase show a reliable decoherence compared to the original frozen Gaussian-based SHXF and the exact quantum mechanical calculation. In particular, the energy conserving phase is crucial for EhXF to reproduce the correct quantum dynamics.

PyUNIxMD: A Python-based excited state molecular dynamics package.

Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.

Real-Space and Real-Time Propagation for Correlated Electron-Nuclear Dynamics Based on Exact Factorization.

We present coupled equations of motion for correlated electron-nuclear dynamics for real-space and real-time propagation with a proper electron-nuclear correlation (ENC) from the exact factorization. Since the original ENC term from the exact factorization is non-Hermitian, the numerical instability arises as we propagate an electronic wave function. In this paper, we propose a Hermitian-type ENC term which depends on the electron density matrix and the nuclear quantum momentum. Moreover, we show that the Hermitian property of the electron-nuclear correlation term can capture quantum (de)coherence with a stable numerical real-space and real-time propagation. As an application, we demonstrate a real-space and real-time propagation of an electronic wave function coupled to trajectory-based nuclear motion for a one-dimensional model Hamiltonian. Our approach can capture nonadiabatic phenomena as well as quantum decoherence in excited state molecular dynamics. In addition, we propose a scheme to extend the current approach to many-body electronic states based on real-time time-dependent density functional theory, testing the nonadiabatic dynamics of a simple molecular system.

Coupled- and Independent-Trajectory Approaches Based on the Exact Factorization Using the PyUNIxMD Package.

We present mixed quantum-classical approaches based on the exact factorization framework. The electron-nuclear correlation term in the exact factorization enables us to deal with quantum coherences by accounting for electronic and nuclear nonadiabatic couplings effectively within classical nuclei approximation. We compare coupled- and independent-trajectory approximations with each other to understand algorithms in description of the bifurcation of nuclear wave packets and the correct spatial distribution of electronic wave functions along with nuclear trajectories. Finally, we show numerical results for comparisons of coupled- and independent-trajectory approaches for the photoisomerization of a protonated Schiff base from excited state molecular dynamics (ESMD) simulations with the recently developed Python-based ESMD code, namely, the PyUNIxMD program.

Atomic transistors based on seamless lateral metal-semiconductor junctions with a sub-1-nm transfer length.

The edge-to-edge connected metal-semiconductor junction (MSJ) for two-dimensional (2D) transistors has the potential to reduce the contact length while improving the performance of the devices. However, typical 2D materials are thermally and chemically unstable, which impedes the reproducible achievement of high-quality edge contacts. Here we present a scalable synthetic strategy to fabricate low-resistance edge contacts to atomic transistors using a thermally stable 2D metal, PtTe2. The use of PtTe2 as an epitaxial template enables the lateral growth of monolayer MoS2 to achieve a PtTe2-MoS2 MSJ with the thinnest possible, seamless atomic interface. The synthesized lateral heterojunction enables the reduced dimensions of Schottky barriers and enhanced carrier injection compared to counterparts composed of a vertical 3D metal contact. Furthermore, facile position-selected growth of PtTe2-MoS2 MSJ arrays using conventional lithography can facilitate the design of device layouts with high processability, while providing low contact resistivity and ultrashort transfer length on wafer scales.

Machine Learning-Assisted Excited State Molecular Dynamics with the State-Interaction State-Averaged Spin-Restricted Ensemble-Referenced Kohn-Sham Approach.

We present a machine learning-assisted excited state molecular dynamics (ML-ESMD) based on the ensemble density functional theory framework. Since we represent a diabatic Hamiltonian in terms of generalized valence bond ansatz within the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS) method, we can avoid singularities near conical intersections, which are crucial in excited state molecular dynamics simulations. We train the diabatic Hamiltonian elements and their analytical gradients with the SchNet architecture to construct machine learning models, while the phase freedom of off-diagonal elements of the Hamiltonian is cured by introducing the phase-less loss function. Our machine learning models show reasonable accuracy with mean absolute errors of ∼0.1 kcal/mol and ∼0.5 kcal/mol/Šfor the diabatic Hamiltonian elements and their gradients, respectively, for penta-2,4-dieniminium cation. Moreover, by exploiting the diabatic representation, our models can predict correct conical intersection structures and their topologies. In addition, our ML-ESMD simulations give almost identical result with a direct dynamics at the same level of theory.

Study of the Decoherence Correction Derived from the Exact Factorization Approach for Nonadiabatic Dynamics.

We present a detailed study of the decoherence correction to surface hopping that was recently derived from the exact factorization approach. Ab initio multiple spawning calculations that use the same initial conditions and the same electronic structure method are used as a reference for three molecules: ethylene, the methaniminium cation, and fulvene, for which nonadiabatic dynamics follows a photoexcitation. A comparison with the Granucci-Persico energy-based decoherence correction and the augmented fewest-switches surface-hopping scheme shows that the three decoherence-corrected methods operate on individual trajectories in a qualitatively different way, but the results averaged over trajectories are similar for these systems.

Surface Hopping Dynamics beyond Nonadiabatic Couplings for Quantum Coherence.

Description of correct electron-nuclear couplings is crucial in modeling of nonadiabatic dynamics. Within traditional semiclassical or mixed quantum-classical dynamics, the coupling between quantum electronic states and classical nuclear trajectories is governed by nonadiabatic coupling vectors coupled to classical nuclear momenta. This enables us to develop a very powerful nonadiabatic dynamics algorithm, namely, surface hopping dynamics, which can describe the splitting of nuclear wave packets and detailed balance. Despite its efficiency and practicality, it suffers from the lack of quantum decoherence due to incorrect accounts for the electron-nuclear coupling. Here we present a new surface hopping algorithm based on the exact electron-nuclear correlation from the exact factorization of molecular wave functions. This algorithm demands comparable computational costs to existing surface hopping methods. Numerical simulations with two-state models and a multidimensional multistate realistic molecule show that the electron-nuclear coupling beyond the nonadiabatic coupling terms can describe the quantum coherence properly.