Second-order Perturbation Method Research Articles

Advancing Our Understanding of the Excited States of the Tantalum Anion.

Our study provides a comprehensive theoretical examination of the energy levels associated with the neutral tantalum atom and its ions in various charge states (Ta, Ta+, and Ta-), employing the multiconfiguration Dirac-Hartree-Fock (MCDHF) method, and relativistic infinite order two-component (IOTC) method with multiconfiguration complete active space self-consistent field (CASSCF) followed by the second-order single-state multireference perturbation (CASPT2) methods. The effect of spin-orbit (SO) coupling is introduced via the restricted active space state interaction (RASSI) method, utilizing atomic mean field SO integrals (AMFI). Through IOTC CASSCF/CASPT2 RASSI calculations, we determined the electron affinity (EA) of the tantalum atom to be 0.321 eV, which stands among the most accurate theoretical values achieved to date. This result closely aligns with the experimental measurement of 0.329 eV. Our investigation highlights potential discrepancies between the predicted symmetry of the excited states of the tantalum anion and experimental observations. Additionally, we calculated the bonding energies for transitions from Ta- to Ta and identified four potential bound or quasi-bound states in the tantalum anion.

XMS-CASPT2//XMS-CASPT2 and XMS-CASPT2//CASSCF at comparison: The impact of dynamic correlation in the excited state optimization of nitronaphthalene.

Here, analytical extended multi-state complete active space second-order perturbation method (XMS-CASPT2) gradients are used to rationalize the decreasing triplet quantum yield trend in 2-nitronaphthalene, 1-nitronaphthalene, and 2-methyl-1-nitronaphthalene, a series of nitro-substituted aromatic compounds. Comparison with the XMS-CASPT2//CASSCF (where CASSCF stands for complete active space self-consistent field method) results highlights the importance of dynamic correlation in geometry optimization and challenges the validity of an XMS-CASPT2//CASSCF approach: XMS-CASPT2 S1 minima leads to planar structures, while CASSCF optimizations trigger a pyramidalization of the nitro group. The XMS-CASPT2 results correlate the reported decreasing triplet quantum yield trend in these species to a decrease in S1 to T2 population transfer and an increase in S1-S0 decay, while no such correlation is observed when using XMS-CASPT2//CASSCF data.

Robust design optimization using a non-intrusive second-order approximation of stochastic moments

This paper presents a new formulation of the second-order fourth-moment method (sometimes referred to as second-order perturbation method or second-order method of moments). The method allows to efficiently predict the stochastic moments of a response function and is therefore often used within robust design optimization. The new approach allows a non-intrusive implementation at the same cost as existing, highly intrusive formulations. Therefore, the new approach can be applied to any objective function without significant implementation effort. It is based on a few finite difference steps into special directions and hence is dependent on the corresponding step sizes. An automatic step size procedure is supplied beside a detailed convergence analysis. The advantages of the new formulation are demonstrated by robust design optimizations of a 2D and a 3D example using the geometrically nonlinear finite element method.

Mechanical analysis for deepwater drilling riser system with structural parameters uncertainty

Uncertain structural parameters increase the mechanical failure risk of the drilling riser system. The effect of uncertain structural parameters on the mechanical characteristics of the drilling riser system is unclear. Therefore, the mechanical characteristics of the drilling riser system with structural parameter uncertainty are analyzed. A mechanical model with uncertain structural parameters is established based on the basic analysis model and the local average theory. The basic analysis model takes into account the combined effects of platform motion, tensioner, riser string, and marine environment. The local average theory is suitable for situations where there is a lack of uncertainty information. A solution approach is proposed using the second-order central perturbation method and the Newmark-β method. The solution method reduces the load by simplifying degrees of freedom and random variables. The correctness of the solution method is verified using the Monte Carlo method. The analysis results indicate that uncertain structural parameters significantly affect the static response, while the dynamic response is less affected. The failure risk due to static axial stress and static shear stress is heightened by uncertain structural parameters. The failure risk arising from static shear stress at the bottom of the riser system experiences a significant increase.

New physical insights into the supporting subspace factorization of XMS-CASPT2 and generalization to multiple spin states via spin-free formulation.

This paper introduces a spin-free formulation of the supporting subspace factorization [C. Song and T. J. Martínez, J. Chem. Phys. 149, 044108 (2018)], enabling a reduction in the computational scaling of the extended multi-state complete active space second-order perturbation (XMS-CASPT2) method for arbitrary spins. Compared to the original formulation that is defined in the spin orbitals and is limited to singlet states, the spin-free formulation in this work treats different spin states equivalently, thus naturally generalizing the idea beyond singlet states. In addition, we will present a new way of deriving the supporting subspace factorization with the purpose of understanding its physical interpretation. In this new derivation, we separate the sources that make CASPT2 difficult into the "same-site interactions" and "inter-site interactions." We will first show how the Kronecker sum can be used to remove the same-site interactions in the absence of inter-site interactions, leading to MP2 energy in dressed orbitals. We will then show how the inter-site interactions can be exactly recovered using Löwdin partition, where the supporting subspace concept will naturally arise. The new spin-free formulation maintains the main advantage of the supporting subspace factorization, i.e., allowing XMS-CASPT2 energies to be computed using highly optimized MP2 energy codes and Fock build codes, thus reducing the scaling of XMS-CASPT2 to the same scaling as MP2. We will present and discuss results that benchmark the accuracy and performance of the new method. To demonstrate how the new method can be useful in studying real photochemical systems, the supporting subspace XMS-CASPT2 is applied to a photoreaction sensitive to magnetic field effects. The new spin-free formulation makes it possible to calculate the doublet and quartet states required in this particular photoreaction mechanism.

Interval strategy-based regularization approach for force reconstruction with multi-source uncertainties

A force reconstruction method for uncertain structures instead of direct measurement is proposed in this paper. Considering the multi-source uncertainties as unknown-but-bounded (UBB) parameters including unknown structural dynamics and measurement errors, this study proposed a novel force reconstruction method using an interval strategy-based regularization approach with truncated singular value decomposition (TSVD). To overcome the drawbacks of incompleted models and limited information, the overall modelling and inversion of the interval force are reconstructed based on the non-probabilistic uncertainty theory. First, The interval uncertainty-based kernel function for force inversion is deduced using the second-order interval perturbation method, in which the uncertainties of the reconstructed force originate from the uncertain singular values and measurement errors. Based on the regularization of the TSVD, the interval bounds of the singular values can be accurately obtained. Combining the two typical regularization approaches, a novel interval strategy-based regularization approach is proposed to optimize the best regularization parameter based on the uncertain interval distribution of singular values. This novel strategy can not only balance the two typical regularization approaches but also reduce the burden of force inversion. Two numerical examples including a truss and space capsule are applied to access the proposed method, and the estimated accurate force can prove the effectiveness.

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories.

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

Comparative study on free vibration analysis of rotating bi-directional functionally graded beams using multiple beam theories with uncertainty considerations

The present study investigates the free vibration behavior of rotating beams made of functionally graded materials (FGMs) with a tapered geometry. The material properties of the beams are characterized by an exponential distribution model. The stiffness and mass matrices of the beams are derived using the principle of virtual energy. These matrices are then evaluated using three different beam theories: Bernoulli–Euler (BE) or Classical Beam Theory (CBT), Timoshenko (T) or First-order Shear Deformation Theory (FSDT), and Reddy (R) or Third-order Shear Deformation Theory (TSDT). Additionally, the study incorporates uncertainties in the model parameters, including rotational velocity, beam material properties, and material distribution. The mean-centered second-order perturbation method is employed to account for the randomness of these properties. To ensure the robustness and accuracy of the probabilistic framework, numerical examples are presented, and the results are compared with those obtained through the Monte Carlo simulation technique. The investigation explores the impact of critical parameters, including material distribution, taper ratios, aspect ratio, hub radius, and rotational speed, on the natural frequencies of the beams is explored within the scope of this investigation. The outcomes are compared not only with previously published research findings but also with the results of 3-Dimensional Finite Element (3D-FE) simulations conducted using ANSYS to validate the model’s effectiveness. The comparisons demonstrate a strong agreement across all evaluations. Specifically, it is observed that for thick beams, the results obtained from FSDT and TSDT exhibit a greater agreement with the 3D-FE simulations compared to CBT. It is shown that the coefficient of variation (C.O.V.) of first mode eigenvalue of TSDT, FSDT and CBT are approximately identical for random rotational velocity and discernible deviations are noted in CBT compared to FSDT and TSDT in the case of random material properties. The findings suggest that TSDT outperforms FSDT by eliminating the need for a shear correction coefficient, thereby establishing its superiority in accurately predicting the natural frequencies of rotating, tapered beams composed of FGMs.

Theoretical Understanding on the Facilitated Photoisomerization of a Carbonyl Supported Borane System.

Boron compound BOMes2 containing an internal B-O bond undergoes highly efficient photoisomerization, followed by sequential structural transformations, resulting in a rare eight-membered B, O-heterocycle (S. Wang, et al. Org. Lett.2019, 21, 5285-5289). In this work, the detailed reaction mechanisms of such a unique carbonyl-supported tetracoordinate boron system in the first excited singlet (S1) state and the ground (S0) state were investigated by using the complete active space self-consistent field and its second-order perturbation (MS-CASPT2//CASSCF) method combined with time-dependent density functional theory (TD-DFT). Moreover, an imine-substituted tetracoordinated organic boron system (BNMes2) was selected for comparative study to explore the intrinsic reasons for the difference in reactivity between the two types of compounds.Steric factorwasfound to influence the photoisomerization activity of BNMes2 and BOMes2. These results rationalize the experimental observations and can provide helpful insights into understanding the excited-state dynamics of heteroatom-doped tetracoordinate organoboron compounds, which facilitates the rational design of boron-based materials with superior photoresponsive performances.

Reactions of 1,2-Azaborinine, a BN-Benzyne, with Organic π Systems.

Ortho-benzyne and 1,2-azaborinine are related by the formal exchange of the CC triple bond by the isoelectronic BN unit. The (2 + 2) and (2 + 4) cycloaddition reactions of 1,2-azaborinine with the different organic π systems (ethene, ethyne, 1,3-butadiene, 1,3-cyclopentadiene, furan, benzene) were examined computationally using density functional, second-order perturbation, and coupled-cluster methods. All reactions of 1,2-azaborinine with the studied substrates are highly exothermic and involve the formation of Lewis acid-base complexes of 1,2-azaborinine and respective π systems. The interaction between the π bond of the substrates and the empty p orbital of the boron atom in these complexes is remarkably strong, resulting in two-step mechanisms for the (2 + 2) and (2 + 4) cycloaddition reactions. Cycloaddition reactions have lower barriers than CH insertion reactions, and (2 + 4) reactions are favored over (2 + 2) cycloadditions.

Emergence of bond-dependent highly anisotropic magnetic interactions in Sr4RhO6 : A theoretical study

The quantum spin liquid states as a natural ground state of the Kitaev model has led to a quest for new materials candidates hosting Kitaev physics. Yet, there are very few material candidates in this category. Using a combination of $ab$ $initio$ and model Hamiltonian methods, we propose that Ruddlesden-Popper compound Sr$_4$RhO$_6$ belongs to this category. With a tight-binding model and exact diagonalization approach, we show that despite substantial trigonal-like distortion, the electronic and magnetic properties of Sr$_4$RhO$_6$ can be well described in terms of pseudo-spin = 1/2 states. Magnetic interactions among pseudo-spins, estimated using the second-order perturbation method are highly bond-dependent anisotropic in nature with two particularly noticeable features, antiferromagnetic Kitaev and Dzyaloshinskii-Moriya interactions. The gaped spin-wave spectra of Sr$_4$RhO$_6$ obtained with linear spin-wave theory is consistent with the underlying magnetic frustration. Additional analysis of the role of individual or a particular combination of magnetic interactions reveals that the spin-wave spectra of Sr$_4$RhO$_6$ is a combined effect of the highly anisotropic interactions and a relatively simpler minimal model may not be plausible in the current case. The crucial insights about coupling between the local structural features and magnetic properties of Sr$_4$RhO$_6$ obtained in this study may be helpful for future studies belonging to this class.

A theoretical study on the photochemical generation of phenylborylene from phenyldiazidoborane.

Organic borylenes are a kind of highly reactive species, which play important roles in a lot of reactions as vigorous intermediates. In this work, we investigated the photochemical generation mechanisms of phenylborylene (PhB) together with the side product N-phenylnitrenoiminoborane (PhNBN) from phenyldiazidoborane (PhBN6) by extrusion of dinitrogen in the two lowest electronic singlet states (S0 and S1) based on the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods combined with time-dependent density functional theory (TD-DFT) calculations. Our results show that the reaction PhBN6 → PhB + 3N2 involves stepwise N2 extrusion three times and the azido region rearrangement. Moreover, we found that the studied photo-induced processes are kinetically feasible because the highest energy barrier is only 0.36 eV and excitation with light of wavelength 254 nm can provide enough excess energy to overcome these energy barriers. Importantly, we revealed that several conical intersections between S1 and S0 states participate and facilitate the studied photochemical processes. Our results not only clarify the experimental observations, (H. F. Bettinger, J. Am. Chem. Soc. 2006, 128, 2534), but also provide valuable insights into borylene chemistry.

Understanding of Photo-Induced Reversible Rearrangement from Borepin to Borirane.

The first reversible photoisomerization between a borepin and a borirane was reported in the photo-induced reactions of B(npy)Ar2 (npy=2-(naphthalen-1-yl) pyridine, Ar=phenyl or electron rich aryl; S. Wang, et al. Angew. Chem. Int. Ed. 2019, 58, 6683-6687). In this work, the detailed mechanisms of the unprecedented reversible photoisomerization between the borepin (compound a) and the borirane (compound b) of B(npy)Ph2 in the first excited singlet (S1 ) state and the ground (S0 ) state were studied by carrying out calculations with the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods combined with time-dependent density functional theory (TD-DFT). The calculation results show that photoexcitation of a-S0 at 365 nm and b-S0 at 450 nm populate their S1 state with evident charge transfer characteristics. The photoisomerization is triggered in the S1 state and ends in the S0 state, at which the intersection points in a (S1 /S0 )x intersection seam participate in and promote phenyl migration and ring-closure processes. Furthermore, we reveal that the not large energy difference (less than 0.6 eV) and similar conjugation properties of π electrons between a-S0 and b-S0 are responsible for their unique photo-reversible reactivity, compared with those of the isomers of the thermally reversible compound B(ppy)Mes2 . Our results contribute to an understanding of the excited-state reactivity of organoboron compounds and will be useful to support the design of new boron-based photo-responsive materials.

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method.

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-cumulative distribution function method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with 92% strong scalability going from 64 to 4096 processors.

Excited-state double proton transfer of 1,8-dihydroxy-2-naphthaldehyde: A MS-CASPT2//CASSCF study

Excited-state double proton transfer (ESDPT) is a controversial issue which has long been plagued with theoretical and experimental communities. Herein, we took 1,8-dihydroxy-2-naphthaldehyde (DHNA) as a prototype and used combined complete active space self-consistent field (CASSCF) and multi-state complete active-space second-order perturbation (MS-CASPT2) methods to investigate ES-DPT and excited-state deactivation pathways of DHNA. Three different tautomer minima of S1-ENOL, S1-KETO-1, and S1-KETO-2 and two crucial conical intersections of S1S0-KETO-1 and S1S0-KETO-2 in.and between the S0 and S1 states were obtained. S1-KETO-1 and S1-KETO-2 should take responsibility for experimentally observing dual-emission bands. In addition, two-dimensional potential energy surfaces (2D-PESs) and linear interpolated internal coordinate paths connecting relevant structures were calculated at the MS-CASPT2//CASSCF level and confirmed a stepwise ESDPT mechanism. Specifically, the first proton transfer from S1-ENOL to S1-KETO-1 is barrierless, whereas the second one from S1-KETO-1 to S1-KETO-2 demands a barrier of ca. 6.0 kcal/mol. The linear interpolated internal coordinate path connecting S1-KETO-1 (S1-KETO-2) and S1S0-KETO-1 (S1S0-KETO-2) is uphill with a barrier of ca. 12.0 kcal/mol, which will trap DHNA in the S1 state while therefore enabling dual-emission bands. On the other hand, the S1/S0 conical intersections would also prompt the S1 system to decay to the S0 state, which could be to certain extent suppressed by locking the rotation of the C5−C8−C9−O10 dihedral angle. These mechanistic insights are not only helpful for understanding ESDPT but also useful for designing novel molecular materials with excellent photoluminescent performances.

Exact-Two-Component Relativistic Multireference Second-Order Perturbation Theory.

As the relativistic corrections become stronger for late-row elements, the fully perturbative treatment of spin-orbit coupling and dynamic correlation may become inadequate for accurate descriptions of chemical properties. In this work, we introduce a determinant-based Kramers-unrestricted exact-two-component multireference second-order perturbation (X2C-MRPT2) method which variationally includes relativistic corrections with a perturbative dynamic correlation. The restricted active space partitioning scheme is employed to provide an adjustable correlation space for the second-order perturbation treatment. The multistate perturbation theory is also developed to improve the descriptions of ground and excited states. Benchmark studies of atomic fine-structure splittings and spectroscopic constants of molecular monohydrides using X2C-MRPT2 are compared to the other perturbative and variational approaches. The results suggest that X2C-MRPT2 is a highly accurate alternative to the fully variational multireference configuration interaction method at only a small fraction of the computational cost.

Substituent and Solvent Effects on the Photoisomerization of Cinnamate Derivatives: An XMS-CASPT2 Study.

Cinnamate derivatives show a variety of photo-induced reactions. Among them is trans-cis photoisomerization, which may involve the nonradiative decay (NRD) process. The extended multistate complete active space second-order perturbation (XMS-CASPT2) method was used in this study as a suitable theory for treating multireference electronic nature, which was frequently manifested in the photoisomerization process. The minimum energy paths of the trans-cis photoisomerization process of cinnamate derivatives were determined, and the activation energies were estimated using the resulting intrinsic reaction coordinate (IRC) paths. Natural orbital analysis revealed that the transition state's (TS) electronic structure is zwitterionic-like, elucidating the solvent and substituent effect on the energy barrier of photoisomerization paths. Furthermore, it was found that the charge on the pyramidalized carbon atom at the TS structure was strongly correlated with the activation energy barrier for the cinnamate derivatives. Thus, it seemingly provided a physical picture of the photoisomerization of cinnamates and was a good descriptor potentially applicable to molecular design for controlling the rate constant of the photoisomerization reaction.

Electronic Structure of Nitrobenzene: A Benchmark Example of the Accuracy of the Multi-State CASPT2 Theory.

The electronic structure of nitrobenzene (C6H5NO2) has been studied by means of the complete active space self-consistent field (CASSCF) and multi-state second-order perturbation (MS-CASPT2) methods. To this end, an active space of 20 electrons distributed in 17 orbitals has been selected to construct the reference wave function. In this work, we have calculated the vertical excitation energies and the energy barrier for the dissociation of the molecule on the ground state into phenyl and nitrogen dioxide. After applying the corresponding vibrational corrections to the electronic energies, it is demonstrated that the MS-CASPT2//CASSCF values obtained in this work yield an excellent agreement between calculated and experimental data. In addition, other active spaces of lower size have been applied to the system in order to check the active space dependence in the results.

Improving Perturbation Theory for Open-Shell Molecules via Self-Consistency

We present an extension of our one-body Møller-Plesset second-order perturbation (OBMP2) method for open-shell systems. We derived the OBMP2 Hamiltonian through the canonical transformation followed by the cumulant approximation to reduce many-body operators into one-body ones. The resulting Hamiltonian consists of an uncorrelated Fock (unperturbed Hamiltonian) and a one-body correlation potential (perturbed Hamiltonian) composed of only double excitations. Molecular orbitals and associated energy levels are then relaxed via self-consistency, similar to Hartree-Fock, in the presence of the correlation at the MP2 level. We demonstrate the OBMP2 performance by considering two examples well-known for requiring orbital optimization: bond breaking and isotropic hyperfine coupling constants. In contrast to noniterative MP2, we show that OBMP2 can yield a smooth transition through the unrestriction point and accurately predict isotropic hyperfine coupling constants.

A new homotopy approach for stochastic static model updating with large uncertain measurement errors

Large measurement errors are a major challenge in structural model updating. Based on the concept of homotopy, a new stochastic static model updating method is proposed to update structural models using uncertain static data with large measurement errors. First, considering the uncertainty of the static measurement, a stochastic model updating equation for element update factors is set up. To solve the stochastic model updating equation, a series of homotopy deformation equations are presented to establish the relationship between the deterministic update factors and the random update factors. Furthermore, the homotopy series expansions of the random update factors can be determined by solving the homotopy deformation equations. Since the measured degrees of freedom of updated structures are usually limited or unavailable, a static condensation technique is used for stochastic model updating. To address the ill-posed problems caused by incomplete measurement information and static measurement errors, the Tikhonov regularization method is used in the process of solving the homotopy deformation equations. Three numerical examples are given to demonstrate the validity of the proposed stochastic model updating method. The numerical results clearly show that unlike the second-order perturbation method, this new method can produce better accuracy in cases with large measurement errors. Compared with the Bayesian method with the Delayed Rejection Adaptive Metropolis sampling technique and even a fast Bayesian method, the proposed method utilizes much less computational time and provides an equivalent accuracy. Finally, static loading experiments on a two-span continuous concrete beam are implemented to validate the proposed model updating method.