Search

Ab initio calculations of electric field gradients (EFGs) in molecular crystals have advanced significantly due to the gauge including projector augmented wave (GIPAW) formalism, which accounts for the infinite periodicity in crystals. However, theoretical accuracies still lag behind experimental ones, making it challenging to distinguish experimentally distinguishable similar structures, a deficiency largely attributed to the limitation of GIPAW codes to generalized gradient approximation (GGA) density functional theory (DFT) functionals. In this study, we investigate whether hybrid DFT functionals can enhance the EFG calculation accuracy and the associated geometry optimization. Using the many-body expansion method, we focus on nitrogen EFGs in amino acids with complex H-bonding, which are often poorly described with GGA functionals. Our results show that both functionals provide highly accurate calculations that surpass current studies and approach experimental precision. The accuracies are also almost three times higher than available GIPAW/GGA calculations in the literature. However, we show that this difference is not due to the GGA functional but rather due to the improper selection of the nitrogen quadrupole moment. [ABSTRACT FROM AUTHOR]

We present an experimental study of the magnetoelectric coupling in the Fe-Ga/Pb[(Mg 1 / 3 Nb 2 / 3 )O 3 ] 0.68 -[PbTiO 3 ] 0.31 thin-film multiferroic composite using x-ray magnetic circular dichroism and ferromagnetic resonance (FMR). Our measurements show evidence for a charge-mediated coupling mechanism, suggested by the asymmetric magnetic remanence (M r e m ) behavior under opposite electric fields (± E) and the asymmetric resonance field (H r) in the FMR measurements. Also, the FMR measurements reveal a perpendicular magnetic anisotropy that can be related to an interface charge effect and it is tunable by the E field. Ab initio calculations support the existence of a charge-mediated coupling at the Fe–Ga/PMN-PT interface. [ABSTRACT FROM AUTHOR]

Peptoids (N-substituted glycines) are a class of sequence-defined synthetic peptidomimetic polymers with applications including drug delivery, catalysis, and biomimicry. Classical molecular simulations have been used to predict and understand the conformational dynamics of single chains and their self-assembly into morphologies including sheets, tubes, spheres, and fibrils. The CGenFF-NTOID model based on the CHARMM General Force Field has demonstrated success in accurate all-atom molecular modeling of peptoid structure and thermodynamics. Extension of this force field to new peptoid side chains has historically required reparameterization of side chain bonded interactions against ab initio data. This fitting protocol improves the accuracy of the force field but is also burdensome and precludes modular extensibility of the model to arbitrary peptoid sequences. In this work, we develop and demonstrate a Modular Side Chain CGenFF-NTOID (MoSiC-CGenFF-NTOID) as an extension of CGenFF-NTOID employing a modular decomposition of the peptoid backbone and side chain parameterizations, wherein arbitrary side chains within the large family of substituted methyl groups (i.e., –CH3, –CH2R, –CHRR′, and –CRR′R″) are directly ported from CGenFF. We validate this approach against ab initio calculations and experimental data to develop a MoSiC-CGenFF-NTOID model for all 20 natural amino acid side chains along with 13 commonly used synthetic side chains and present an extensible paradigm to efficiently determine whether a novel side chain can be directly incorporated into the model or whether refitting of the CGenFF parameters is warranted. We make the model freely available to the community along with a tool to perform automated initial structure generation. [ABSTRACT FROM AUTHOR]

Precise material design and surface engineering play a crucial role in enhancing the performance of optoelectronic devices. These efforts are undertaken to particularly control the optoelectronic properties and regulate charge carrier dynamics at the surface and interface. In this study, we used ultrafast scanning electron microscopy (USEM), which is a powerful and highly sensitive surface tool that provides unique information about the photoactive charge dynamics of material surfaces selectively and spontaneously in real time and space in high spatial and temporal resolution. Here, time-resolved images of CdTe (110), CdSe (100), GaAs (110), and other semiconductors revealed that the presence of oxide layers on the surface of materials leads to an increase in the work function (WF) and trap state densities upon optical excitation, leading to the formation of dark image contrast in USEM experiments. These findings were further supported by ab initio calculations, which confirmed the reliability of the observed changes in the excited surface WFs. Besides enhancing our understanding of surface charge dynamics, it also offers valuable insights into manipulating these properties. This study paves the way for precise control and the potential to design highly efficient light-harvesting devices. [ABSTRACT FROM AUTHOR]

This research paper presents a new fundamental approach for evaluating accurate ab initio quartic, sextic, and octic centrifugal distortion parameters of A-reduced rotational effective Hamiltonians of asymmetric top molecules. In this framework, the original Watson Hamiltonian, expanded up to sextic terms of kinetic and potential energies, is subjected to a series of vibrational and rotational operator unitary transformations, leading to reduced Watson effective Hamiltonians for the equilibrium configuration, ground state, and weakly perturbed vibrationally excited states. The proposed scheme is based on a numerical-analytic implementation of the sixth-order Van Vleck operator perturbation theory with the systematic normal ordering of vibrational rising and lowering operators (a†, a) and cylindrical angular momentum operators (Jz, J+, J−). The efficiency of the developed theoretical model is demonstrated by the juxtaposition of predicted centrifugal distortion parameters for several three to eight atomic molecules, including H2S, CH2O, C2H4, CH2D2, CH2F2, CH2Cl2, and B2H6, using the coupled-cluster single double triple/quadruple-ζ level of quantum chemistry. In comparison with the values derived using the customary analytic expressions, the calculated quartic and sextic parameters may improve by an order of magnitude in the fourth and sixth orders, respectively, reaching an accuracy of about 1%. Predicted octic constants can serve as an excellent starting point for fitting to experimental spectra. [ABSTRACT FROM AUTHOR]

Sodium bismuth titanate (NBT) reveals a rhombohedral (R3c) phase at room temperature. Ferroelectricity reduces with the advent of a tetragonal (P4bm) phase at the depolarization temperature, Td ∼ 456 K. AC conductivity (σac) studies exposed a small-to-large polaron transition at Td. Barrier energy (WH) was ∼1.60 eV at T < Td for the small polarons in the R3c phase, which drastically reduced to ∼0.043 eV with the advent of the P4bm phase for the large polarons for T > Td. This is associated with the sharp rise in conductivity for T > Td. Ab initio calculations consider the electronic distortion due to oxygen vacancies, which creates trap states in the band structure. The energy gap (ΔE) between the trap states and the conduction band was ∼1.4 eV (R3c) and ∼0.2 eV (P4bm). These values are comparable to the experimental WH. The P4bm phase is more distorted than the R3c phase from charge density and structural distortion calculations. This indicates the formation of large polarons in the P4bm phase, compared to that of small polarons in R3c. The formation energy of the polaron (Epolaron) was calculated from the structural distortion and electron localization energies. The P4bm phase shows lower Epolaron (−0.26 eV) than R3c (−0.36 eV), indicating higher conductivity for the P4bm phase. NBT was chemically modified by adding BCZT to validate the small to large polaronic crossover at Td. This is discussed in light of σac measurements. WH decreased with BCZT incorporation, thereby increasing the conductivity. This is a consequence of the increased lattice distortion due to BCZT incorporation. [ABSTRACT FROM AUTHOR]

The short-range order and intermediate-range order in GeO2 glass are investigated by molecular dynamics using machine-learning interatomic potential trained on ab initio calculation data and compared with the reverse Monte Carlo fitting of neutron diffraction data. To characterize the structural differences in each model, the total/partial structure factors, coordination number, ring size and shape distributions, and persistent homology analysis were performed. These results show that although the two approaches yield similar two-body correlations, they can lead to three-dimensional models with different short- and intermediate-range ordering. A clear difference was observed especially in the ring distributions; RMC models exhibit a broad distribution in the ring size distribution, while neural network potential molecular dynamics yield much narrower ring distributions. This confirms that the density functional approximation in the ab initio calculations determines the preferred network assembly more strictly than RMC with simple coordination constraints even when using multiple diffraction data. [ABSTRACT FROM AUTHOR]

We report a potential energy surface (PES) development for the F− + SiH3I system to study its gas-phase reactions through quasi-classical dynamics simulations. The PES is represented by a full-dimensional permutationally invariant polynomial fitted to composite coupled cluster energy points obtained at the ManyHF-[CCSD-F12b + BCCD(T) − BCCD]/aug-cc-pVTZ(-PP) level of theory. The development was automated by Robosurfer, which samples the configurational space, manages ab initio calculations, and iteratively extends the fitting set. When selecting the ab initio method, we address two types of electronic structure calculation issues: first, the gold standard CCSD(T)-F12b is prone to occasional breakdown due to the perturbative (T) contribution, whereas CCSD-F12b + BCCD(T) − BCCD, with the Brueckner (T) term, is more robust; second, the underlying Hartree–Fock calculation may not always converge to the global minimum, resulting in highly erroneous energies. To mitigate this, we employed ManyHF, configuring the Hartree–Fock calculations with multiple initial guess orbitals and selecting the solution with the lowest energy. According to the simulations, the title system exhibits exceptionally high and diverse reactivity. We observe two dominant product formations: SN2 and proton abstraction. Moreover, SiH2F− + HI, SiHFI− + H2, SiH2FI + H−, SiH2 + FHI−, SiH2 + HF + I−, and SiHF + H2 + I− formations are found at lower probabilities. We differentiated inversion and retention for SN2, both being significant throughout the entire collision energy range. Opacity- and excitation functions are reported, and the details of the atomistic dynamics are visually examined via trajectory animations. [ABSTRACT FROM AUTHOR]

We report on the application of the recently developed Fano-ADC(2,2) method to compute total and partial Auger decay widths of molecular core–hole states, including explicit evaluation of double Auger decay branching ratios. The method utilizes the fast-convergent intermediate state representation to construct many-electron wave functions and is readily applicable to atoms, molecules, and clusters. The ADC(2,2) scheme describes the initial and final states of the normal Auger decay consistently up to the second order of perturbation theory. In addition, excitations with two electrons in the continuum provide access to three-electron decay modes. The method yields decay widths and the Auger electron spectra in excellent agreement with the experiment, demonstrating the high accuracy of partial widths. The average relative error of double Auger decay branching ratios compared to available experimental data is about 30%, which should be evaluated as an excellent result considering the universality of the method, the complexity of the double decay process, and the neglection of nuclear motion in this study. [ABSTRACT FROM AUTHOR]

In order to understand the dispersion interactions between molecules and to provide information about the potential energy surface of geometry evolutions, NbN12− and N2·NbN12− complexes were investigated by using photoelectron spectroscopy and ab initio calculations. The experimental adiabatic detachment energy (ADE) and vertical detachment energy (VDE) of NbN12− were both measured to be 2.129 ± 0.030 eV. The experimental ADE and VDE of N2·NbN12− were measured to be 2.17 ± 0.05 and 2.23 ± 0.05 eV, respectively, which are slightly higher than those of NbN12−. The structures of NbN12−/0 were confirmed to be hexacoordinated octahedrons. The investigation of N2·NbN12− structures shows that it is stable for N2 to bind to the face or vertex site of octahedron NbN12−; the face-side-on structure has the lowest energy. The calculations based on symmetry-adapted perturbation theory suggest that the dispersion term is predominant and leads to the stability of N2·NbN12− complexes. [ABSTRACT FROM AUTHOR]

CrN-based alloy thin films are of interest as thermoelectric materials for energy harvesting. Ab initio calculations show that dilute alloying of CrN with 3 at. % W substituting Cr induce flat electronic bands and push the Fermi level EF into the conduction band while retaining dispersive Cr 3d bands. These features are conducive for both high electrical conductivity σ and high Seebeck coefficient α and, hence, a high thermoelectric power factor α2σ. To investigate this possibility, epitaxial CrWxNz films were grown on c-sapphire by dc-magnetron sputtering. However, even films with the lowest W content (x = 0.03) in our study contained metallic h-Cr2N, which is not conducive for a high α. Nevertheless, the films exhibit a sizeable power factor of α2σ ∼ 4.7 × 10−4 W m−1 K−2 due to high σ ∼ 700 S cm−1, and a moderate α ∼ − 25 μV/K. Increasing h-Cr2N fractions in the 0.03 < x ≤ 0.19 range monotonically increases σ, but severely diminishes α leading to two orders of magnitude decrease in α2σ. This trend continues with x > 0.19 due to W precipitation. These findings indicate that dilute W additions below its solubility limit in CrN are important for realizing a high thermoelectric power factor in CrWxNz films. [ABSTRACT FROM AUTHOR]

Anion photoelectron spectra of ThSO− and ThS2− were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO−, ThS2, and ThS2−, as well as the oxygen congeners ThO2 and ThO2−. The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO2 that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Overall the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol. [ABSTRACT FROM AUTHOR]

The Hugoniot is the equation of state of a shock-compressed material and is a key part of high-pressure physics. One way of calculating it is via the Hugoniostat that has significant computational advantages over direct calculation via non-equilibrium molecular dynamics. We introduce a number of improvements to the Hugoniostat, which significantly reduce the run time and the number of atoms required for converged results. Consequently, ab initio Hugoniot calculations are tractable. We illustrate the benefits through simple model potentials and with density functional theory calculations of argon. [ABSTRACT FROM AUTHOR]

Anion photoelectron spectra of ThSO− and ThS2− were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO−, ThS2, and ThS2−, as well as the oxygen congeners ThO2 and ThO2−. The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO2 that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Overall the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol. [ABSTRACT FROM AUTHOR]

Pressure-induced polymorphism has recently been demonstrated in several high entropy alloys. This offers a new window into the much-debated issue of phase selection and stability in these systems. Here, we examine the effect of cryogenic temperatures on the pressure-induced transition from face centered cubic to hexagonal close-packed structures of the prototype CoCrFeMnNi (Cantor) alloy. We observe a reduction in the critical pressure for the onset of the polymorphic transition as the temperature decreases, confirming the progressive stabilization of the hexagonal phase with decreasing temperature previously predicted by ab initio calculations accounting for magnetic interactions. We argue that in situ high-pressure experiments at cryogenic temperatures, which suppress time-dependent transformation triggered at higher temperatures, present a unique opportunity to significantly improve our understanding of these complex alloys. [ABSTRACT FROM AUTHOR]

Very recently, open-cage metallo-azafullerenes PbC100N4H4 and Pb2C100N4H4 containing one Pb–N4–C moiety have been synthesized via the electron beam. Herein, we utilized density functional theory calculations in combination with ab initio molecular dynamics (AIMD) simulations to study the geometric and electronic structures, bonding properties, thermodynamic stability, and catalytic performance of MC100N4H4 and M2C100N4H4 (M = Ge, Sn, Pb). Metal–nitrogen distances and metal–metal distances increase along with the metal radius while the metal atom is positively charged. Energy decomposition analysis revealed that the bonding interactions between M and the C100N4H4 fragment could be described as the donor–acceptor interaction between M(ns0(n−1)d10np4) and C100N4H4 fragment, in which the orbital interactions terms contribute more than the electrostatic interactions. AIMD simulations demonstrate that those metallo-azafullerenes exhibit thermodynamic stability at room temperature. These metallo-azafullerenes, which could serve as typical carbon-supported single-atom catalysts, possess enhanced catalytic performance toward the oxygen reduction reaction (ORR) compared to the planar catalysts, which is attributed to the curvature of metallo-azafullerenes. GeC100N4H4 and SnC100N4H4 exhibit high catalytic performance in the 4e-ORR pathway to H2O, whereas only PbC100N4H4 is suitable for the 2e-ORR reaction pathway because of the difficulty in obtaining electrons. All M2C100N4H4 favors the 4e-reaction pathway due to the presence of the axial metal atom. Our finding of open-cage metallo-azafullerenes as efficient single-atom catalysts holds profound implications for both fundamental research in catalysis and practical applications in fuel cells and other electrochemical devices. [ABSTRACT FROM AUTHOR]

Global diabatic potential energy surfaces (PESs) of CH2+ are constructed using the neural network method with a specific function based on 18 213 ab initio points. The multi-reference configuration interaction method with the aug-cc-pVQZ basis set is adopted to perform the ab initio calculations. The topographical properties of the diabatic PESs are examined in detail. In general, the diabatic PESs provide an accurate quasi-diabatic representation. To validate the diabatic PESs, the dynamics studies of the C+(2P1/2, 3/2) + H2 (v0 = 0, j0 = 0) → H + CH+(X1Σ+) reaction are performed using the time-dependent wave packet method. The reaction probabilities, integral cross sections, differential cross sections, and rate constants are calculated and compared with the experimental and theoretical results. Non-adiabatic dynamics results are in good agreement with experimental data. In addition, the non-adiabatic effect in the C+(2P1/2, 3/2) + H2 reaction is significant due to the non-adiabatic results being obviously larger than adiabatic values. The reasonable non-adiabatic dynamics results indicate that present diabatic PESs can be recommended for any type of dynamics study. [ABSTRACT FROM AUTHOR]

Global diabatic potential energy surfaces (PESs) of CH2+ are constructed using the neural network method with a specific function based on 18 213 ab initio points. The multi-reference configuration interaction method with the aug-cc-pVQZ basis set is adopted to perform the ab initio calculations. The topographical properties of the diabatic PESs are examined in detail. In general, the diabatic PESs provide an accurate quasi-diabatic representation. To validate the diabatic PESs, the dynamics studies of the C+(2P1/2, 3/2) + H2 (v0 = 0, j0 = 0) → H + CH+(X1Σ+) reaction are performed using the time-dependent wave packet method. The reaction probabilities, integral cross sections, differential cross sections, and rate constants are calculated and compared with the experimental and theoretical results. Non-adiabatic dynamics results are in good agreement with experimental data. In addition, the non-adiabatic effect in the C+(2P1/2, 3/2) + H2 reaction is significant due to the non-adiabatic results being obviously larger than adiabatic values. The reasonable non-adiabatic dynamics results indicate that present diabatic PESs can be recommended for any type of dynamics study. [ABSTRACT FROM AUTHOR]

Conducting a comprehensive molecular-level evaluation of a photoacid generator (PAG) and its subsequent impact on lithography performance can facilitate the rational design of a promising 193 nm photoresist tailored to specific requirements. In this study, we integrated spectroscopy and computational techniques to meticulously investigate the pivotal factors of three prototypical PAG anions, p-toluenesulfonate (pTS−), 2-(trifluoromethyl)benzene-1-sulfonate (TFMBS−), and triflate (TF−), in the lithography process. Our findings reveal a significant redshift in the absorption spectra caused by specific PAG anions, attributed to their involvement in electronic transition processes, thereby enhancing the transparency of the standard PAG cation, triphenylsulfonium (TPS+), particularly at ∼193 nm. Furthermore, the electronic stability of PAG anions can be enhanced by solvent effects with varying degrees of strength. We observed the lowest vertical detachment energy of 6.6 eV of pTS− in PGMEA solution based on the polarizable continuum model, which prevents anion loss at 193 nm lithography. In addition, our findings indicate gas-phase proton affinity values of 316.4 kcal/mol for pTS−, 308.1 kcal/mol for TFMBS−, and 303.2 kcal/mol for TF−, which suggest the increasing acidity strength, yet even the weakest acid pTS− is still stronger than strong acid HBr. The photolysis of TPS+-based PAG, TPS+·pTS−, generated an excited state leading to homolysis bond cleavage with the lowest reaction energy of 83 kcal/mol. Overall, the PAG anion pTS− displayed moderate acidity, possessed the lowest photolysis reaction energy, and demonstrated an appropriate redshift. These properties collectively render it a promising candidate for an effective acid producer. [ABSTRACT FROM AUTHOR]

We report the x-ray absorption spectrum (XAS) of the tert-butyl radical, C4H9. The radical was generated pyrolytically from azo-tert-butane, and the XAS of the pure radical was obtained by subtraction of spectra recorded at different temperatures. The bands in the XAS were assigned by ab initio calculations that are in very good agreement with the experimental data. The lowest energy signal in the XAS is assigned to the C1s electron transition from the central carbon atom to the singly occupied molecular orbital (SOMO), while higher transitions correspond to C1s excitations from terminal carbon atoms. Furthermore, we investigated the fragmentation of the radical following resonant C1s excitation by electron–ion-coincidence spectroscopy. Several fragmentation channels were identified. The C1s excitation of the terminal carbons is associated with a stronger fragmentation tendency compared to the lowest C1s excitation of the central carbon into the SOMO. For this core excited state, we still observe an intact parent ion, C 4 H 9 + , and a comparatively higher tendency to dissociate into CH 3 + + C 3 H 6 + . [ABSTRACT FROM AUTHOR]

In this work, we performed static density functional theory calculations and ab initio metadynamics simulations to systematically investigate the association mechanisms and dynamic structures of four kinds of ion pairs that could be formed before the nucleation of CaCO3. For Ca2+– H C O 3 − and Ca2+– C O 3 2 − pairs, the arrangement of ligands around Ca2+ evolves between the six-coordinated octahedral structure and the seven-coordinated pentagonal bipyramidal structure. The formation of ion pairs follows an associative ligand substitution mechanism. Compared with H C O 3 − , C O 3 2 − exhibits a stronger affinity to Ca2+, leading to the formation of a more stable precursor phase in the prenucleation stage, which promotes the subsequent CaCO3 nucleation. In alkaline environments, excessive OH− ions decrease the coordination preference of Ca2+. In this case, the formation of Ca(OH)+– C O 3 2 − and Ca(OH)2– C O 3 2 − pairs favors the dissociative ligand substitution mechanism. The inhibiting effects of OH− ion on the CaCO3 association can be interpreted from two aspects, i.e., (1) OH− neutralizes positive charges on Ca2+, decreases the electrostatic interactions between Ca2+ and C O 3 2 − , and thus hinders the formation of the CaCO3 monomer, and (2) OH− decreases the capacity of Ca2+ for accommodating O, making it easier to separate Ca2+ and C O 3 2 − ions. Our findings on the ion association behaviors in the initial stage of CaCO3 formation not only help scientists evaluate the impact of ocean acidification on biomineralization but also provide theoretical support for the discovery and development of more effective approaches to manage undesirable scaling issues. [ABSTRACT FROM AUTHOR]

Improved two-body and three-body potentials for helium have been used to calculate from first principles the third density and acoustic virial coefficients for both 4He and 3He. For the third density virial coefficient C(T), uncertainties have been reduced by a factor of 4–5 compared to the previous state of the art; the accuracy of first-principles C(T) now exceeds that of the best experiments by more than two orders of magnitude. The range of calculations has been extended to temperatures as low as 0.5 K. For the third acoustic virial coefficient γa(T), we applied the Schlessinger point method, which can calculate γa and its uncertainty based on the C(T) data, overcoming some limitations of direct path-integral calculation. The resulting γa are calculated at temperatures down to 0.5 K; they are consistent with available experimental data but have much smaller uncertainties. The first-principles data presented here will enable improvement of primary temperature and pressure metrology based on gas properties. [ABSTRACT FROM AUTHOR]

The x-ray restrained wave function (XRW) method is a quantum crystallographic technique to extract wave functions compatible with experimental x-ray diffraction data. The approach looks for wave functions that minimize the energies of the investigated systems and also reproduce sets of x-ray structure factors. Given the strict relationship between x-ray structure factors and electron distributions, the strategy practically allows determining wave functions that correspond to given (usually experimental) electron densities. In this work, the capabilities of the XRW approach were further tested. The aim was to evaluate whether the XRW technique could serve as a tool for suggesting new exchange-correlation functionals for density functional theory or refining existing ones. Additionally, the ability of the method to address the influences of the crystalline environment was also assessed. The outcomes of XRW computations were thus compared to those of traditional gas-phase, embedding quantum mechanics/molecular mechanics, and fully periodic calculations. The results revealed that, irrespective of the initial conditions, the XRW computations practically yield a consensus electron density, in contrast to the currently employed density functional approximations (DFAs), which tend to give a too large range of electron distributions. This is encouraging in view of exploiting the XRW technique to develop improved functionals. Conversely, the calculations also emphasized that the XRW method is limited in its ability to effectively address the influences of the crystalline environment. This underscores the need for a periodic XRW technique, which would allow further untangling the shortcomings of DFAs from those inherent to the XRW approach. [ABSTRACT FROM AUTHOR]

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 − 3 to 3 × 10 3 s − 1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s − 1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well. [ABSTRACT FROM AUTHOR]

Quantum Dot Cellular Automata (QCA) is an emerging trend in the field of nanoelectronics, and computing can be regarded as an alternative to the traditional complementary metal–oxide–semiconductor technology. The paper is devoted to the study of the key functional properties of the cells for molecular QCA based on mixed valence molecules. The theoretical results for the heat dissipation under the conditions of the fast nonadiabatic switching event and cell–cell response function are obtained in the framework of the quantum-mechanical vibronic approach. These results are parameterized using the previous reliable ab initio calculations performed for oxidized norbornadiene and its polycyclic derivatives with variable lengths of the bridge. The comparative analysis of the dimeric and bidimeric molecular cells composed of these compounds is given. It is underlined that the conditions of a strong non-linear response and a low heat release are contradictory. However, despite this problem, a parametric regime is proposed, which provides a low heat release in combination with a strong nonlinear response of the working cell to the electric field induced by the polarized driver cell. [ABSTRACT FROM AUTHOR]

Acetaldehyde is a very relevant atmospheric species whose photodissociation has been extensively studied in the first absorption band both experimentally and theoretically. Very few works have been reported on acetaldehyde photodissociation at higher excitation energies. In this work, the photodissociation dynamics of acetaldehyde is investigated by means of high-level multireference configuration interaction ab initio calculations. Five different fragmentation pathways of acetaldehyde are explored by calculating the potential-energy curves of the ground and several excited electronic states along the corresponding dissociating bond distances. The excitation energy range covered in the study is up to 10 eV, nearly the ionization energy of acetaldehyde. We intend to rationalize the available experimental results and, in particular, to elucidate why some of the studied fragmentation pathways are experimentally observed in the different excitation energy regions and some others are not. Based on the shape of the calculated potential curves, we are able to explain the main findings of the available experiments, also suggesting possible dynamical dissociation mechanisms in the different energy regions. Thus, the reported potential curves are envisioned as a useful tool to interpret the currently available experiments as well as future ones on acetaldehyde photodissociation at excitation wavelengths in the range studied here. [ABSTRACT FROM AUTHOR]

Ethane ( C 2 H 6 ) is anticipated to be the most stable compound within the carbon–hydrogen system under the 100 GPa pressure range. Nevertheless, the properties of ethane under pressure are still poorly documented. Here, we present a comprehensive study of the structural and vibrational properties of C 2 H 6 in a diamond anvil cell at pressures up to 150 GPa. To obtain detailed data, ethane single-crystal was grown in a helium pressure-transmitting medium. Utilizing single-crystal x-ray diffraction, the distortion mechanism between the tetragonal and monoclinic phases, occurring over the 3.2–5.2 GPa pressure range, is disclosed. Subsequently, no phase transition is observed up to 150 GPa. The accurately measured compression curve is compared to various computational approximations. The vibrational modes measured by Raman spectroscopy and infrared absorption are well identified, and their evolution is well reproduced by ab initio calculations. In particular, an unusual anticrossing phenomenon occurs near 40 GPa between a rocking and a stretching mode, likely attributable to intermolecular interactions through hydrogen bonding. [ABSTRACT FROM AUTHOR]

The isotope shifts in electron affinities of Pb were measured by Walter et al. [Phys. Rev. A 106, L010801 (2022)] to be −0.002(4) meV for 207–208Pb and −0.003(4) meV for 206–208Pb by scanning the threshold of the photodetachment channel Pb− ( S 3 / 2 ◦ 4 ) − Pb (3P0), while Chen and Ning reported 0.015(25) and −0.050(22) meV for the isotope shifts on the binding energies measured relative to 3P2 using the SEVI method [J. Chem. Phys. 145, 084303 (2016)]. Here we revisited these isotope shifts by using our second-generation SEVI spectrometer and obtained −0.001(15) meV for 207–208Pb and −0.001(14) meV for 206–208Pb, respectively. In order to aid the experiment by theory, we performed the first ab initio theoretical calculations of isotope shifts in electron affinities and binding energies of Pb, as well as the hyperfine structure of 207Pb−, by using the MCDHF and RCI methods. The isotope shifts in electron affinities of 207–208Pb and 206–208Pb are −0.0023(8) and −0.0037(13) meV for the 3P0 channel, respectively, in good agreement with Walter et al.'s measurements. The isotope shifts in binding energies relative to 3P1,2, −0.0015(8) and −0.0026(13) meV for 207–208Pb and 206–208Pb, respectively, are compatible with the present measurements. The hyperfine constant for the ground state of 207Pb− obtained by the present calculations, A ( S 3 / 2 ◦ 4 ) = − 1118 MHz, differs by a factor of 3 from the previous estimation by Bresteau et al. [J. Phys. B: At., Mol. Opt. Phys. 52, 065001 (2019)]. The reliability is supported by the good agreement between the theoretical and experimental hyperfine parameters of 209Bi. [ABSTRACT FROM AUTHOR]

The isotope shifts in electron affinities of Pb were measured by Walter et al. [Phys. Rev. A 106, L010801 (2022)] to be −0.002(4) meV for 207–208Pb and −0.003(4) meV for 206–208Pb by scanning the threshold of the photodetachment channel Pb− ( S 3 / 2 ◦ 4 ) − Pb (3P0), while Chen and Ning reported 0.015(25) and −0.050(22) meV for the isotope shifts on the binding energies measured relative to 3P2 using the SEVI method [J. Chem. Phys. 145, 084303 (2016)]. Here we revisited these isotope shifts by using our second-generation SEVI spectrometer and obtained −0.001(15) meV for 207–208Pb and −0.001(14) meV for 206–208Pb, respectively. In order to aid the experiment by theory, we performed the first ab initio theoretical calculations of isotope shifts in electron affinities and binding energies of Pb, as well as the hyperfine structure of 207Pb−, by using the MCDHF and RCI methods. The isotope shifts in electron affinities of 207–208Pb and 206–208Pb are −0.0023(8) and −0.0037(13) meV for the 3P0 channel, respectively, in good agreement with Walter et al.'s measurements. The isotope shifts in binding energies relative to 3P1,2, −0.0015(8) and −0.0026(13) meV for 207–208Pb and 206–208Pb, respectively, are compatible with the present measurements. The hyperfine constant for the ground state of 207Pb− obtained by the present calculations, A ( S 3 / 2 ◦ 4 ) = − 1118 MHz, differs by a factor of 3 from the previous estimation by Bresteau et al. [J. Phys. B: At., Mol. Opt. Phys. 52, 065001 (2019)]. The reliability is supported by the good agreement between the theoretical and experimental hyperfine parameters of 209Bi. [ABSTRACT FROM AUTHOR]

The in-plane thermal conductivity of hexagonal boron nitride (h-BN) with varying thicknesses is a key property that affects the performance of various applications from electronics to optoelectronics. However, the transition of the thermal conductivity from two-dimensional (2D) to three-dimensional (3D) h-BN remains elusive. To answer this question, we have developed a machine learning interatomic potential within the neuroevolution potential (NEP) framework for h-BN, achieving a high accuracy akin to ab initio calculations in predicting its thermal conductivity and phonon transport from monolayer to multilayers and bulk. Utilizing molecular dynamics simulations based on the NEP, we predict the thermal conductivity of h-BN with a thickness up to ∼100 nm, demonstrating that its thermal conductivity quickly decreases from the monolayer and saturates to the bulk value above four layers. The saturation of its thermal conductivity is attributed to the little change in phonon group velocity and lifetime as the thickness increases beyond four layers. In particular, the weak thickness dependence of phonon lifetime in h-BN with a nanoscale thickness results from its extremely high phonon focusing along the in-plane direction. This research bridges the knowledge gap of phonon transport between 2D and 3D h-BN and will benefit the thermal design and performance optimization of relevant applications. [ABSTRACT FROM AUTHOR]

This study was focused on the photochemistry of OAlOH and three possible pathways, which were studied with high-level multireference configuration interaction ab initio calculations. We computed cuts of the six-dimensional potential energy surfaces for the ground, the lowest singlet and triplet excited states, and probed the photodissociation mechanisms and the stabilities. The OAlOH electronic spectrum, with an energy reaching 7.15 eV, contained four prominent peaks. Photodissociation to AlO, OH, and AlOH constituted a plausible mechanism within the deep-UV range (λ = 250.4 nm). Our data indicated the photostability of OAlOH in the near-UV‒Vis region, so detection with laser-induced fluorescence is possible. Fluorescence and phosphorescence may occur upon excitation at 363.5 nm. The roles of OAlOH in the photochemical reactions of Al-bearing molecules in the upper atmosphere and VY Canis Majoris are discussed. [ABSTRACT FROM AUTHOR]

Face-to-face annealed and sputter-deposited aluminum nitride (FFA Sp-AlN) has potential in deep-ultraviolet light-emitting devices. Herein, the effects of the substrate off-cut angle (θsub) from an r-plane sapphire toward the c-axis projection direction and sputtering temperature (Tsp) on the crystallinity and surface morphology of a-plane AlN films are investigated. Increasing θsub in the minus-off direction, which occurs when the substrate surface approaches the sapphire c-plane, and lowering Tsp suppress the mixing of anomalous non-a-direction oriented domains. This reduced mixing enhances the surface flatness and crystallinity of a-plane FFA Sp-AlN. Moreover, the c-axis direction of the a-plane AlN film is inverted depending on the substrate off-cut angle. Ab initio calculations indicate that the interface stability between the r-plane sapphire substrate and the a-plane AlN film can explain the dependence of the crystallinity and c-axis orientation of FFA Sp-AlN on the surface off-cut of the r-plane sapphire. [ABSTRACT FROM AUTHOR]

Phosphate-based glasses (PBGs) in the CaO–Na2O–P2O5 system have diverse applications as biomaterials due to their unique dissolution properties. A reactive force field (ReaxFF) has been developed to simulate these materials at the atomic level. The ReaxFF parameters of PBGs, including the interaction between phosphorus and calcium atoms, have been developed using a published code based on genetic algorithms. The training data, including the atomic charges, atomic forces, bond and angle parameters, and different differential energies, are chosen and measured from static quantum-mechanical calculations and ab initio molecular dynamics of PBGs. We did different short- and medium-range structural analyses on the bulk simulated PBGs with different compositions to validate the developed potential. Radial and angular distribution functions and coordination numbers of network formers and modifiers, as well as the network connectivity of the glass, are in agreement with experimental and previous simulations using both shell-model classical force fields and ab initio simulated data; for example, the coordination number of phosphorus is 4.0. This successful development of ReaxFF parameters being able to describe the bulk PBGs enables us to work on the dissolution behavior of the glasses, including the interaction of water molecules with PBGs, in future works. [ABSTRACT FROM AUTHOR]

C H 3 + , a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by using the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction ab initio calculations are done, and full-dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full-dimensional quantum wave packet method in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids, allowing us to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion bar, result in a lesser destruction of C H 3 + compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry. [ABSTRACT FROM AUTHOR]

An accurate description of the long-range (LR) interaction is essential for understanding the collision between cold or ultracold molecules. However, to our best knowledge, there lacks a general approach to construct the intermolecular potential energy surface (IPES) between two arbitrary molecules and/or atoms in the LR region. In this work, we derived analytical expressions of the LR interaction energy, using the multipole expansion of the electrostatic interaction Hamiltonian and the non-degenerate perturbation theory. To make these formulae practical, we also derived the independent Cartesian components of the electrostatic properties, including the multipole moments and polarizabilities, of the monomer for a given symmetry using the properties of these components and the group-theoretical methods. Based on these newly derived formulae, we developed a FORTRAN program, namely ABLRI, which is capable of calculating the interaction energy between two arbitrary monomers both in their non-degenerate electronic ground states at large separations. To test the reliability of this newly developed program, we constructed IPESs for the electronic ground state of H2O–H2 and O2–H systems in the LR region. The interaction energy computed by our program agreed well with the ab initio calculation, which shows the validity of this program. [ABSTRACT FROM AUTHOR]

Thorium nitrides have been the topic of intense studies due to their prospective applications as advanced nuclear fuels. The phase diagram of the Th–N scheme, however, continues unknown at low temperatures and extremely high pressures. In this article, we examine the Th–N system's phase diagram up to 300 GPa from the first-principle approach using universal structure predictor: evolutionary Xtallography (USPEX) method. Apart from the experimentally observed phase (ThN, Th2N3, and Th3N4), there are several unique chemical stoichiometries, i.e., ThN3, ThN4, ThN6, ThN8, ThN10, and ThN12 are found to have stability fields on the Th–N phase diagram at pressure of 3.0, 32, 100, 42, 28, and 236 GPa along with previously predicted composition ThN2 at 3.5 GPa. The structural stability of the predicted compositions is further assessed by evaluating the elastic and dynamic stability. Out of all above mentioned compositions, ThN3 is possibly a metastable one at 0 GPa. Electronic structure calculations predict that all newly discovered compositions are metallic except ThN10, which is semi-metallic at high pressures. Further, we predict that ThN4 and ThN6 have high electron–phonon coupling constant of 1.874 and 0.894 with Tc around 21.22 and 25.02 K, respectively, at 100 GPa. [ABSTRACT FROM AUTHOR]

Identifying the oxidation structure of two-dimensional interfaces is crucial to improve surface chemistry and electronic properties. Beyond graphene with only phenyl rings, a novel carbon-nitrogen material, C3N, presents an intrinsic heterogeneous surface morphology where each phenyl ring is encircled by six nitrogen atoms, yet its atomistic oxidation structure remains unclear. Here, combining a series of density functional theory calculations and ab initio molecular dynamics simulations, we demonstrate that thermodynamically favorable oxidation loci are confined to the phenyl ring, and kinetic transformations of oxidation structures are feasible along the phenyl ring, whereas those toward nitrogen atoms are proven to be extremely difficult. These results are attributed to the lower barrier of oxygen atom migration along the phenyl ring, while the significantly high barriers toward nitrogen atoms are due to the heterogeneous potential energy surface for oxygen–C3N interaction. This work highlights the significance of surface morphology on the characteristics of oxidation structure, offering insights into tunable electronic properties via confined interfacial oxidation. [ABSTRACT FROM AUTHOR]