-
Penn Bredahl posted an update 6 months, 1 week ago
For the TADF molecules considered, the DF-QDPT2 method provides a mean absolute error (MAE) of 0.13 eV, while the MAE values of DF-SA-CASSCF and SS-SR-CASPT2 are 0.65 and 0.74 eV, respectively. The performances of B3LYP and PBE0 are slightly better than that of DF-QDPT2, while M06-2X and ωB97X-D provide noticeably higher errors compared with DF-QDPT2. Furthermore, the standard CASSCF without state-averaging yields dramatic errors with an MAE value of 3.0 eV. Our results demonstrate that eigenvalues of the DF-QDPT2-effective Hamiltonian can be reliably used for the prediction of singlet-triplet transition energies, while eigenvalues of DF-CASSCF/DF-SA-CASSCF fail to provide accurate predictions. Overall, we conclude that the DF-QDPT2 method emerges as a very useful tool for the computation of excited-state properties.The ion-molecule reaction H3+ + CO → H2 + HCO+/HOC+, which initiates the formation of crucial organic molecules, plays a key role in interstellar and circumstellar environments. In this work, the quasi-classical trajectory method is employed to study the reaction dynamics on a recently developed full-dimensional global potential energy surface (PES). The calculated product internal energy distributions and relative internal excited fractions agree reasonably well with the experimental measurements. For the two reaction channels, most of the available energy flows into the vibrational modes of HCO+ or HOC+ at low collision energies, followed by the translational mode and the rotational modes of HCO+ or HOC+. As the collision energy increases, the proportion of the product translational energy increases while the proportion of the product vibrational energy decreases. Furthermore, the CH and CO stretching modes and their combination bands are effectively excited for the product HCO+ while the bending mode is remarkably excited for the product HOC+.Using molecular dynamics (MD) simulations, we study the mechanism of stress corrosion cracking in graphene. Two sets of modelings are conducted. In the first one, large graphene sheets with cracks in the armchair and zigzag directions are exposed to oxygen molecules. The crack growth as a result of chemical reactions between carbon radicals and oxygen molecules at different mechanical tensile stress levels is studied. In the second set of simulations, MD simulations are combined with the density functional-based tight bonding method to enhance the accuracy. This set of modelings focuses on a smaller zone in the vicinity of the crack tip. The impact of initial crack orientation on corrosion is studied by investigating corrosion of cracks in both armchair and zigzag directions. We investigate the subcritical crack propagation occurring as a result of the combined effects of both mechanical loading and chemical reactions. Our results show that cracks in graphene can grow due to chemical reactions with the environmental molecules. The MD modelings also predict that reaction of carbon atoms with oxygen molecules might lead to a stress relaxation at the crack tip, hence preventing further crack propagation. The results show that subcritical crack growth can happen by two mechanisms, which include the failure of C-C bonds or by removing the carbon atoms from graphene sheets in the form of CO or CO2 molecules.We introduce the deep post Hartree-Fock (DeePHF) method, a machine learning-based scheme for constructing accurate and transferable models for the ground-state energy of electronic structure problems. DeePHF predicts the energy difference between results of highly accurate models such as the coupled cluster method and low accuracy models such as the Hartree-Fock (HF) method, using the ground-state electronic orbitals as the input. It preserves all the symmetries of the original high accuracy model. The added computational cost is less than that of the reference HF or DFT and scales linearly with respect to system size. We examine the performance of DeePHF on organic molecular systems using publicly available data sets and obtain the state-of-art performance, particularly on large data sets.The photophysical relaxation pathways of tzA, tzG, and tzI luminescent nucleobases were investigated with the MS-CASPT2 quantum-chemical method and double-ζ basis sets (cc-pVDZ) in gas and condensed phases (1,4-dioxane and water) with the sequential Monte Carlo/CASPT2 and free energy gradient (FEG) methods. Solvation shell structures, in the ground and excited states, were examined with the pairwise radial distribution function (G(r)) and solute-solvent hydrogen-bond networks. Site-specific hydrogen bonding analysis evidenced relevant changes between both electronic states. The three luminescent nucleobases share a common photophysical pattern, summarized as the lowest-lying 1(ππ*) bright state that is populated directly after the absorption of radiation and evolves barrierless to the minimum energy structure, from where the excess of energy is released by fluorescence. AUZ454 From the 1(ππ*)min region, the conical intersection with the ground state ((ππ*/GS)CI) is not accessible due to the presence of high energetic barriers. By combining the present results with those reported earlier by us for the pyrimidine fluorescent nucleobases, we present a comprehensive description of the photophysical properties of this important class of new fluorescent nucleosides.A 15-dimensional analytical form for the potential energy and dipole moment surfaces of the SF6 molecule in the ground electronic state is obtained using ab initio methods. In order to calculate the equilibrium S-F distance, we applied the coupled cluster CCSD(T) method and several versions of the correlation-consistent basis sets from valence triple-zeta (VTZ) and augmented valence triple-zeta (AVTZ) to core-valence quadruple-zeta (CVQZ) with Douglas-Kroll (DK) relativistic corrections that provided good agreement with an empirical equilibrium value. Ab initio electronic energies on 15D grids of nuclear geometries are computed using the CCSD(T) method with VTZ and CVQZ-DK basis sets. The analytical representation of the potential energy surface is determined through an expansion in symmetry-adapted products of nonlinear coordinates up to the 5th order. The influence of additional redundant coordinates on the quality of the fit was investigated. Parameters of full-dimensional dipole moment surfaces are determined up to the 4th order expansion in normal mode coordinates.