Second-order Perturbation Theory Research Articles

Emergence of Exceptional Points in Periodic Metastructures With Hidden Parity-Time Symmetric Defects

Abstract We study the elastodynamics of a periodic metastructure incorporating a defect pair that enforces a parity-time (PT) symmetry due to judiciously engineered imaginary impedance elements—one having energy amplification (gain) and the other having an equivalent attenuation (loss) mechanism. We show that their presence affects the initial band structure of the periodic Hermitian metastructure and leads to the formation of numerous exceptional points (EPs) which are mainly located at the band edges where the local density of modes is higher. The spatial location of the PT-symmetric defect serves as an additional control over the number of emerging EPs in the corresponding spectra as well as the critical non-Hermitian (gain/loss) strength required to create the first EP—a specific defect location minimizes the critical non-Hermitian strength. We use both finite element and coupled-mode-theory-based models to investigate these metastructures and use a time-independent second-order perturbation theory to further demonstrate the influence of the size of the metastructure and the PT-symmetric defect location on the minimum non-Hermitian strength required to create the first EP in a band. Our findings motivate feasible designs for the experimental realization of EPs in elastodynamic metastructures.

Microsolvation of electrons by a handful of ammonia molecules.

Microsolvation of electrons in ammonia is studied here via anionic NH3 n - clusters with n = 2-6. Intensive samplings of the corresponding configurational spaces using second-order perturbation theory with extended basis sets uncover rich and complex energy landscapes, heavily populated by many local minima in tight energy windows as calculated from highly correlated coupled cluster methods. There is a marked energetical preference for structures that place the excess electron external to the molecular frame, effectively coordinating it with the three protons from a single ammonia molecule. Overall, as the clusters grow in size, the lowest energy dimer serves as the basic motif over which additional ammonia molecules are attached via unusually strong charge-assisted hydrogen bonds. This is apriori quite unexpected because, on electrostatic grounds, the excess electron would be expected to be in contact with as many protons as possible. Accordingly, a full quantum mechanical treatment of the bonding interactions under the tools provided by the quantum theory of atoms in molecules is carried out in order to dissect and understand the nature of intermolecular contacts. Vertical detachment energies reveal bound electrons even for n = 2.

E and B modes of the CMB y-type distortions: polarised kinetic Sunyaev-Zeldovich effect from the reionisation and post-reionisation eras

We study the E and B mode polarisation of the cosmic microwave background (CMB) originating from the transverse peculiar velocity of free electrons, at second order in perturbation theory, during the reionisation and post-reionisation eras. Interestingly, the spectrum of this polarised kinetic Sunyaev-Zeldovich (SZ) effect can be decomposed into a blackbody part and a y-type distortion. The y-distortion part is distinguishable from the primary E and B modes and also the lensing B modes. Furthermore, it is also differentiable from the other y-type signals, such as the thermal SZ effect, which are unpolarised. We show that this signal is sensitive to the reionisation history, in particular to how fast reionisation happens. The E and B modes of y-type distortion provide a way to beat the cosmic variance of primary CMB anisotropies and are an independent probe of the cosmological parameters. The blackbody component of the pkSZ effect would be an important foreground for the primordial tensor modes for tensor to scalar ratio r ≲ 3 × 10-5.

Biharmonic (ω, 2ω) ionization of atoms by elliptically-polarized light. Carving the photoelectron angular distributions

The biharmonic (ω, 2ω) photoionization of atomic inner-shell electrons opens up new perspectives for studying nonlinear light–atom interactions at intensities in the transition regime from weak to strong-field physics. In particular, the control of the frequency and polarization of biharmonic beams enables one to carve the photoelectron angular distribution and to enhance the resolution of ionization measurements by the (simultaneous) absorption of photons. Apart from its quite obvious polarization dependence, the photoelectron angular distributions are sensitive also to the (relative) intensity, the phase difference and the temporal structure of the incoming beam components, both at resonant and nonresonant frequencies. Here, we describe and analyze several characteristic features of biharmonic ionization in the framework of second-order perturbation theory and (so-called) ionization pathways, as they are readily derived from the interaction of inner-shell electrons with the electric-dipole field of the incident beam. We show how the photoelectron angular distribution and elliptical dichroism can be shaped in rather an unprecedented way by just tuning the properties of the biharmonic field. Since such fields are nowadays accessible from high-harmonic sources or free-electron lasers, these and further investigations might help extract photoionization amplitudes or the phase difference of incoming beams.

Theoretical Spectroscopic Study of Two Ketones of Atmospheric Interest: Methyl Glyoxal (CH3COCHO) and Methyl Vinyl Ketone (CH3COCH═CH2).

Two ketones of atmospheric interest, methyl glyoxal and methyl vinyl ketone, are studied using explicitly correlated coupled cluster theory and core–valence correlation-consistent basis sets. The work focuses on the far-infrared region. At the employed level of theory, the rotational constants can be determined to within a few megahertz of the experimental data. Both molecules present two conformers, trans/cis and antiperiplanar (Ap)/synperiplanar (Sp), respectively. trans-Methyl glyoxal and Ap-methyl vinyl ketone are the preferred structures. cis-Methyl glyoxal is a secondary minimum of very low stability, which justifies the unavailability of experimental data in this form. In methyl vinyl ketone, the two conformers are almost isoenergetic, but the interconversion implies a relatively high torsional barrier of 1798 cm–1. A very low methyl torsional barrier was estimated for trans-methyl glyoxal (V3 = 273.6 cm–1). Barriers of 429.6 and 380.7 cm–1 were computed for Ap- and Sp-methyl vinyl ketone. Vibrational second-order perturbation theory was applied to determine the rovibrational parameters. The far-infrared region was explored using a variational procedure of reduced dimensionality. For trans-methyl glyoxal, the ground vibrational state was estimated to split by 0.067 cm–1, and the two low excited energy levels (1 0) and (0 1) were found to lie at 89.588 cm–1/88.683 cm–1 (A2/E) and 124.636 cm–1/123.785 cm–1 (A2/E). For Ap- and Sp-methyl vinyl ketone, the ground vibrational state splittings were estimated to be 0.008 and 0.017 cm–1, respectively.

Electronic states and transitions of PrO and PrO+ probed by threshold ionization spectroscopy and spin-orbit multiconfiguration perturbation theory.

The precise ionization energy of praseodymium oxide (PrO) seeded in supersonic molecular beams is measured with mass-analyzed threshold ionization (MATI) spectroscopy. A total of 33 spin-orbit (SO) states of PrO and 23 SO states of PrO+ are predicted by second-order multiconfigurational quasi-degenerate perturbation (MCQDPT2) theory. Electronic transitions from four low-energy SO levels of the neutral molecule to the ground state of the singly charged cation are identified by combining the MATI spectroscopic measurements with the MCQDPT2 calculations. The precise ionization energy is used to reassess the ionization energies and the reaction enthalpies of the Pr+ O → PrO+ + e- chemi-ionization reaction reported in the literature. An empirical formula that uses atomic electronic parameters is proposed to predict the ionization energies of lanthanide monoxides, and the empirical calculations match well with available precise experimental measurements.

Covariant transverse-traceless projection for secondary gravitational waves

Second-order tensor modes induced by nonlinear gravity are a key component of the cosmological background of gravitational waves. A detection of this background would allow us to probe the primordial power spectrum at otherwise inaccessible scales. Usually, the energy density of these gravitational waves is studied within perturbation theory in a particular gauge -- a connection between our physical spacetime and a fictitious background. It is a widely recognized issue that the second-order, scalar-induced gravitational waves are gauge dependent. This issue arises because they are not well-defined as tensors in the physical spacetime at second-order and are thus unphysical. In this paper, we propose the covariant transverse-traceless projection of the extrinsic curvature to study cosmological gravitational waves on a spatial hypersurface. We define a new energy density which is based purely on spacetime tensors, independent of perturbation theory, and thus is gauge invariant by definition. We show that, in the context of second-order perturbation theory, this new energy density contains only propagating modes in the constant-time hypersurface in the Newtonian gauge. We further show that we can recover the same gravitational waves after a transformation to the synchronous gauge, so long as we correctly identify the Newtonian hypersurface.

Large-Scale Benchmarking of Multireference Vertical-Excitation Calculations via Automated Active-Space Selection.

We have calculated state-averaged complete-active-space self-consistent-field (SA-CASSCF), multiconfiguration pair-density functional theory (MC-PDFT), hybrid MC-PDFT (HMC-PDFT), and n-electron valence state second-order perturbation theory (NEVPT2) excitation energies with the approximate pair coefficient (APC) automated active-space selection scheme for the QUESTDB benchmark database of 542 vertical excitation energies. We eliminated poor active spaces (20–40% of calculations) by applying a threshold to the SA-CASSCF absolute error. With the remaining calculations, we find that NEVPT2 performance is significantly impacted by the size of the basis set the wave functions are converged in, regardless of the quality of their description, which is a problem absent in MC-PDFT. Additionally, we find that HMC-PDFT is a significant improvement over MC-PDFT with the translated PBE (tPBE) density functional and that it performs about as well as NEVPT2 and second-order coupled cluster on a set of 373 excitations in the QUESTDB database. We optimized the percentage of SA-CASSCF energy to include in HMC-PDFT when using the tPBE on-top functional, and we find the 25% value used in tPBE0 to be optimal. This work is by far the largest benchmarking of MC-PDFT and HMC-PDFT to date, and the data produced in this work are useful as a validation of HMC-PDFT and of the APC active-space selection scheme. We have made all the wave functions produced in this work (orbitals and CI vectors) available to the public and encourage the community to utilize this data as a tool in the development of further multireference model chemistries.

Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database.

Double-hybrid density functional theory (DHDFT) offers a pathway to accuracy approaching composite wavefunction approaches such as G4 theory. However, the Görling–Levy second-order perturbation theory (GLPT2) term causes them to partially inherit the slow ∝L–3 (with L the maximum angular momentum) basis set convergence of correlated wavefunction methods. This could potentially be remedied by introducing F12 explicit correlation: we investigate the basis set convergence of both DHDFT and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large and chemically diverse general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) benchmark suite. The B2GP-PLYP-D3(BJ) and revDSD-PBEP86-D4 DHDFs are investigated as test cases, together with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets as large as cc-pVQZ-F12. We show that F12 greatly accelerates basis set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12 basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPPD in orbital-based approaches, and in fact comparable in quality to cc-pV(5+d)Z. Somewhat surprisingly, aug-cc-pVDZ-F12 is not required even for the anionic subsets. In conclusion, DHDF-F12/VDZ-F12 eliminates concerns about basis set convergence in both the development and applications of double-hybrid functionals. Mass storage and I/O bottlenecks for larger systems can be circumvented by localized pair natural orbital approximations, which also exhibit much gentler system size scaling.

Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit.

Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange–correlation functionals. If the exact weak-interaction expansion from the second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work, we fill this gap by presenting a fully self-consistent-field (SCF) implementation of some popular ACM functionals. While using second-order perturbation theory at weak interactions, we have also introduced new generalized gradient approximations (GGAs), beyond the usual point-charge-plus-continuum model, for the first two leading terms at strong interactions, which are crucial to ensure robustness and reliability. We then assess the SCF–ACM functionals for molecular systems and for prototypical strong-correlation problems. We find that they perform well for both the total energy and the electronic density and that the impact of SCF orbitals is directly connected to the accuracy of the ACM functional form. For the H2 dissociation, the SCF–ACM functionals yield significant improvements with respect to standard functionals also thanks to the use of the new GGAs for the strong-coupling functionals.

Electronic absorption spectra from off-diagonal quantum master equations.

Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrödinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrödinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from the electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion method.

Accurate property prediction by second order perturbation theory: The REMP and OO-REMP hybrids.

The prediction of molecular properties such as equilibrium structures or vibrational wavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing Retaining the Excitation Degree - Møller-Plesset (REMP) and Orbital Optimized REMP (OO-REMP) hybrid perturbation theories, showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers and dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations. OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 and 0.5pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20cm-1. Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3%. Across all investigated properties, it turns out that a 20%:80% Møller-Plesset:Retaining the Excitation Degree mixing ratio consistently provides the best results. This is in line with our previous findings, featuring closed and open shell reaction energies.

On the origin of bonding in metals: lithium as a case study

The bonding in metals is analysed within the framework of the PATMOS (Perturbed AToms in MOlecules and Solids) model. The electronic binding energy per atom is written as a sum of a distortion energy of the atoms and the partitioned interaction energy comprising Coulombic, exchange and correlation terms. The adopted physical model of the infinite system, is spherical embedding of the atoms of the reference unit cell. Correlation energies are calculated by second-order Møller-Plesset and second-order Epstein-Nesbet perturbation theory. The binding energy of lithium solid is calculated for 16 nearest neighbour distances from 4.0 to 10.0 Bohr. Electron correlation is of paramount importance for the binding energy. The calculated cohesive energy is ( 0.0571 ± 0.004 ) Hartree comparing with experimental value 0.0599 Hartree. The bonding picture is characterised by slightly expanded atomic orbitals.

Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths

Currently, there is a growing interest in the development of a new hierarchy of methods based on the concept of seniority, which has been introduced quite recently in quantum chemistry. Despite the enormous potential of these methods, the accurate description of both dynamical and static correlation effects within a single and in-principle-exact approach remains a challenge. In this work, we propose an alternative formulation of reduced density-matrix functional theory (RDMFT) where the (one-electron reduced) density matrix is mapped onto an ab initio seniority-zero wave function. In this theory, the exact natural orbitals and their occupancies are determined self-consistently from an effective seniority-zero calculation. The latter involves a universal higher-seniority density matrix functional for which an adiabatic connection (AC) formula is derived and implemented under specific constraints that are related to the density matrix. The pronounced curvature of the (constrained) AC integrand, which is numerically observed in prototypical hydrogen chains and the Helium dimer, indicates that a description of higher-seniority correlations within second-order perturbation theory is inadequate in this context. Applying multiple linear interpolations along the AC or connecting second-order perturbation theory to a full-seniority treatment via Pad\'{e} approximants are better strategies. Such information is expected to serve as a guide in the future design of higher-seniority density-matrix functional approximations.

Multiconfiguration Pair-Density Functional Theory for Chromium(IV) Molecular Qubits.

Pseudotetrahedral organometallic complexes containing chromium(IV) and aryl ligands have been experimentally identified as promising molecular qubit candidates. Here we present a computational protocol based on multiconfiguration pair-density functional theory for computing singlet–triplet gaps and zero-field splitting (ZFS) parameters in Cr(IV) aryl complexes. We find that two multireference methods, multistate complete active space second-order perturbation theory (MS-CASPT2) and hybrid multistate pair-density functional theory (HMS-PDFT), perform better than Kohn–Sham density functional theory for singlet–triplet gaps. Despite the very small values of the ZFS parameters, both multireference methods performed qualitatively well. MS-CASPT2 and HMS-PDFT performed particularly well for predicting the trend in the ratio of the rhombic and axial ZFS parameters, |E/D|. We have also investigated the dependence and sensitivity of the calculated ZFS parameters on the active space and the molecular geometry. The methodologies outlined here can guide future prediction of ZFS parameters in molecular qubit candidates.

Multireference Study of Optically Addressable Vanadium-Based Molecular Qubit Candidates.

Molecular electron spin qubits with optical manipulation schemes are some of the most promising candidates for modern quantum technologies. Key values that determine a compound's viability for optical-spin initialization and readout include its singlet-triplet gap and zero-field splitting (ZFS) parameters. Generally, these values are very small in magnitude and are thus difficult to reproduce with theoretical methods. Here, we study a previously identified optically addressable molecular qubit, (C6F5)3trenVCNtBu (tren = tris(2-aminoethyl)amine), using the complete active space self-consistent field (CASSCF) and post-CASSCF methods: complete active space second-order perturbation theory (CASPT2), multiconfiguration pair-density functional theory (MC-PDFT), and hybrid MC-PDFT (HMC-PDFT). Of those methods, we successfully reproduce the singlet-triplet gap and ZFS parameters with reasonable accuracy using 0.5 HMC-PDFT and CASPT2. Four additional V3+ complexes with differing ligands were also investigated. We found that the ligands have minimal effect on the spin properties of the molecule and propose them to be optically addressable qubit candidates. These potential qubits are further analyzed in terms of ab initio ligand field theory (AILFT) to understand the influence of the ligands on the singlet-triplet gap and ZFS parameters.

Numerical implementation of the Cubic Galileon model in pinocchio

ABSTRACT We present a perturbative treatment of non-linear galaxy clustering in the context of the cubic Galileon modified gravity model, in terms of second-order Lagrangian Perturbation theory and an extension of ellipsoidal collapse that includes Vainshtein screening. We numerically implement such prescriptions in the approximate pinocchio code, and use it to generate realizations of the matter density field and halo catalogues with different prescriptions for ellipsoidal collapse. We investigate the impact of three different approximations in the computation of collapse times on the halo mass function, halo bias, and matter power spectrum. In the halo mass function, both the modified gravity effect and the screening effect are significant in the high-mass end, similar to what is found for other MG models. We perform a comparison with N-body simulations to assess the validity of our approach, and show that we can reproduce the same trend observed in simulations for all quantities considered. With a simple modification to the grouping algorithm of pinocchio to take into account the gravity modification, and without the need to re-calibrate the algorithm, we show that we can reproduce the linear halo bias and the mildly non-linear matter power spectrum of simulations with good accuracy, especially for the implementation with Vainshtein screening. We stress that, while approximate, our method is orders of magnitude faster than a full N-body simulation, making it an optimal tool for the quick generation of large sets of halo catalogues for cosmological observables.

Resolving Ultrafast Photoinitiated Dynamics of the Hachimoji 5-aza-7-deazaguanine Nucleobase: Impact of Synthetically Expanding the Genetic Alphabet†.

The guanine derivative, 5-aza-7-deazaguanine (5N7C G) has recently been proposed as one of four unnatural bases, termed Hachimoji (8-letter) to expand the genetic code. We apply steady-state and time-resolved spectroscopy to investigate its electronic relaxation mechanism and probe the effect of atom substitution on the relaxation mechanism in polar protic and polar aprotic solvents. Mapping of the excited state potential energy surfaces is performed, from which the critical points are optimized by using the state-of-art extended multi-state complete active space second-order perturbation theory. It is demonstrated that excitation to the lowest energy 1 ππ* state of 5N7C G results in complex dynamics leading to ca. 10- to 30-fold slower relaxation (depending on solvent) compared with guanine. A significant conformational change occurs at the S1 minimum, resulting in a 10-fold greater fluorescence quantum yield compared with guanine. The fluorescence quantum yield and S1 decay lifetime increase going from water to acetonitrile to propanol. The solvent-dependent results are supported by the quantum chemical calculations showing an increase in the energy barrier between the S1 minimum and the S1 /S0 conical intersection going from water to propanol. The longer lifetimes might make 5N7C G more photochemically active to adjacent nucleobases than guanine or other nucleobases within DNA.

A local pair-natural orbital-based complete-active space perturbation theory using orthogonal localized virtual molecular orbitals.

The multireference second-order perturbation theory (CASPT2) is known to deliver a quantitative description of various complex electronic states. Despite its near-size-consistent nature, the applicability of the CASPT2 method to large, real-life systems is mostly hindered by large computational and storage costs for the two-external tensors, such as two-electron integrals, amplitudes, and residuum. To this end, Menezes and co-workers developed a reduced-scaling CASPT2 scheme by incorporating the local pair-natural orbital (PNO) representation of the many-body wave functions using non-orthonormal projected atomic orbitals (PAOs) into the CASPT theory [F. Menezes et al., J. Chem. Phys. 145, 124115 (2016)]. Alternatively, in this paper, we develop a new PNO-based CASPT2 scheme using the orthonormal localized virtual molecular orbitals (LVMOs) and assess its performance and accuracy in comparison with the conventional PAO-based counterpart. Albeit the compactness, the LVMOs were considered to perform somewhat poorly compared to PAOs in the local correlation framework because they caused enormously large orbital domains. In this work, we show that the size of LVMO domains can be rendered comparable to or even smaller than that of PAOs by the use of the differential overlap integrals for domain construction. Optimality of the MOs from the CASSCF treatment is a key to reducing the LVMO domain size for the multireference case. Due to the augmented Hessian-based localization algorithm, an additional computational cost for obtaining the LVMOs is relatively minor. We demonstrate that the LVMO-based PNO-CASPT2 method is routinely applicable to large, real-life molecules such as Menshutkin SN2 reaction in a single-walled carbon nanotube reaction field.

Phase control of localization in the nonlinear two-mode system from harmonic mixing driving: Perturbative analysis and symmetry consideration

In this paper, we present a rigorous analysis of symmetry and underlying physics of the nonlinear two-mode system driven by a harmonic mixing field, by means of multiple scale asymptotic analysis method. The effective description in the framework of the second-order perturbative theory provides an accurate picture for understanding the Floquet eigenspectrum and dynamical features of the nonlinear two-mode system, showing full agreement with the prediction of symmetry considerations. We find that two types of symmetries play significant role in the dynamical features of this model, the mechanism behind which can be interpreted in terms of the effective description. The results are of relevance for the phase control of the atomic localization in Bose–Einstein condensates or switch of the optical signals in nonlinear mediums.