New diabatic potential energy surfaces of the NaH2 system and dynamics studies for the Na(3p) + H2 → NaH + H reaction

A diabatization method based upon integrating the diabatic potential gradient difference

In this work we develop a new scheme to construct a diabatic potential energy matrix (DPEM). We propose a diabatization method which is based on integrating the diabatic potential gradient difference to diabatize adiabatic ab initio energies. This method is capable of performing high-precision adiabatic-to-diabatic transformations, with a unique advantage in effectively handling the significant fluctuations in derivative-couplings caused by conical intersection (CI) seams. The above scheme is applied to the DPEM construction of the Na(3p) + H2 → NaH + H reaction. The fitting data including adiabatic energies, energy gradients and derivative-couplings obtained from a previous benchmark DPEM are diabatized and fitted using a general neural network fitting procedure to generate the DPEM. The produced DPEM can effectively describe nonadiabatic processes involving different electronic states. We further perform quantum dynamical calculations on the new DPEM and the previous benchmark DPEM, and the obtained results demonstrate the effectiveness of our scheme.

Constructing Diabatic Potential Energy Matrices with Neural Networks Based on Adiabatic Energies and Physical Considerations: Toward Quantum Dynamic Accuracy

A permutation invariant polynomial-neural network (PIP-NN) approach for constructing the global diabatic potential energy matrices (PEMs) of the coupled states of molecules is proposed. Specifically, the diabatization scheme is based merely on the adiabatic energy data of the system, which is ideally a most convenient way due to not requiring additional ab initio calculations for the data of the derivative coupling or any other physical properties of the molecule. Considering the permutation and coupling characteristics of the system, particularly in the presence of conical intersections, some vital treatments for the off-diagonal terms in diabatic PEM are essentially needed. Taking the photodissociation of H₂O(X~/B~)/NH₃(X~/A~) and nonadiabatic reaction Na(3p) + H₂ → NaH(Σ⁺) + H for example, this PIP-NN method is shown to build up the global diabatic PEMs effectively and accurately. The root-mean-square errors of the adiabatic potential energies in the fitting for three different systems are all small (<10 meV). Further quantum dynamic calculations show that the absorption spectra and product branching ratios in both H₂O(X~/B~) and NH₃(X~/A~) nonadiabatic photodissociation are well reproduced on the new diabatic PEMs, and the nonadiabatic reaction probability of Na(3p) + H₂ → NaH(Σ⁺) + H obtained on the new diabatic PEMs of the 1²A₁ and 1²B₂ states is in reasonably good agreement with previous theoretical result as well, validating this new PIP-NN method.

Unusual presentation of lepidic adenocarcinoma in a healthy female

Background

Lepidic adenocarcinoma represents a histologic pattern of non-small cell lung cancer that characteristically arises in the lung periphery with tracking alongside pre-existing alveolar walls. Noninvasive and invasive variants of lepidic adenocarcinoma are dependent on parenchymal destruction, vascular, or pleural invasion. The lepidic-predominant lung malignancies are collectively recognized as slow growing with rare metastasis and excellent prognosis. The World Health Organization classification of lung malignancies depends on molecular and histopathological findings. CT findings most commonly include ground-glass characteristics, commonly mistaken for inflammatory or infectious etiology. These tumors are generally surgically resectable and associated with better survival given infrequent nodal and extrathoracic involvement. Rarely these tumors present with diffuse pneumonic-type involvement associated with worse outcomes despite lack of nodal and distant metastases.

Case presentation

A 63-year-old Caucasian athletic immunocompetent female presented with 2 months of progressive shortness of breath, fatigue, loss of appetite and 15 pound weight loss. History was only notable for well controlled essential hypertension and hypothyroidism. Contrast computed tomography angiogram and positron emission tomography revealed diffuse hypermetabolic interstitial and airspace abnormalities of the lungs without lymphadenopathy (or distant involvement) in addition to right hydropneumothorax and left pleural effusion. Baseline laboratory testing was unremarkable, and extensive bacterial and fungal testing returned negative. Bronchoscopy and video-assisted thoracoscopic surgery was subsequently performed with pleural fluid cytology, lung and pleural biopsies returning positive for lepidic adenocarcinoma with 2% programmed death ligand 1 expression and genomic testing positive for PTEN gene deletion. Prior to treatment, the patient perished on day 15 of admission.

Conclusion

We present a rare case of lepidic predominant adenocarcinoma with extensive bilateral aerogenous spread in the context of no lymphovascular invasion in a healthy, low risk patient. This case presentation may add to the body of knowledge regarding the different behavior patterns of lepidic predominant adenocarcinomas.

Wave Packet Approach to Adiabatic and Nonadiabatic Dynamics of Cold Inelastic Scatterings

Due to the extremely large de Broglie wavelength of cold molecules, cold inelastic scattering is always characterized by the time-independent close-coupling (TICC) method. However, the TICC method is difficult to apply to collisions of large molecular systems. Here, we present a new strategy for characterizing cold inelastic scattering using wave packet (WP) method. In order to deal with the long de Broglie wavelength of cold molecules, the total wave function is divided into interaction, asymptotic and long-range regions (IALR). The three regions use different numbers of ro-vibrational basis functions, especially the long-range region, which uses only one function corresponding to the initial ro-vibrational state. Thus, a very large grid range can be used to characterize long de Broglie wavelengths in scattering coordinates. Due to its better numerical scaling law, the IALR-WP method has great potential in studying the inelastic scatterings of larger collision systems at cold and ultracold regimes.

Stereodynamics-Controlled Product Branching in the Nonadiabatic H + NaD → Na(3s, 3p) + HD Reaction at Low Temperatures

Nonadiabatic processes play an important role at energies near or higher than conical intersection of adiabatic potential energy surfaces in chemical reactions. In this work, dynamics of the nonadiabatic H + NaD reaction at low temperatures are studied by using the quantum wave packet method based on an improved L-shaped grid. The nonadiabatic H + NaD reaction has two exothermic reaction channels: Na(3s) + HD and Na(3p) + HD; the latter can only occur via nonadiabatic transition. The dynamics results show that the product branching of the H + NaD reaction at collision energies ranging from 20 to 80 cm^-1 is controlled by stereodynamics. The Na(3s) and Na(3p) reaction channels occur through collinear collision and side-on collision, respectively. When the collision energy is lower than 20 cm^-1, the resonance-mediated reaction mechanism is dominant in both the Na(3s) and Na(3p) reaction channels.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Quantum Dynamics Studies of the Significant Intramolecular Isotope Effects on the Nonadiabatic Be+(2P) + HD → BeH+/BeD+ + D/H Reaction

Quantum time-dependent wave packet dynamics studies on the nonadiabatic Be⁺(²P) + HD → BeH⁺/BeD⁺ + D/H reaction are performed for the first time employing recently constructed diabatic potential energy surfaces. Strong intramolecular isotope effects and unusual results are presented, which are attributed to the dynamic effects of shallow wells induced by avoided crossing on the diagonal V₂₂^d surface. The BeH⁺ + D and BeD⁺ + H channels are dominated by high-J and low-J partial waves, respectively. The BeD⁺/BeH⁺ branching ratio is larger than 10 at low energy and gradually decreases with increasing collision energy. The BeH⁺ product is primarily distributed at low vibrational states, whereas there exists an obvious population inversion of vibrational states on the BeD⁺ product. The results of differential cross sections suggest that the formation of the BeH⁺ + D channel favors a direct reaction process, while the BeD⁺ + H channel is mainly generated by the complex-forming mechanism.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Quantum dynamics studies of isotope effects in the Mg+(3p) + HD → MgH+/MgD+ + D/H insertion reaction

The time-dependent wave packet quantum dynamics studies for the Mg+(3p) + HD → MgH+/MgD+ + D/H diabatic reaction are carried out for the first time on recently developed diabatic YHWCH potential energy surfaces [Phys. Chem. Chem. Phys., 2018, 20, 6638-6647]. The results of reaction probabilities and total integral cross sections show a dramatic preference to the formation of MgD+ over MgH+ owing to the insertion reaction mechanism in the title reaction. The MgD+/MgH+ branching ratio witnesses a monotonic decrease from 10.58 to 3.88 at collision energy range of 0.01 to 0.20  eV, and at the collision energy of 0.114  eV, it is close to the experimental value of 5. The rovibrational state-resolved ICSs of the two channels show the products MgD+ have higher vibrational excitation and hotter rotational state distributions. The opacity function P(J) suggests that the MgH+ + D channel and MgD+ + H channel are dominated by high-b and low-b collisions, respectively. Both forward and backward scattering peaks are found in the differential cross section curves, whereas the angle distributions of products are not strictly forward-backward symmetric because of the short lifetime of the complex in the reaction.

Dynamics study on the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction: insertion-abstraction mechanism

Time-dependent wave packet calculations are carried out for two reaction channels of the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction. The potential well on the excited state potential energy surface makes the reaction preferable to proceed through the insertion reaction path. The dominance of the NaD + H reaction channel and product rotational state distributions are found to be in agreement with the characteristics of typical adiabatic insertion reactions. However, significant forward scattering peaks in the differential cross sections (DCS) are found to be inconsistent with the forward-backward symmetric scattering characteristic of typical adiabatic insertion reactions, which indicate that the Na(3p) + HD reaction is dominated by a direct reaction mechanism. The comparison between adiabatic and non-adiabatic calculated DCSs reveals that the non-adiabatic couplings in the reaction could reduce the lifetime of the intermediate complex. Finally, the insertion-abstraction mechanism is put forward for the non-adiabatic Na(3p) + HD reaction.