Hartree-Fock Research Articles

Energy levels, oscillator strengths, transition probabilities and lifetimes of the O-like Cl X ion

Energy levels, lifetimes, oscillator strengths and transition probabilities for the multicharged oxygen like Cl X ion have been calculated with the configuration expansion: 2s22p4, 2p6, 2s22p33p, 2s2p5, 2s22p33s, 2s22p34s, 2s22p33d, 2s22p34d, 2s2p43s and 2s2p44s. We used two methods in the calculations: the Hartree–Fock pseudo-relativistic approach and the Thomas–Fermi–Dirac–Amaldi potential approach using the Cowan and the AUTOSTRUCTURE atomic structure codes respectively. Results have been compared with available experimental data from the National Institute of Standards and Technology Atomic Spectra Database (NIST-ASD) and with results from other available calculations.There is a great lack of atomic structure data of Cl X and obtained new data will be important for physics applications and astrophysical modeling.

Studying shape phase transition in even-proton nuclei via mass measurements

The transition between the spherical and deformed shapes of the nuclei (shape phase transition) is an interesting property of nuclei that enriches our understanding of the nuclear structure. One method for studying it is through measurement of a wide range of atomic masses and derivation of the two-neutron separation energies (S2n). Using a multiple-reflection time-of-flight mass spectrometer at the FRS Ion Catcher, GSI, atomic masses in the region of A ∼ 150 were measured. This allowed an wide study of phase shape transitions in this region of the nuclear chart. In addition, comparison with the Skyrme Hartree–Fock plus BCS approach are done; showing a remarkable good agreement between theory and experiment.

Excited State-Specific CASSCF Theory for the Torsion of Ethylene.

State-specific complete active space self-consistent field (SS-CASSCF) theory has emerged as a promising route to accurately predict electronically excited energy surfaces away from molecular equilibria. However, its accuracy and practicality for chemical systems of photochemical interest have yet to be fully determined. We investigate the performance of the SS-CASSCF theory for the low-lying ground and excited states in the double bond rotation of ethylene. We show that state-specific approximations with a minimal (2e,2o) active space provide comparable accuracy to state-averaged calculations with much larger active spaces, while optimizing the orbitals for each excited state significantly improves the spatial diffusivity of the wave function. However, the incorrect ordering of state-specific solutions causes excited state solutions to coalesce and disappear, creating unphysical discontinuities in the potential energy surface. Our findings highlight the theoretical challenges that must be overcome to realize practical applications of state-specific electronic structure theory for computational photochemistry.

Nature of partial sigma bond.

This study investigates the formation of partial sigma (σ) covalent bonds in experimentally synthesizable biradicals formed from hydrogenated and fluorinated C8, C20, and C60 cage structures, by assessing their stability, geometry, and bonding character in singlet and triplet states using restricted B3LYP-D3/6-31+G(d,p) theory, natural bond orbital (NBO) analysis, and complete active space self-consistent field (CASSCF) method. The results show that these partial σCC bonds have Wiberg bond orders of 0.38 to 0.48 and bond lengths ranging from 2.62 Å to 5.93 Å. Cage size influences the characteristics of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), with electrons favoring more antibonding orbitals in smaller cages where electrons reside more on the exterior of the cage and favoring bonding orbitals in larger ones where electrons are more in the interior. Fluorination enhances electron density on bonding orbitals. The analysis further clarified that the differentiation between antibonding and bonding features of HOMOs and LUMOs extends beyond merely electron transfer from s- to p-atomic orbitals, also noting possible interactions of the same symmetry repel. The study also introduces hyperconjugation from α-position CH bonds as a factor in stabilizing partial σ-bond formation. The results also caution against the use of broken symmetry methodology in unrestricted SCF wavefunctions for biradicals, such as those in this study as it may cause large spin contamination and thus errors in the calculated electronic properties results.

Numerical Study of Two Opposing Weak Polyelectrolyte Brushes by the Self-consistent Field Theory

Numerical Study of Two Opposing Weak Polyelectrolyte Brushes by the Self-consistent Field Theory

Superexchange Mechanism in Coupled Triangulenes Forming Spin-1 Chains.

We show that the origin of the antiferromagnetic coupling in spin-1 triangulene chains, which were recently synthesized and measured by Mishra et al. ( Nature 2021, 598, 287-292), originates from a superexchange mechanism. This process, mediated by intertriangulene states, opens the possibility to control parameters in the effective bilinear-biquadratic spin model. We start from the derivation of an effective tight-binding model for triangulene chains using a combination of tight-binding and Hartree-Fock methods fitted to hybrid density functional theory results. Next, correlation effects are investigated within the configuration interaction method. Our low-energy many-body spectrum for NTr = 2 and NTr = 4 triangulene chains agree well with the bilinear-biquadratic spin-1 chain antiferromagnetic model when indirect coupling processes and superexchange coupling between triangulene spins are taken into account.

VSCF/VCI theory based on the Podolsky Hamiltonian.

While the vibrational spectra of semi-rigid molecules can be computed on approaches relying on the Watson Hamiltonian, floppy molecules or molecular clusters are better described by Hamiltonians, which are capable of dealing with any curvilinear coordinates. It is the kinetic energy operator (KEO) of these Hamiltonians, which render the correlated calculations relying on them rather costly. Novel implementation of vibrational self-consistent field theory and vibrational configuration interaction theory on the basis of the Podolsky Hamiltonian are reported, in which the inverse of the metric tensor, i.e., the G matrix, is represented by an n-mode expansion expressed in terms of polynomials. An analysis of the importance of the individual terms of the KEO with respect to the truncation orders of the n-mode expansion is provided. Benchmark calculations have been performed for the cis-HOPO and methanimine, H2CNH, molecules and are compared to experimental data and to calculations based on the Watson Hamiltonian and the internal coordinate path Hamiltonian.

Predicting the One-Particle Density Matrix with Machine Learning.

Two of the most widely used electronic-structure theory methods, namely, Hartree-Fock and Kohn-Sham density functional theory, require the iterative solution of a set of Schrödinger-like equations. The speed of convergence of such a process depends on the complexity of the system under investigation, the self-consistent-field algorithm employed, and the initial guess for the density matrix. An initial density matrix close to the ground-state matrix will effectively allow one to cut out many of the self-consistent steps necessary to achieve convergence. Here, we predict the density matrix of Kohn-Sham density functional theory by constructing a neural network that uses only the atomic positions as information. Such a neural network provides an initial guess for the density matrix far superior to that of any other recipes available. Furthermore, the quality of such a neural-network density matrix is good enough for the evaluation of interatomic forces. This allows us to run accelerated ab initio molecular dynamics with little to no self-consistent steps.

Vibrational Self-Consistent Field (VSCF) and Post-VSCF Method Calculations Combined with the Reference Interaction Site Model Self-Consistent Field Method Coupled with the Constrained Spatial Electron Density Distribution: Applications to NaHCOO in Aqueous Phase.

Investigating vibrational behavior in solution is crucial for understanding molecular dynamics within a solvent environment. Notably, the analysis of Raman spectra for molecules in solution is important owing to its ability to unveil intricate solute-solvent interactions. Previous studies have effectively employed frequency calculations utilizing the reference interaction site model self-consistent field method in conjunction with constrained spatial electron density distribution (RISM-SCF-cSED) to understand molecular vibrations in solution, primarily focusing on fundamental vibrational modes. However, the oversight of overtones and combination tones in these studies prompted us to combine the vibrational self-consistent field (VSCF) and vibrational second-order Mo̷ller-Plesset perturbation (VMP2) methods with RISM-SCF-cSED to address these aspects theoretically. Illustrating the efficacy of this integrated approach, we computed the Raman spectra of sodium formate (NaHCOO) in water, revealing the necessity of accounting for molecular anharmonicity in solution vibrational analysis. Our findings underscore the potency of VSCF and VMP2 in conjunction with RISM-SCF-cSED as a robust theoretical framework for such calculations.

From many-body quantum dynamics to the Hartree–Fock and Vlasov equations with singular potentials

In a combined mean-field and semiclassical regime, we consider the time evolution of N fermions interacting through singular pair interaction potentials of the form \pm |x-y|^{-a} , which includes the Coulomb and gravitational interactions. We prove that the many-body dynamics of mixed states are well approximated by solutions of the Hartree–Fock and Vlasov equations in terms of Schatten norms. The errors in these approximations are expressed in terms of the expected number of particles, N , and the Planck constant, h . For cases where a\in(0,1/2) , we obtain local-in-time results when N^{-1/2}\ll h \leq N^{-1/3} . Notably, this leads to the derivation of the Vlasov equation with singular potentials. For cases where a\in[1/2,1] , our results hold only within a small time scale or require an N -dependent cut-off. A fundamental ingredient in our analysis is the propagation of regularity for solutions to the Hartree–Fock equation uniformly in the Planck constant, which holds for a\in(0,1] .

MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces.

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals.

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

Ground state properties of the screened helium atom under harmonic confinement

In this contribution, we show electronic property results for the ground state of the helium atom screened by a Debye plasma under an isotropic harmonic confinement as a function of the screening length (λD) and the angular frequency (ω0). The associate Hamiltonian is solved by expanding the wave function in the Hylleraas basis to incorporate the electron correlation in a simple manner. Furthermore, in order to compare with other findings the Hartree–Fock approach is implemented using a numerical grid method. The main results show that the total energy reaches the ionization threshold by the action of the electron coupling with the plasma medium and the strength of the harmonic confinement. The electron–electron average distance shows the compression of the system as the Debye plasma and the harmonic interactions become strongest (low λD and large ω0 values). We found acceptable results for the electron correlation energy when separate interactions are considered. In the absence of both interactions, the findings are in excellent agreement with other theoretical results.

Automatic State Interaction with Large Localized Active Spaces for Multimetallic Systems.

The localized active space self-consistent field method factorizes a complete active space wave function into an antisymmetrized product of localized active space wave function fragments. Correlation between fragments is then reintroduced through localized active space state interaction (LASSI), in which the Hamiltonian is diagonalized in a model space of LAS states. However, the optimal procedure for defining the LAS fragments and LASSI model space is unknown. We here present an automated framework to explore systematically convergent sets of model spaces, which we call LASSI[r, q]. This method requires the user to select only r, the number of electron hops from one fragment to another, and q, the number of fragment basis functions per Hilbert space, which converges to CASCI in the limit of r, q → ∞. Numerical tests of this method on the trimetal oxo-centered complexes [Fe(III)Al(III)Fe(II)(μ3-O)(HCOO)6] and [Fe(III)2Fe(II)(μ3-O)(HCOO)6] show efficient convergence to the CASCI limit with 4-10 orders of magnitude fewer states than CASCI.

Assessment of Iron-Based Spin-Orbit Coupling Effects in Pt-Fe Heterobimetallic Lantern Complexes via 57Fe Mössbauer Spectroscopy.

A series of heterobimetallic lantern complexes, [PtFe(SOCR)4(pyX)] where R = Me, X = H (1), X = NH2 (2), X = SMe (3); R = Ph, X = H (4), X = NH2 (5), X = SMe (6), have been prepared and characterized spectroscopically. Compounds 1, 4, and 5 are reported herein for the first time. The high-spin iron(II) sites of 1-6 have been investigated using 57Fe Mössbauer spectroscopy. Although the isomer shift of these species is nearly identical, their quadrupole splitting exhibits a much larger variation. Moreover, the zero-field Mössbauer spectra of 3-5 show surprising changes over time which are likely indicative of small structural distortions. The field dependent Mössbauer study of 1 and 6 revealed a zero field splitting (ZFS) characterized by a relatively large and positive D value. The combined Density Functional Theory (DFT) and ab initio Complete Active Space Self-Consistent Field (CASSCF) investigation of 1-6 indicates that their ground state is best described using a linear combination of {|xz⟩, |yz⟩} states. Our theoretical analysis suggests that the ZFSs and magnitude of the quadrupole splitting of 1-6 are determined by the spin-orbit coupling of the three lowest orbital states which have a T2g parentage.

Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes.

Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.

Pseudo-Spin Symmetry and the Hints for Unstable and Superheavy Nuclei

The pseudo-spin symmetry (PSS) provides an important angle to understand nuclear microscopic structure and the novel phenomena found in unstable nuclei. The relativistic Hartree–Fock (RHF) theory, that takes the important degrees of freedom associated with the π-meson and ρ-tensor (ρ-T) couplings into account, provides an appropriate description of the PSS restoration in realistic nuclei, particularly for the pseudo-spin (PS) doublets with high angular momenta (l˜). The investigations of the PSS within the RHF theory are recalled in this paper by focusing on the effects of the Fock terms. Aiming at common artificial shell closures appearing in previous relativistic mean-field calculations, the mechanism responsible for the PSS restoration of high-l˜ orbits is stressed, revealing the manifestation of nuclear in-medium effects on the PSS, and thus, providing qualitative guidance on modeling the in-medium balance between nuclear attractions and repulsions. Moreover, the essential role played by the ρ-T coupling, that contributes mainly via the Fock terms, is introduced as combined with the relations between the PSS and various nuclear phenomena, including the shell structure and the evolution, novel halo and bubble-like phenomena, and the superheavy magicity. As the consequences of the nuclear force in complicated nuclear many-body systems, the PSS itself and the mechanism therein can not only deepen our understanding of nuclear microscopic structure and relevant phenomena, but also provide special insight into the nature of the nuclear force, which can further enrich our knowledge of nuclear physics.

Limitations of mean-field approximations in describing shift-current and injection-current in materials

We theoretically and computationally investigate bulk photovoltaic effects, with a specific focus on shift-current and injection-current. Initially, we perform a numerical analysis of the direct current (dc) induced by a laser pulse with a one-dimensional model, utilizing mean-field theories such as time-dependent Hartree–Fock and time-dependent Hartree methods. Our numerical results, obtained with mean-field theories, reveal that the dc component of the current, as a second-order nonlinear effect, exists even after irradiation with linearly polarized light as a second-order nonlinear effect, indicating the generation of injection-current. Conversely, when we employ the independent-particle approximation, no injection-current is generated by linearly polarized light. To develop the microscopic understanding of injection-current within the mean-field approximation, we further analyze the dc component of the current with the perturbation theory, employing the mean-field approximations, the independent-particle approximation, and the exact solution of the many-body Schrödinger equation. The perturbation analysis clarifies that the injection-current induced by linearly polarized light under the mean-field approximations is an artifact caused by population imbalance, created through quantum interference from unphysical self-excitation pathways. Therefore, investigation of many-body effects on the bulk photovoltaic effects have to be carefully conducted in mean-field schemes due to potential contamination by unphysical dc current. Additionally, we perform the first-principles electron dynamics calculation for BaTiO3 based on the time-dependent density functional theory, and we confirm that the above findings from the one-dimensional model calculation and the perturbation analysis apply to realistic systems. Published by the American Physical Society 2024

Nuclear Matter Equation of State in the Brueckner–Hartree–Fock Approach and Standard Skyrme Energy Density Functionals

The equation of state of asymmetric nuclear matter as well as the neutron and proton effective masses and their partial-wave and spin–isospin decomposition are analyzed within the Brueckner–Hartree–Fock approach. Theoretical uncertainties for all these quantities are estimated by using several phase-shift-equivalent nucleon–nucleon forces together with two types of three-nucleon forces, phenomenological and microscopic. It is shown that the choice of the three-nucleon force plays an important role above saturation density, leading to different density dependencies of the energy per particle. These results are compared to the standard form of the Skyrme energy density functional, and we find that it is not possible to reproduce the BHF predictions in the (S,T) channels in symmetric and neutron matter above saturation density, already at the level of the two-body interaction, and even more including the three-body interaction.

Quantum chemical, vibrational, microhardness, thermal and antibacterial studies on polymorphic drug crystal: Sulfanilamide (β-form)

The widely known pharmaceutical drug compound, Sulfanilamide, crystallized in the β-form of its polymorphic phases by slow evaporation of the parent solution at room temperature. The single crystal X-ray diffraction study provided information about the cell parameters, confirming the β-phase of the compound. Optimization of the molecule by quantum chemical computations, viz., Hartree-Fock (HF) Theory and Density Functional Theory (DFT) using the B3LYP function with 6–311++G(d,p) basis set for both in restricted methods were carried out. The Frontier Molecular Orbital studies gave the HOMO and LUMO values, and Mulliken charge analysis showed the possible charge delocalization and donor-acceptor sites. Vibrational bands of the functional groups present in the β-Sulfanilamide compound were identified and investigated using FT-IR and FT-Raman spectroscopy in the range of 4000–400 cm−1. The grown crystal was transparent in the entire UV region with a lower cut off wavelength of 241 nm as investigated by the UV–Vis-NIR analysis. Thermogravimetric and differential thermal analyses of the compound were carried out and the results were plotted as TGA/DTA graph. Microhardness tests were used to determine the mechanical parameters of this pharmaceutical compound. Sulfanilamide, a well-known antibacterial agent and used to treat a variety of bacterial infections, was tested using the zone inhibition method and the antibacterial activity was confirmed.