Calculations Of Potential Energy Surfaces Research Articles

Effects of triaxial deformation on the fission barrier in the Z = 118 − 120 nuclei* *Supported by the National Natural Science Foundation of China (Grant Nos. 12205076 and 12047504), the China Postdoctoral Science Foundation (Grant No. 2020M670012), and the Launching Fund of Henan University of Technology (Grant No. 2021BS047).

By using potential energy surface (PES) calculations in the three-dimensional space (β 2, γ, β 4) within the framework of the macroscopic-microscopic model, the fission trajectory and fission barrier for Z = 118(Og), 119, 120 nuclei has been systematically investigated. The calculated PES includes macroscopic liquid-drop energy, microscopic shell correction and pairing correction. Taking the 294Og176 nucleus as an example, we discuss the next closed shell after Z = 82 and N = 126 with the calculated Woods–Saxon single-particle levels. Then, the results of PES in 294Og is illustrated from the (X, Y) scale to the (β 2, γ) scale. The γ degree of freedom reveals the shape evolution clearly during the fission process. The structure near the minimum and saddle point of the PES in the Z = 118, 119, 120 nuclei is demonstrated simultaneously. Based on the potential energy curves, general trends of the evolution of the fission barrier heights and widths are also studied. The triaxial deformation in these superheavy mass regions plays a vital role in the first fission barrier, showing a significant reduction in both triaxial paths. In addition, the model-dependent fission barriers of proton-rich nuclei 295Og, 296119, and 297120 are analyzed briefly. Our studies could be valuable for synthesizing the superheavy new elements in the forthcoming HIAF and other facilities.

From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra.

Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.

Inelastic scattering of NO(A2Σ+) + CO2: rotation-rotation pair-correlated differential cross sections.

A crossed beam velocity-map ion-imaging apparatus has been used to determine differential cross sections (DCSs) for the rotationally inelastic scattering of NO(A2Σ+, v = 0, j = 0.5) with CO2, as a function of both NO(A, v = 0, N') final state and the coincident final rotational energy of the CO2. The DCSs are dominated by forward-peaked scattering for all N', with significant rotational excitation of CO2, and a small backward scattered peak is also observed for all final N'. However, no rotational rainbow scattering is observed and there is no evidence for significant product rotational angular momentum polarization. New ab initio potential energy surface calculations at the PNO-CCSD(T)-F12b level of theory report strong attractive forces at long ranges with significant anisotropy relative to both NO and CO2. The absence of rotational rainbow scattering is consistent with removal of low-impact-parameter collisions via electronic quenching, in agreement with the literature quenching rates of NO(A) by CO2 and recent electronic structure calculations. We propose that high-impact-parameter collisions, that do not lead to quenching, experience strong anisotropic attractive forces that lead to significant rotational excitation in both NO and CO2, depolarizing product angular momentum while leading to forward and backward glory scattering.

Theoretical study on the insertion reaction of CH(X2) into the oh bond in n-C4H9OH

CH radicals play an important role in the combustion of hydrocarbon. The insertion mechanism of a CH radical into the O-H bond of n-C4H9OH is investigated theoretically by a detailed potential energy surface calculation at the BHandHLYP/6-311++G(3df,2p) and CCSD(T)/6-311++G(d,p) (single-point) levels. Our results show that the CH radical attacks into the oxygen atom in n-C4H9OH to form a prereaction complex (COMP) to be followed by an insertion of the CH radical into the O-H bond of the n-C4H9OH molecule to form the low-lying intermediate IS1 (CH2OCH2CH2CH2CH3). This intermediate can isomerize to form IS2 (CH3OCH2CHCH2CH3), IS3 (CH3CH2CH2CH2CH2O), and IS4 (CH3CH2CH2CH2CHOH). These intermediates can decompose to yield 9 products (PR1-PR9) in which major ones are PR1 (CH2CH2CH2CH3 + CH2O), PR2 (CH2CHCH2CH3 + CH3O) and PR3 (CH2CHCH2OCH3 + CH3).

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2).

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Development of a scalar-based geometric parameterization approach for the crystal structure landscape of dithienylethene-based crystalline solids.

Dithienylethenes (DTEs) are a promising class of organic photoswitches that can be used to create crystalline solids with properties controlled by light. However, the ability of DTEs to adopt multiple conformations, only one of which is photoactive, complicates the rational design of these materials. Herein, the synthesis and structural characterization of 19 crystalline solids containing a single DTE molecule are described. A novel D-D analysis of the molecular geometries obtained from rotational potential energy surface calculations and the ensemble of experimental structures were used to construct a crystal landscape for DTE. Of the 19 crystal structures, 17 contained photoinactive DTE rotamers and only 2 were photoactive. These results highlight the challenges associated with the design of these materials. Overall, the D-D analysis described herein provides rapid, effective and intuitive means of linking the molecular structure to photoactivity that could be applied more broadly to afford a general strategy for producing photoactive diarylethene-based crystalline solids.

A computational characterization of the reaction mechanisms for the reactions N(2D) + CH3CN and HC3N and implications for the nitrogen-rich organic chemistry of Titan

The reactions of atomic nitrogen in its first electronically excited state, N(2D), with acetonitrile (CH3CN) and cyanoacetylene (HC3N) have been investigated by performing electronic structure calculations of the doublet potential energy surfaces and RRKM estimates of the product branching fractions. According to our results, the insertion of N(2D) into one of the sigma CH bonds of acetonitrile leading to the formation of cyanomethanimine (also known as iminoacetonitrile) and the N(2D) addition to the CC triple bond of cyanoacetylene leading to the formation of dicyanocarbene are the main reaction pathways under the conditions typical of the upper atmosphere of Titan, the massive moon of Saturn. Other molecular products originating from other reaction pathways only give a minor contribution.Implications for the atmospheric chemistry of Titan, as well as implications in prebiotic chemistry and in the chemistry of the interstellar medium, are also noted.

Excited-state dynamics of deuterated indigo

Indigo, a rich blue dye, is an incredibly photostable molecule that has survived in ancient art for centuries. It is also unique in that it can undergo both an excited-state hydrogen and proton transfer on the picosecond timescale followed by a ground-state back transfer. Previously, we performed gas phase excited-state lifetime studies on indigo to study these processes in a solvent-free environment, combined with excited-state calculations. We found two decay pathways, a fast sub-nanosecond decay and a slow decay on the order of 10 ns. Calculations of the excited-state potential energy surface found that both hydrogen and proton transfer are nearly isoenergetic separated by a 0.1 eV barrier. To further elucidate these dynamics, we now report a study with deuterated indigo, using resonance-enhanced multi-photon ionization and pump-probe spectroscopy with mass spectrometric isotopomer selection. From new calculations of the excited-state potential energy surface, we find sequential double-proton or hydrogen transfer, whereby the trajectory to the second transfer passes a second barrier and then encounters a conical intersection that leads back to the ground state. We find that deuteration only increases the excited-state lifetimes of the fast decay channel, suggesting tunneling through the first barrier, while the slower channel is not affected and may involve a different intermediate state.Graphical abstract

Cluster aspects of collinear cluster tri-partition (CCT)

New results demonstrating the so-called Ni-bump being the most populated mode of the collinear cluster tripartition (CCT) of 252Cf(sf) are presented. The physical scenario of this effect is discussed. It is tested by calculations of potential energy surfaces for the fission of the intermediate fragments formed after first rupture of the mother nucleus. Fission barriers are extracted and mass asymmetries at saddles are compared with the masses of the fission fragments that take part in the Ni-bump.

Density Functional Theory Study of Synergistic Gas Sensing Using an Electrically Conductive Mixed Ligand Two-Dimensional Metal-Organic Framework.

Two-dimensional conductive metal-organic frameworks (2D-cMOFs) have been adopted in electrochemical sensing applications owing to their superior electrical conductivity and large surface area. Here, we performed a density functional theory (DFT) analysis to study the synergistic impact of introducing a secondary organic ligand to the 2D-cMOF system. In this study, cobalt-hexaiminobenzene (Co-HIB) and cobalt-2,3,6,7,10,11-hexaiminotriphenylene (Co-HITP) were combined to form a mixed ligand MOF named, Co-HIB-HITP. A DFT-level comparative study was designed to access stability, synergistic gas adsorption capability, and gas adsorption mechanism, important factors in sensing material development. A potential energy surface calculation predicted the structural stability of Co-HIB-HITP at larger interlayer displacements around 3.6-4.2 Å regions along the ab-plane than its unmixed states, Co-HIB and Co-HITP, indicating the tunability of the stacking mode using the mixed ligand system. Furthermore, the adsorption capabilities toward toxic gases, NH3, H2S, NO, and NO2, were investigated, and Co-HIB-HITP revealed superiority over unmixed 2D-cMOFs in H2S and NH3 gas adsorption energies by showing 158 and 170% improvement, respectively. Finally, an electron charge density analysis revealed Co-HIB-HITP's unique stacking mode and Co-metal density as contributing factors to its gas-selective synergy effect. The AB stacked layers and an intermediate metal density (5.25%) significantly improved the electrostatic interactions with H2S and NH3 by inducing a change in the chemical environment of the gas binding sites. This work proposes the dual-ligand 2D-cMOF as the promising design strategy for the next-generation sensing material.

Accelerating variational quantum eigensolver convergence using parameter transfer

One impediment to the useful application of variational quantum algorithms in quantum chemistry is slow convergence with large numbers of classical optimization parameters. In this work, we evaluate a quantum computational warm-start approach for potential energy surface calculations. Our approach, which is inspired by conventional computational methods, is evaluated using simulations of the variational quantum eigensolver. Significant speedup is demonstrated relative to calculations that rely on a Hartree–Fock initial state, both for ideal and sampled simulations. The general approach of transferring parameters between similar problems is promising for accelerating current and near-term quantum chemistry calculations on quantum hardware, and is likely applicable beyond the tested algorithm and use case.

Gaussian process regression adaptive density-guided approach: Toward calculations of potential energy surfaces for larger molecules.

We present a new program implementation of the Gaussian process regression adaptive density-guided approach [Schmitz et al., J. Chem. Phys. 153, 064105 (2020)] for automatic and cost-efficient potential energy surface construction in the MidasCpp program. A number of technical and methodological improvements made allowed us to extend this approach toward calculations of larger molecular systems than those previously accessible and maintain the very high accuracy of constructed potential energy surfaces. On the methodological side, improvements were made by using a Δ-learning approach, predicting the difference against a fully harmonic potential, and employing a computationally more efficient hyperparameter optimization procedure. We demonstrate the performance of this method on a test set of molecules of growing size and show that up to 80% of single point calculations could be avoided, introducing a root mean square deviation in fundamental excitations of about 3 cm-1. A much higher accuracy with errors below 1 cm-1 could be achieved with tighter convergence thresholds still reducing the number of single point computations by up to 68%. We further support our findings with a detailed analysis of wall times measured while employing different electronic structure methods. Our results demonstrate that GPR-ADGA is an effective tool, which could be applied for cost-efficient calculations of potential energy surfaces suitable for highly accurate vibrational spectra simulations.

RETRACTED:Exploring the intermolecular interactions in carbon disulfide dimer: An ab initio study using an improved Lennard–Jones potential energy surface for physical insights

AbstractCarbon disulfide dimer (CS2)2 is a model system that has been widely studied in the field of computational chemistry because of its relevance to a variety of chemical and biological processes. The (CS2)2 dimer is a relatively simple molecular system composed of two carbon disulfide (CS2) molecules interacting with each other through intermolecular forces. Despite its apparent simplicity, the (CS2)2 dimer exhibits a rich array of structural and dynamical properties that are of great interest to researchers. In this research, we present an ab initio study of the intermolecular interactions in the carbon disulfide dimer (CS2)2 using an improved Lennard–Jones (ILJ) potential with CCSD(T)/QZVPP calculations. The potential energy surface of (CS2)2 is calculated using high‐level quantum mechanical calculations based on the CCSD(T)/def2‐qzvpp method, which accurately accounts for electron correlation effects. The resulting potential energy surface is then fitted to an ILJ potential energy function, which includes both long‐range dipole–dipole interactions and short‐range repulsive interactions. The calculated potential energy surface reveals a rich variety of structural and dynamical properties of (CS2)2, including multiple minima and saddle points, which are sensitive to the relative orientation of the two CS2 molecules. It is essential to use extended basis sets to accurately incorporate the significant quadrupole moment of CS2, which we have calculated to be 2.44 a.u. The results of this study demonstrate the importance of using high‐level ab initio methods for the accurate calculation of potential energy surfaces in complex molecular systems such as (CS2)2. The use of an ILJ potential, which takes into account both dipole–dipole interactions and short‐range repulsive interactions, provides a more accurate and efficient approach for modeling intermolecular interactions in (CS2)2 and other similar systems. The results of this study will be useful for understanding the behavior of carbon disulfide dimers in different environments and for the development of new materials and chemical processes.

Experimental and Kinetic Study on the Gas-Phase Pyrolysis of Syringol

Syringyl (3,5-dimethoxy-4-hydroxyphenyl) units are abundant in hardwood lignin, which explains the large amount of syringol present in bio-oils derived from these feeds. To improve understanding of the decomposition chemistry of these syringyl units, the gas-phase pyrolysis kinetics of syringol have been studied experimentally and theoretically under typical fast pyrolysis conditions. The experiments were performed at 673–923 K in a micropyrolyzer unit hyphenated with a gas chromatography (GC) × gas chromatography-flame ionization detection/time-of-flight mass spectrometry and a customized gas chromatograph, which allows to detect and quantify products with boiling points up to 823 K next to permanent gases and water. New obtained experimental data show that 2,3-dihydroxybenzaldehyde is the most important phenolic product, and a large amount of pyrolytic water was formed. The proposed decomposition pathways of syringol to these two products show two new features: (1) 2,3-dihydroxybenzaldehyde is formed from a unimolecular decomposition of syringol without any stable intermediates and (2) both pathways leading to 2,3-dihydroxybenzaldehyde and water involve one more H-atom migration step to form the more stable radicals. Potential energy surface calculations at the CBS-QB3 level of theory indicate that the new features are kinetically favorable. Furthermore, a detailed kinetic model for syringol pyrolysis including 361 reactions and 93 species was constructed. Thermodynamic data for important species and kinetics for crucial reactions were determined with ab initio calculations at the CBS-QB3 level of theory. The model predicts well the conversion of syringol and the formation of major products. A rate of production analysis shows that the new proposed pathway is the major decomposition pathway of syringol. These findings shed new light on the pyrolysis and oxidation chemistry of syringyl species.

Cooling dynamics of energized naphthalene and azulene radical cations.

Naphthalene and azulene are isomeric polycyclic aromatic hydrocarbons (PAHs) and are topical in the context of astrochemistry due to the recent discovery of substituted naphthalenes in the Taurus Molecular Cloud-1 (TMC-1). Here, the thermal- and photo-induced isomerization, dissociation, and radiative cooling dynamics of energized (vibrationally hot) naphthalene (Np+) and azulene (Az+) radical cations, occurring over the microsecond to seconds timescale, are investigated using a cryogenic electrostatic ion storage ring, affording "molecular cloud in a box" conditions. Measurement of the cooling dynamics and kinetic energy release distributions for neutrals formed through dissociation, until several seconds after hot ion formation, are consistent with the establishment of a rapid (sub-microsecond) Np+ ⇌ Az+ quasi-equilibrium. Consequently, dissociation by C2H2-elimination proceeds predominantly through common Az+ decomposition pathways. Simulation of the isomerization, dissociation, recurrent fluorescence, and infrared cooling dynamics using a coupled master equation combined with high-level potential energy surface calculations [CCSD(T)/cc-pVTZ], reproduce the trends in the measurements. The data show that radiative cooling via recurrent fluorescence, predominately through the Np+ D0 ← D2 transition, efficiently quenches dissociation for vibrational energies up to ≈1eV above dissociation thresholds. Our measurements support the suggestion that small cations, such as naphthalene, may be more abundant in space than previously thought. The strategy presented in this work could be extended to fingerprint the cooling dynamics of other PAH ions for which isomerization is predicted to precede dissociation.

Infrared Spectra of the Water-CO2 Complex in the 4.3-3.6 μm Region and Determination of the Ground State Tunneling Splitting for HDO-CO2.

Spectra of water─CO2 dimers are studied using a tunable mid-infrared source to probe a pulsed slit jet supersonic expansion. H2O-CO2 and D2O-CO2 are observed in the CO2 ν3 fundamental region (≈2350 cm-1), D2O-CO2 is also observed in the D2O ν3 fundamental region (≈2790 cm-1), and HDO-CO2 is observed in the HDO O-D stretch fundamental region (≈2720 cm-1), all for the first time in these regions. Analysis of the spectra yields excited state rotational parameters and vibrational shifts. They also yield the first experimental values of the ground state internal rotation tunneling splittings for D2O-CO2 (0.003 cm-1) and HDO-CO2 (0.0234 cm-1). The latter value is a direct determination made possible by the reduced symmetry of HDO-CO2. These results provide stringent and easily interpreted tests for theoretical water-CO2 potential energy surface calculations.

Photodissociation of water molecule at short photon wavelengths: Dynamical studies

In our last study [J. Phys. B At. Mol. Opt. Phys. 54, 125,102 (2021).], we reported the ab initio calculation of the full-dimensional potential energy surfaces of water molecule including 9 A’ and 9 A” states in Cs symmetry. In this study, we performed additional non-adiabatic semi-classical studies based on the potential energy surfaces. Our simulation successfully repeated the near picosecond lifetime of the F∼1A′ state measured by time resolved photo-electron spectra experiment [Chinese J. Chem. Phys. 32, 53 (2019)]. We also determined the dissociation branching ratio including H + OH(X, A), H + H + O and H2+O channels. In addition, the reaction path corresponding to H2+O (1S) channel is clearly marked out, which is found in recent free-electron laser experiment [Nat. Commun. 12, 6,303 (2021)].

Microscopic study of neutron-induced fission process of 239Pu via zero- and finite-temperature density functional theory* *Supported by the National Natural Science Foundation of China (11790325, 12275081, 11790320, 11790321, 11961131010, 11605054, 12105369, 12147219, 12047568) and the Continuous Basic Scientific Research Project (WDJC-2019-09)

To study the neutron-induced fission of Pu, potential energy surface (PES) calculations were performed using zero and finite-temperature density functional theory (FT-DFT) with the Skyrme force. The energy of the incident neutron was simulated by the temperature of the FT-DFT. The variations of the least-energy fission path, fission barrier, total kinetic energy, scission line, and mass distribution of fission fragments with the incident neutron energy were analyzed. It was learned that an increase in the temperature lowers the barrier height, the isomeric-state energy, and the ridge between symmetric and asymmetric fission valleys. Additionally, the gaps of the single particle levels become smaller with an increase in the temperature. As the temperature increases, the pre-fission region shrinks, and the scission occurs at smaller deformation around the symmetric fission channel. At low temperatures, the pairing correlations in the collective space are similar to those in zero-temperature DFT, and when the temperature is 0.3 MeV, the pairing gaps decrease rapidly. Two different methods were used to calculate the fission yields of the neutron-induced fission Pu (n, f) with different incident neutron energies, in the framework of time-dependent generator coordinate method (TDGCM). One way to calculate the fission yield of Pu (n, f) is to solve the collective equation of the TDGCM by using the PES from the FT-DFT with the corresponding temperature. The other involves using the PES from the zero-temperature DFT and adjusting the initial collective energy of the wave packet in the TDGCM according to the incident neutron energy. For the cases of the lower incident neutron energies, these two methods gave similar results and reproduced the experimental peak and width of fission fragment distribution. However, for the highest incident neutron energy considered in this study, the results from the TDGCM using the PES from zero-temperature DFT deviated explicitly from the experimental data, whereas those obtained by using the PES from FT-DFT remained close to the experimental data. This indicated that, with the increase in the incident neutron energy, the shell structure of the compound nuclei changed explicitly; thus, it may not be effective to use the PES from zero-temperature to perform the fission dynamic calculation.

An experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction with implications for the photochemical models of Titan.

We report on a combined experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction, which is of relevance in the aromatic chemistry of the atmosphere of Titan. Experimentally, the reaction was studied (i) under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (Ec) of 31.8 kJ mol-1 to determine the primary products, their branching fractions (BFs), and the reaction micromechanism, and (ii) in a continuous supersonic flow reactor to determine the rate constant as a function of temperature from 50 K to 296 K. Theoretically, electronic structure calculations of the doublet C6H6N potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction mechanism. The reaction is found to proceed via barrierless addition of N(2D) to the aromatic ring of C6H6, followed by formation of several cyclic (five-, six-, and seven-membered ring) and linear isomeric C6H6N intermediates that can undergo unimolecular decomposition to bimolecular products. Statistical estimates of product BFs on the theoretical PES were carried out under the conditions of the CMB experiments and at the temperatures relevant for Titan's atmosphere. In all conditions the ring-contraction channel leading to C5H5 (cyclopentadienyl) + HCN is dominant, while minor contributions come from the channels leading to o-C6H5N (o-N-cycloheptatriene radical) + H, C4H4N (pyrrolyl) + C2H2 (acetylene), C5H5CN (cyano-cyclopentadiene) + H, and p-C6H5N + H. Rate constants (which are close to the gas kinetic limit at all temperatures, with the recommended value of 2.19 ± 0.30 × 10-10 cm3 s-1 over the 50-296 K range) and BFs have been used in a photochemical model of Titan's atmosphere to simulate the effect of the title reaction on the species abundances as a function of the altitude.

Molecular potential energy computation via graph edge aggregate attention neural network

Accurate potential energy surface (PES) calculation is the basis of molecular dynamics research. Using deep learning (DL) methods can improve the speed of PES calculation while achieving competitive accuracy to ab initio methods. However, the performance of DL model is extremely sensitive to the distribution of training data. Without sufficient training data, the DL model suffers from overfitting issues that lead to catastrophic performance degradation on unseen samples. To solve this problem, based on the message passing paradigm of graph neural networks, we innovatively propose an edge-aggregate-attention mechanism, which specifies the weight based on node and edge information. Experiments on MDI7 and QM9 datasets show that our model not only achieves higher PES calculation accuracy but also has better generalization ability compared with Schnet, which demonstrates that edge-aggregate-attention can better capture the inherent features of equilibrium and non-equilibrium molecular conformations.