Search for: All records

The role of the starting configurations in nonadiabatic dynamics is investigated, and it is found for the first time that the dissociation times of the surface hopping trajectories initializing with the reactant are much greater than those starting from the transition state; however, the two sets of chemiexcitation yields are almost the same within standard errors in luminol chemiluminescence. The significant singlet chemiexcitation yield is explained by the difference between the number of S0 to S1 excitations and S1 to S0 de-excitations, and the negligible triplet can be explained similarly. The substantial singlet chemiexcitation yield obtained via the trajectory surface hopping method is attributed to the existence of the chemiexcitation region, through which the intrinsic reaction path does not pass. 

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. The methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT. 

The energy gaps, spin-orbit coupling (SOC) and admixture coefficients over a series of the configurations are evaluated by the SA-CASSCF/6-31G, SA-CASSCF/6-31G*, SA-CASSCF/ANO-RCC-VDZP and MS-CASPT2/ANO-RCC-VDZP to reveal the extent of the inaccuracy of the SA-CASSCF. By comparing the mean absolute errors for the energy gaps and the admixture coefficient magnitudes (ACMs) measured between the SA-CASSCF/6-31G, SA-CASSCF/6-31G* or SA-CASSCF/ANO-RCC-VDZP and the MS-CASPT2/ANO-RCC-VDZP, the SA-CASSCF/6-31G is selected as the electronic structure method in the nonadiabatic molecular dynamics simulation. The major components of the ACMs of the SA-CASSCF/6-31G and MS-CASPT2/ANO-RCC-VDZP are identified and compared, we find that the ACMs are underestimated by the SA-CASSCF/6-31G, which is verified by the reasonable triplet quantum yield simulated by the trajectory surface hopping and the calibrated SA-CASSCF/6-31G. The magnitude of the singlet-triplet mixing positively correlates to the hopping probability between the mixed singlet and triplet states, which is confirmed by the computed S-T transition probability. 

We propose the generalized initial velocity sampling algorithm at a transition state, in which the total initial kinetic energy and extra positive initial velocity along the reaction coordinate have been introduced to improve the accuracy and efficiency of nonadiabatic dynamics simulation. This sampling algorithm is very useful for chemical reactions with multiple transition states, e.g., two transition states in the thermolysis of four-membered heterocyclic peroxides. The dependence of the chemiexcitation yields, dissociation times, and other quantities on the total initial kinetic energy and extra positive initial velocity has been investigated. By taking different CASPT2 corrections and additional positive initial velocities into account, we found that the resulting triplet quantum yield matches the experimental result perfectly. In most ensembles, the secondary primary intersystem crossing, i.e., the first singlet excited state to the first triplet state, is a major direct production channel for the first triplet state product trajectories. 

Creating analytic representations of multiple potential energy surfaces for modeling electronically nonadiabatic processes is a major challenge being addressed in various ways by the chemical dynamics community. In this work, we introduce a new method that can achieve convenient learning of multiple potential energy surfaces (PESs) and their gradients (negatives of the forces) for a polyatomic system. This new method, called compatibilization by deep neural network (CDNN), is demonstrated to be accurate and, even more importantly, to be automatic. The only required input is a database with geometries and potential energies. The method produces a matrix, called the compatible potential energy matrix (CPEM), that may be interpreted as the electronic Hamiltonian in an implicit nonadiabatic basis, and the analytic adiabatic potential energy surfaces and their gradients are obtained by diagonalization and automatic differentiation. We show that the CPEM, which is neither adiabatic nor necessarily diabatic, can be discovered automatically during the learning procedure by the special design of a CDNN architecture. We believe that the CDNN method will be very useful in practice for learning coupled PESs for polyatomic systems because it is accurate and fully automatic. 

The thermolysis of trans-3,4-dimethyl-1,2-dioxetane is studied by trajectory surface hopping. The significant difference between long and short dissociation times is rationalized by frustrated dissociations and the time spent in triplet states. If the C−C bond breaks through an excited state channel, then the trajectory passes over a ridge of the potential energy surface of that state. The calculated triplet quantum yields match the experimental results. The dissociation half-times and quantum yields follow the same ascending order as per the product states, justifying the conjecture that the longer dissociation time leads to a higher quantum yield, proposed in the context of the methylation effect. The populations of the molecular Coulomb Hamiltonian and diagonal states reach equilibrium, but the triplet populations with different Sz components fluctuate indefinitely. Certain initial velocities, leading the trajectories to given product states, can be identified as the most characteristic features for sorting trajectories according to their product states. 

Deep learning (DL) models have demonstrated state-of-the-art performance in the classification of diagnostic imaging in oncology. However, DL models for medical images can be compromised by adversarial images, where pixel values of input images are manipulated to deceive the DL model. To address this limitation, our study investigates the detectability of adversarial images in oncology using multiple detection schemes. Experiments were conducted on thoracic computed tomography (CT) scans, mammography, and brain magnetic resonance imaging (MRI). For each dataset we trained a convolutional neural network to classify the presence or absence of malignancy. We trained five DL and machine learning (ML)-based detection models and tested their performance in detecting adversarial images. Adversarial images generated using projected gradient descent (PGD) with a perturbation size of 0.004 were detected by the ResNet detection model with an accuracy of 100% for CT, 100% for mammogram, and 90.0% for MRI. Overall, adversarial images were detected with high accuracy in settings where adversarial perturbation was above set thresholds. Adversarial detection should be considered alongside adversarial training as a defense technique to protect DL models for cancer imaging classification from the threat of adversarial images.