Extensive Ab Initio Research Articles

Vibronic coupling in the energetically six lowest electronic states of oxirane radical cation.

Multi-dimensional quantum mechanical simulations are carried out to understand the multi-state and multi-mode vibronic interactions in the first six low-lying viz., X̃2B1, Ã2A1, B̃2B2, C̃2A2, D̃2A1, and Ẽ2B1 electronic states of c-C2H4O·+. Vibronic coupling theory is applied to study interactions among electronic states using symmetry selection rules. A model 6 × 6 diabatic electronic Hamiltonian is constructed. The parameters of the diabatic Hamiltonian are estimated by performing extensive abinitio electronic structure calculations, using the EOM-IP-CCSD method. The nuclear dynamics calculations are performed with both time-independent and time-dependent quantum mechanical methods. The calculated vibronic structures of six electronic states are found to be in excellent agreement with the available experimental findings. Progressions found in the theoretical spectrum are assigned in terms of vibrational modes. It is found that extremely strong vibronic interactions among the X̃2B1-Ã2A1, B̃2B2-C̃2A2, and D̃2A1-Ẽ2B1 electronic states results into highly overlapping vibronic bands due to multiple multi-state conical intersections. The impact of associated nonadiabatic effects on the vibronic structure and dynamics of the mentioned electronic states is examined at length. Interesting comparison is made with the results obtained for the isomeric acetaldehyde radical cation.

Electronic density response of warm dense hydrogen on the nanoscale.

The properties of hydrogen at warm dense matter (WDM) conditions are of high importance for the understanding of astrophysical objects and technological applications such as inertial confinement fusion. In this work, we present extensive ab initio path integral Monte Carlo results for the electronic properties in the Coulomb potential of a fixed ionic configuration. This gives us unique insights into the complex interplay between the electronic localization around the protons with their density response to an external harmonic perturbation. We find qualitative agreement between our simulation data and a heuristic model based on the assumption of a local uniform electron gas model, but important trends are not captured by this simplification. In addition to being interesting in their own right, we are convinced that our results will be of high value for future projects, such as the rigorous benchmarking of approximate theories for the simulation of WDM, most notably density functional theory.

Analysing the dynamic structure of warm dense matter in the imaginary-time domain: theoretical models and simulations.

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [Nat. Commun. 13, 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density-density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Alkali Halide and MIBC Interaction at Typical Flotation Interfaces in Saline Water as Determined by Molecular Dynamics Simulations

The molecular structure of the liquid–vapor interfaces of aqueous solutions of alkali metal halides and methyl isobutyl carbinol (MIBC, (CH3)2CHCH2COCH3) is determined by using molecular dynamics simulations with polarizable force fields for the first time. The salts are chlorides, and iodides, some of which are found in raw and partially desalinated seawater increasingly used in flotation operations in regions affected by severe and prolonged drought. The density profiles at the interfaces show that all ions prefer the interface; however, with MIBC, non-polarizable ions, generally small ones, are increasingly pushed into the liquid bulk. A few ions of comparatively less ionic NaCl than KCl and CsCl, persist at the interface, consistent with spectroscopy observations. On the other hand, strongly polarizable ions such as I− always share the interface with MIBC. In the presence of chlorides, the frother chains at the interface stretch slightly more toward vapor than in freshwater; however, in the presence of iodides, the chains stretch so much that they become orthogonal to the interface, giving rise to a well-packed monolayer, which is the most effective configuration. The dominant water configurations at the interface are double donor and single donor, with hydrogen atoms pointing toward the liquid, consistent with studies with sum-frequency generation experiments and extensive ab initio simulations. This picture changes radically in the presence of MIBC and salts. Depending on the halide and MIBC concentration, the different molecular configurations at the interface lead to very different surface tensions. The structure and properties of these new salt-rich interfaces and their impact on the location and arrangement of frother molecules should serve the flotation practitioner, especially in the search for the best frother and dosing in poor-quality water.

The uniform electron gas at high temperatures: ab initio path integral Monte Carlo simulations and analytical theory

We present extensive new Ab initio path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) in the high-temperature regime, 8≤θ=kBT/EF≤128. This allows us to study the convergence of different properties towards the classical limit. In particular, we investigate the classical relation between the static structure factor S(q) and the static local field correction G(q), which is only fulfilled at low densities. Moreover, we compare our new results for the interaction energy to the parametrization of the UEG by Groth et al. (2017), which interpolates between PIMC results for θ≤8 and the Debye–Hückel limit, and to higher order analytical virial expansions. Finally, we consider the momentum distribution function n(q) and find an interaction-induced increase in the occupation of the zero-momentum state even for θ≳32. All PIMC data are freely available online, and can be used as input for improved parametrizations and as a rigorous benchmark for approximate methods.

Spin-resolved density response of the warm dense electron gas

We present extensive ab initio path integral Monte Carlo (PIMC) results for the spin-resolved density response of the uniform electron gas at warm dense matter conditions. This allows us to unambiguously assess the accuracy of previous theoretical approximations, thereby providing valuable insights for the future development of dielectric schemes. From a physical perspective, we observe a nontrivial manifestation of an effective electron-electron attraction that emerges in the spin-offdiagonal static density response function at strong coupling, ${r}_{s}\ensuremath{\gtrsim}5$. All PIMC results are freely available online and can be used to benchmark approximations and simulation schemes.

Nonlinear interaction of external perturbations in warm dense matter

AbstractWe present extensive new ab initio path integral Monte Carlo (PIMC) results for an electron gas at warm dense matter conditions that is subject to multiple harmonic perturbations. In addition to the previously investigated nonlinear effects at the original wave number (Dornheim et al., Phys. Rev. Lett., 2020, 125, 085001) and the excitation of higher harmonics (Dornheim et al., Phys. Rev. Res., 2021, 3, 033231), the presence of multiple external potentials leads to mode‐coupling effects, which constitute the dominant nonlinear effect and lead to a substantially more complicated density response compared to linear response theory. One possibility to estimate mode‐coupling effects from a PIMC simulation of the unperturbed system is given in terms of generalized imaginary‐time correlation functions that have been recently introduced by Dornheim et al. (J. Chem. Phys., 2021, 155, 054110). In addition, we extend our previous analytical theory of the nonlinear density response of the electron gas in terms of the static local field correction (Dornheim et al., 2020, Phys. Rev. Lett., 125, 235001), which allows for a highly accurate description of the PIMC results with negligible computational cost.

An accurate many-body expansion potential energy surface for AlH2 (22A') and quantum dynamics in Al(3P) + H2 (v0 = 0-3, j0 = 0, 2, 4, 6) collisions.

An accurate potential energy surface is constructed for the excited state of AlH2 by fitting extensive ab initio points calculated at the multi-reference configuration interaction level based on aug-cc-pV(Q+d)Z and aug-cc-pV(5+d)Z basis sets. All the calculated energies are corrected via the many-body expansion method and extrapolated to the complete basis set limit. The various topographic features of the new potential energy surface are investigated to demonstrate the correct behavior of Al(3P) + H2(X1Σg+) and AlH(a3Π) + H(2S) dissociation limits. By employing the time-dependent wave packet approach, the integral scattering cross-sections obtained from the Coriolis coupling calculation and the centrifugal sudden approximation, respectively, are compared in detail and show that the former has a higher effect on the reaction. Moreover, the thermal rate constants for Al(3P) + H2 (v0 = 0-3, j0 = 0, 2, 4, 6) in the temperature range of 0-5000 K are calculated, thereby providing insights into the influence of ro-vibrational quantum numbers on the thermal rate constants.

Density response of the warm dense electron gas beyond linear response theory: Excitation of harmonics

In a recent Letter, Dornheim et al. [PRL 125, 085001 (2020)] have investigated the nonlinear density response of the uniform electron gas in the warm dense matter regime. More specifically, they have studied the cubic response function at the first harmonic, which cannot be neglected in many situations of experimental relevance. In this work, we go one step further and study the full spectrum of excitations at the higher harmonics of the original perturbation based on extensive new ab initio path integral Monte Carlo (PIMC) simulations. We find that the dominant contribution to the density response beyond linear response theory is given by the quadratic response function at the second harmonic in the moderately nonlinear regime. Furthermore, we show that the nonlinear density response is highly sensitive to exchange-correlation effects, which makes it a potentially valuable new tool of diagnostics. To this end, we present a new theoretical description of the nonlinear electronic density response based on the recent effective static approximation to the local field correction [PRL 125, 235001 (2020)], which accurately reproduces our PIMC data with negligible computational cost.

Spectral control via multi-species effects in PW-class laser-ion acceleration

Laser-ion acceleration with ultra-short pulse, petawatt-class lasers is dominated by non-thermal, intra-pulse plasma dynamics. The presence of multiple ion species or multiple charge states in targets leads to characteristic modulations and even mono-energetic features, depending on the choice of target material. As spectral signatures of generated ion beams are frequently used to characterize underlying acceleration mechanisms, thermal, multi-fluid descriptions require revision for predictive capabilities and control in next-generation particle beam sources. We present an analytical model with explicit inter-species interactions, supported by extensive ab initio simulations. This enables us to derive important ensemble properties from the spectral distribution resulting from these multi-species effects for arbitrary mixtures. We further propose a potential experimental implementation with a novel cryogenic target, delivering jets with variable mixtures of hydrogen and deuterium. Free from contaminants and without strong influence of hardly controllable processes such as ionization dynamics, this would allow a systematic realization of our predictions for the multi-species effect.

Ab initio spectroscopy of the aluminum methylene (AlCH2) free radical.

Extensive ab initio investigations of the ground and electronic excited states of the AlCH2 free radical have been carried out in order to predict the spectroscopic properties of this, as yet, undetected species. Difficulties with erratic predictions of the ground state vibrational frequencies, both in the literature and in the present work, have been traced to serious broken-symmetry instabilities in the unrestricted Hartree-Fock orbitals at the ground state equilibrium geometry. The use of restricted open-shell Hartree-Fock or complete active space self consistent field orbitals avoids these problems and leads to consistent and realistic sets of vibrational frequencies for the ground state. Using the internally contracted multireference configuration interaction method with aug-cc-pV(T+d)Z basis sets, we have calculated the geometries, energies, dipole moments, and vibrational frequencies of eight electronic states of AlCH2 and AlCD2. In addition, we have generated Franck-Condon simulations of the expected vibronic structure of the Ã-X̃, B̃-X̃, C̃-X̃, and C̃-Ã band systems, which will be useful in searches for the electronic spectra of the radical. We have also simulated the expected rotational structure of the 0-0 absorption bands of these transitions at modest resolution under supersonic expansion cooled conditions. Our conclusion is that if AlCH2 can be generated in sufficient concentrations in the gas phase, it is most likely detectable through the B̃2A2-X̃2B1 or C̃2A1-X̃2B1 electronic transitions at 515 nm and 372 nm, respectively. Both band systems have vibrational and rotational signatures, even at modest resolution, that are diagnostic of the aluminum methylene free radical.

Ab initio studies of propene oxide formation at gold nanocatalysts supported on anatase-TiO2

Small gold particles supported on TiO2 are promising catalysts for the direct epoxidation of propene. Some experimental studies suggest a critical role of the oxide/gold interfacial region, where a bidentate propoxy species is formed during the reaction. However, it is not yet clear that such complex is a true reaction intermediate, and not a mere spectator or even a deactivating species. In order to clarify this issue, we have performed extensive ab initio simulations of the propene oxide formation on a model catalyst formed by small gold clusters (Au4, Au5 and Au8) supported on the anatase-TiO2(101) surface. In a variety of situations (with or without the supported gold cluster, with a varying amount of coadsorbed water or hydroxyl groups, etc.) we have investigated the stability of the bidentate propoxy species adsorbed at TiO2, as well as the barrier for its desorption into propene oxide. In all cases, the bidentate propoxy bonded at two surface Ti sites is highly stable, with a binding energy of the order of 1.5 eV. Such stability leads to very high energy barriers for propene oxide formation, of almost 2 eV. This suggests that such species is likely to be an spectator. Alternative reaction pathways, with the propoxy species formed in direct contact with the gold catalyst, are a much more probable candidate for a true reaction intermediate towards propene oxide formation.

Phosphonate-assisted tetranuclear lanthanide assemblies: observation of the toroidic ground state in the TbIII analogue.

The reaction of LnCl3·6H2O with a multidentate flexible Schiff base ligand (LH4), H2O3PtBu and trifluoroacetic acid (tfaH) afforded a series of homometallic tetranuclear complexes, [Ln4(LH2)2(O3PtBu)2(μ2-η1η1tfa)2][2Cl] (Ln = DyIII (1), TbIII (2) and GdIII (3)). The tetranuclear lanthanide core contains two structurally different lanthanide centres, one being in a distorted trigonal dodecahedron geometry and the other in a distorted trigonal prism. Complexes 1-3 were investigated via direct and alternating current (DC and AC) magnetic susceptibility measurements. Only 1 revealed a weak single-molecule magnet (SMM) behaviour. Alternating current (ac) magnetic susceptibility measurements on 1 reveal a frequency-dependent out-of-phase signal. However, the absence of distinct maxima in the χ'' peak (within the temperature/frequency range of our experiments) prevented deduction of the experimental energy barrier for magnetization reversal (Ueff) and the relaxation time. We have carried out extensive ab initio (CASSCF + RASSI-SO + SINGLE_ANISO + POLY_ANISO) calculations on complexes 1-2 to gain deeper insights into the nature of magnetic anisotropy. Our calculations yielded only one exchange coupling parameter between the two LnIII centres bridged by the ligand (neglecting the exchange between the LnIII centres that are not proximal wrt each other). All the extracted J values indicate a weakly antiferromagnetic coupling between the metal centres (J = -0.025 cm-1 for 1 and J = -0.015 cm-1 for 2). Calculated exchange coupled Ucal values of ∼5 and ∼1 cm-1 in 1 and 2 respectively nicely corroborated the experimental observations regarding weak and no SMM characteristics. Our calculations indicated the presence of a net single-molecule toroidal (SMT) behaviour in complex 2. On the other hand, fitting the magnetic data (susceptibility and magnetization) in the isotropic cluster 3 revealed weak AFM exchange couplings of J1 = 0.025 cm-1 and J2 = -0.020 cm-1 which are consistent with those for GdIII ions.

Nanoconfined Water within Graphene Slit Pores Adopts Distinct Confinement-Dependent Regimes.

In view of the increasing importance of nanoconfined aqueous solutions for various technological applications, it has become necessary to understand how strong confinement affects the properties of water at the level of molecular and even electronic structure. By performing extensive ab initio simulations of two-dimensionally nanoconfined water lamellae between graphene sheets subject to different interlayer spacings, we find new regimes at interlayer distances of 10 Å and less where water can be described neither to behave like interfacial water nor to be bulklike at the level of its H-bonding characteristics and electronic structure properties. It is expected that this finding will offer new opportunities to tune both diffusive and reactive processes taking place in aqueous environments that are strongly confined by chemically inert hard walls.

Material Properties for the Interiors of Massive Giant Planets and Brown Dwarfs

Abstract We present thermodynamic material and transport properties for the extreme conditions prevalent in the interiors of massive giant planets and brown dwarfs. They are obtained from extensive ab initio simulations of hydrogen–helium mixtures along the isentropes of three representative objects. In particular, we determine the heat capacities, the thermal expansion coefficient, the isothermal compressibility, and the sound velocity. Important transport properties such as the electrical and thermal conductivity, opacity, and shear viscosity are also calculated. Further results for associated quantities, including magnetic and thermal diffusivity, kinematic shear viscosity, as well as the static Love number k 2 and the equidistance, are presented. In comparison to Jupiter-mass planets, the behavior inside massive giant planets and brown dwarfs is stronger dominated by degenerate matter. We discuss the implications on possible dynamics and magnetic fields of those massive objects. The consistent data set compiled here may serve as a starting point to obtain material and transport properties for other substellar H–He objects with masses above one Jovian mass and finally may be used as input for dynamo simulations.

Acid Dissociation in HCl-Water Clusters is Temperature Dependent and Cannot be Detected Based on Dipole Moments.

The dissociation of acids in aqueous environments at low temperatures in the presence of a limited amount of water is underlying a wealth of processes from atmospheric to interstellar science. For the paradigmatic case of HCl(H_{2}O)_{n} clusters, our extensive abinitio path integral simulations quantify in terms of free energy differences and barriers that n=4 water molecules are indeed required to dissociate HCl at low temperatures. Increasing the temperature, however, reverses the process and thus counteracts dissociation by fluctuation-driven recombination. The size of the electric dipole moment is shown to not correlate with the acid being in its dissociated or molecular state, thus rendering its measurement as a function of n unable to detect the dissociation transition.

Systematic Study on Hydrated Arginine: Clear Theoretical Evidence for the Canonical-to-Zwitterionic Structure Transition.

Extensive ab initio investigations have been performed to characterize the stable conformers of hydrated arginine (Arg-H2O). Many new low-energy canonical Arg-H2O conformers were identified and they are more stable than previous results. The large energy differences (more than 5.00 kcal mol-1) between the canonical and zwitterionic Arg-H2O isomers calculated by the composite CBS-QB3 method confirmed the dominance of the zwitterions. The micro effects of corrections of the zero-point energy and the basis set superposition error on the stability of hydrated isomers were carefully examined for the first time. The infrared (IR) spectra were simulated at a recommended temperature and the notable spectral differences enable the unambiguous identification of the different hydrated forms. Further transition state calculations revealed that the canonical Arg-H2O can be transformed to the zwitterions using the amino group as a bridge. Our study thus shows valuable insights into the hydration of large flexible molecules and provides solid theoretical evidence for the canonical-to-zwitterionic structure transition of hydrated arginine.

Identification of the protonation site of gaseous triglycine: the cis-peptide bond conformation as the global minimum.

Extensive ab initio investigations have been performed to characterize stable conformers of protonated triglycine (GGGH) in the gas phase. Calculations using the composite CBS-QB3 method confirmed that the most favorable site of protonation on triglycine at 298 K is still the traditional amino nitrogen, rather than the more-recently reported amide oxygen. Furthermore, a non-proline cis-peptide bond conformer is identified for the first time as the global minimum of GGGH. Further transition state calculations considering the temperature effects explained why the previous experimental infrared multiple photon dissociation (IRMPD) spectrum contains a combination of two local minima, rather than a global one. First-principles simulations have been performed for near-edge X-ray absorption fine-structure (NEXAFS) spectra and X-ray photoelectron spectra (XPS) at the C, N and O K-edges to identify the notable spectral differences that enable the unambiguous identification of different protonated forms. The calculated proton affinity (PA) and gas basicity (GB) of triglycine are in excellent agreement with the experimental values. Our study thus provides valuable insights into the protonation of short peptides and illustrates the competition between cis and trans peptide bonds.

Chemistry in nanoconfined water.

Nanoconfined liquids have extremely different properties from the bulk, which profoundly affects chemical reactions taking place in nanosolvation. Here, we present extensive ab initio simulations of a vast set of chemical reactions within a water lamella that is nanoconfined by mineral surfaces, which might be relevant to prebiotic peptide formation in aqueous environments. Our results disclose a rich interplay of distinct effects, from steric factors typical of reactions occurring in small spaces to a charge-stabilization effect in nanoconfined water at extreme conditions similar to that observed in bulk water when changing from extreme to ambient conditions. These effects are found to modify significantly not only the energetics but also the mechanisms of reactions happening in nanoconfined water in comparison to the corresponding bulk regime.

Ab initio study of intrinsic profiles of liquid metals and their reflectivity

The free surfaces of liquid metals are known to exhibit a stratified profile that, in favourable cases, shows up in experiments as a peak in the ratio between the reflectivity function and that of an ideal step-like profile. This peak is located at a wave-vector related to the distance between the layers of the profile. In fact the surface roughness produced by thermally induced capillary waves causes a depletion of the previous so called intrinsic reflectivity by a damping factor that may hinder the observation of the peak. The behaviour of the intrinsic reflectivity below the layering peak is however far from being universal, with systems as Ga or In where the reflectiviy falls uniformly towards the q → 0 value, others like Sn or Bi where a shoulder appears at intermediate wavevectors, and others like Hg which show a minimum. We have performed extensive ab initio simulations of the free liquid surfaces of Bi, Pb and Hg, that yield direct information on the structure of the profiles and found that the macroscopic capillary wave theory usually employed in order to remove the capillary wave components fails badly in some cases for the typical sample sizes affordable in ab initio simulations. However, a microscopic method for the determination of the intrinsic profile is shown to be succesful in obtaining meaningful intrinsic profiles and corresponding reflectivities which reproduce correctly the qualitative behaviour observed experimentally.