Electric Field Gradient Calculations for Ice VIII and IX Using Polarizable Embedding: A Comparative Study on Classical Computers and Quantum Simulators

We test the performance of the polarizable embedding variational quantum eigensolver self-consistent field (PE-VQE-SCF) model for computing electric field gradients with comparisons to conventional complete active space self-consistent-field (CASSCF) calculations and experimental results. We compute quadrupole coupling constants for ice VIII and ice IX. We find close agreement of the quantum-computing PE-VQE-SCF results with the results from the classical PE-CASSCF calculations and with experiment. Furthermore, we observe that the inclusion of the environment is crucial for obtaining results that match the experimental data. The calculations for ice VIII are within the experimental uncertainty for both CASSCF and VQE-SCF for oxygen and lie close to the experimental value for ice IX as well. With the VQE-SCF, which is based on an adaptive derivative-assembled problem-tailored (ADAPT) ansatz, we find that the inclusion of the environment and the size of the different basis sets do not directly affect the gate counts. However, by including an explicit environment, the wavefunction and therefore the optimization problem become more complicated, which usually results in the need to include more operators from the operator pool, thereby increasing the depth of the circuit.

Electric Field Gradient in Chiral and Tetrahedral Molecules within High-Order LRESC Formalism

In this work, we present the electric field gradient (EFG) given by the linear response elimination of the small component (LRESC) scheme up to the 1/c4 order (c is the speed of light in vacuum) in CHFClX (X = Br, I, At) chiral molecules, together with CHF2Br and CH2FX (X = Br, I, At) tetrahedral systems. The former could be good candidates for further parity violation studies, especially when heavy atoms are surrounding. In this context, the LRESC scheme demonstrates effective applicability to large tetrahedral and chiral molecules that incorporate heavy elements, with relativistic effects playing a crucial role. The LRESC results of EFG exhibit an excellent agreement with those calculated at the four-component level, giving differences of only hundredths order in a.u. (atomic units) for the bromine nucleus and less than 0.1 a.u. for the iodine nucleus. Regarding the other nuclei, for the chiral molecules, there is a heavy atom effect on the light atom (HALA) for chlorine and fluorine atoms as the substituent halogen atom becomes heavier. Furthermore, the electronic part of the EFG for the central carbon and the fluorine nuclei presents an important dependence with the environment in the molecules under study. With accurate calculations of the EFG and tabulated nuclear quadrupole moment, the nuclear quadrupole coupling constant is obtained within the LRESC scheme, including for the first time correlation effects on the spin-dependent corrections with this methodology, providing results close to the experimental ones for Cl, Br, and I atoms. At the Hartree-Fock level, the differences are around 6% for Br and I nuclei, and at the density functional theory level with the LDA and PBE0 functionals, the differences are no more than 2%.

Calculation of electric field gradients with the exact two-component (X2C) quasi-relativistic method and its local approximations

When calculating electric field gradients (EFGs), relativistic and electron correlation effects are crucial for obtaining accurate results, and the commonly used density functional methods produce unsatisfactory results, especially for heavy elements and/or strongly correlated systems. In this work, a stand-alone program is presented, which enables calculation of EFGs from the molecular orbitals supplied by an external high accuracy quantum chemical calculation and includes relativistic effects through the exact two-component (X2C) formalism and efficient local approximations to it. Application to BiN and BiP molecules shows that a high precision can be achieved in the calculation of nuclear quadrupole coupling constants of 209Bi by combining advanced ab initio methods with the X2C approach. For seventeen iron compounds, the Mössbauer nuclear quadrupole splittings (NQS) of 57Fe calculated using a double-hybrid functional method are in very good agreement with the experimental values. It is shown that, for strongly correlated molecules, the double-hybrid functionals are much more accurate than the commonly used hybrid functionals. The computer program developed in this study furnishes a useful utility for obtaining EFGs and related nuclear properties with high accuracy.

Coordination of Mercury(II) in Water Promoted over Hydrolysis in Solvated Clusters [Hg(H2O)1-6](aq)2+: Insights from Relativistic Effects and Free Energy Analysis

Understanding the nature of the interaction between mercury(II) ions, Hg2+, and water molecules is crucial to describe the stability and chemical behavior of structures formed during solvation, as well as the conditions that favor the Hg2+ coordination or inducing water hydrolysis. In our study, we explored exhaustively the potential energy surface of Hg2+ with up to six water molecules. We analyzed electronic and Gibbs free energies, binding, and nuclear magnetic resonance parameters. We used the zeroth-order regular approximation Hamiltonian, including scalar and spin-orbit relativistic corrections for free energy calculations and geometry optimizations to explore the interplay between electron correlation and relativistic effects. We analyzed intermolecular interactions with energy decomposition analysis, quantum theory of atoms in molecules, and natural bond orbital. Additionally, we used the four-component Dirac Hamiltonian to compute solvent effect on the magnetic shielding and J-coupling constants. Our results revealed that the water hydrolysis by Hg2+ requires a minimum of three water molecules. We found that the interaction between Hg2+ and water molecules is an orbital interaction due to relativistic effects and the most stable structures are opened-shape clusters, reducing the number of oxygen-mercury contacts and maximizing the formation of hydrogen bonds among water molecules. In these types of clusters, Hg2+ promotes the water hydrolysis over coordination with oxygen atoms. However, when we considered the change associated with the transfer of a cluster from the ideal gas to a solvated system, our solvation free energy analysis revealed that closed-shape clusters are more favorable, maximizing the number of oxygen-mercury contacts and reducing the formation of hydrogen bonds among water molecules. This finding suggests that, under room conditions, the coordination of Hg2+ is more favorable than hydrolysis. Our results have significant implications for understanding Hg2+ behavior in water, helping to develop targeted strategies for mercury remediation and management, and contributing to advancements in the broader field of environmental chemistry.

High order relativistic corrections on the electric field gradient within the LRESC formalism

In this work, we present relativistic corrections to the electric field gradient (EFG) given by the Linear Response Elimination of the Small Component (LRESC) scheme at 1/c² order and including for the first time spin-dependent (SD) corrections at 1/c⁴ order. We show that these new terms improve the performance of LRESC as results with this methodology are very close to those calculated at the four-component Dirac-Hartree-Fock (4c-DHF) level. We assess the new corrections in BrY and AtY di-halogen (Y = F, Cl, Br, I, and At) and XZY bi-linear molecules (Z = Zn, Cd, and Hg; X, Y = F, Cl, Br, I, and At). At the 4c-DHF level, we analyze the contributions coming from the large and small components of the relativistic 4c wave function to the electronic part of EFG and compare them with the LRESC corrections to find their electronic origin. For the HgX₂ (X = Cl, Br, and I) subset, when the SD correcting terms are included, LRESC calculations match very well with 4c-DHF ones and those from the literature, with differences less than 1% for molecules containing up to three heavy atoms. We show that LRESC gives accurate values of EFG, allowing the analysis of the electronic origin of relativistic effects in terms of well-known nonrelativistic operators.

Free Molecule Studies by Perturbed γ-γ Angular Correlation: A New Path to Accurate Nuclear Quadrupole Moments

Accurate nuclear quadrupole moment values are essential as benchmarks for nuclear structure models and for the interpretation of experimentally determined nuclear quadrupole interactions in terms of electronic and molecular structure. Here, we present a novel route to such data by combining perturbed γ-γ angular correlation measurements on free small linear molecules, realized for the first time within this work, with state-of-the-art ab initio electronic structure calculations of the electric field gradient at the probe site. This approach, also feasible for a series of other cases, is applied to Hg and Cd halides, resulting in Q(^{199}Hg,5/2^{-})=+0.674(17)  b and Q(^{111}Cd,5/2^{+})=+0.664(7)  b.

Enhancing NMR Quantum Computation by Exploring Heavy Metal Complexes as Multiqubit Systems: A Theoretical Investigation

Assembled together with the most common qubits used in nuclear resonance magnetic (NMR) quantum computation experiments, spin-1/2 nuclei, such as ¹¹³Cd, ¹⁹⁹Hg, ¹²⁵Te, and ⁷⁷Se, could leverage the prospective scalable quantum computer architectures, enabling many and heteronuclear qubits for NMR quantum information processing (QIP) implementations. A computational design strategy for prescreening recently synthesized complexes of cadmium, mercury, tellurium, selenium, and phosphorus (called MRE complexes) as suitable qubit molecules for NMR QIP is reported. Chemical shifts and spin-spin coupling constants (SSCCs) in five MRE complexes were examined using the spin-orbit zeroth order regular approximation (ZORA) at the density functional theory level and the four-component relativistic Dirac-Kohn-Sham approach. In particular, the influence of different conformers, basis sets, exchange-correlation functionals, and methods to treat the relativistic as well as solvent effects were studied. The differences in the chemical shifts and SSCCs between different low energy conformers of the studied complexes were found to be very small. The TZ2P basis set was found to be the optimum choice for the studied chemical shifts, while the TZ2P-J basis set was the best for the couplings studied in this work. The PBE0 exchange-correlation functional exhibited the best performance for the studied MRE complexes. The addition of solvent effects has not improved on the gas phase results in comparison to the experiment, with the exception of the phosphorus chemical shift. The use of MRE complexes as qubit molecules for NMR QIP could face the challenges in single qubit control and multiqubit operations. They exhibit chemical shifts appropriately dispersed, allowing qubit addressability and exceptionally large spin-spin couplings, which could reduce the time of quantum gate operations and likely preserve the coherence.

Hyperfine Structure Constants on the Relativistic Coupled Cluster Level with Associated Uncertainties

Accurate predictions of hyperfine structure (HFS) constants are important in many areas of chemistry and physics, from the determination of nuclear electric and magnetic moments to benchmarking of new theoretical methods. We present a detailed investigation of the performance of the relativistic coupled cluster method for calculating HFS constants within the finite-field scheme. The two selected test systems are 133Cs and 137BaF. Special attention has been paid to construct a theoretical uncertainty estimate based on investigations on basis set, electron correlation and relativistic effects. The largest contribution to the uncertainty estimate comes from higher order correlation contributions. Our conservative uncertainty estimate for the calculated HFS constants is ∼5.5%, while the actual deviation of our results from experimental values is <1% in all cases.

First example of a high-level correlated calculation of the indirect spin-spin coupling constants involving tellurium: tellurophene and divinyl telluride

This paper documents the very first example of a high-level correlated calculation of spin-spin coupling constants involving tellurium taking into account relativistic effects, vibrational corrections and solvent effects for medium sized organotellurium molecules. The (125)Te-(1)H spin-spin coupling constants of tellurophene and divinyl telluride were calculated at the SOPPA and DFT levels, in good agreement with experimental data. A new full-electron basis set, av3z-J, for tellurium derived from the "relativistic" Dyall's basis set, dyall.av3z, and specifically optimized for the correlated calculations of spin-spin coupling constants involving tellurium was developed. The SOPPA method shows a much better performance compared to DFT, if relativistic effects calculated within the ZORA scheme are taken into account. Vibrational and solvent corrections are next to negligible, while conformational averaging is of prime importance in the calculation of (125)Te-(1)H spin-spin couplings. Based on the performed calculations at the SOPPA(CCSD) level, a marked stereospecificity of geminal and vicinal (125)Te-(1)H spin-spin coupling constants originating in the orientational lone pair effect of tellurium has been established, which opens a new guideline in organotellurium stereochemistry.