Search

We performed high-level ab initio quantum chemical calculations, incorporating higher-order excitations, spin–orbit coupling (SOC), and the Gaunt interaction, to calculate the electron affinities (EAs) of alkaline earth (AE) metal atoms (Ca, Sr, Ba, and Ra), which are notably small. The coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method is insufficient to accurately calculate the EAs of AE metal atoms. Higher-order excitations proved crucial, with the coupled-cluster singles, doubles, and triples with perturbative quadruples [CCSDT(2)Q] method effectively capturing dynamic electron correlation effects. The contributions of SOC (ΔESOs) to the EAs calculated using the multireference configuration interaction method with the Davidson correction, including SOC, positively enhance the EAs; however, these contributions are overestimated. The Dirac–Hartree–Fock (DHF)-CCSD(T) method addresses this overestimation and provides reasonable values for ΔESO (ΔESO−D). Employing additional sets of diffuse and core–valence correlation basis sets is critical for accurately calculating the EAs of AE metal atoms. The contributions of the Gaunt interaction (ΔEGaunt) to the EAs of AE metal atoms are negligible. Notably, the CCSDT(2)Q with the complete basis set limit + ΔESO−D + ΔEGaunt produced EA values for Ca, Sr, and Ba that closely aligned with experimental data and achieved accuracy exceeding the chemical accuracy. Based on our findings, the accurately proposed EA for Ra is 9.88 kJ/mol. [ABSTRACT FROM AUTHOR]

The authors implanted P atoms at 50 and 140-keV energies with respective fluences of 1 × 10 14 and 2 × 10 14 /cm 2 into diamond thin films synthesized by chemical vapor deposition on the type-Ib diamond substrate formed by high-pressure and high-temperature conditions. The occupational sites of implanted P atoms were determined in each processing stage of P implantation, ion-beam-induced epitaxial crystallization (IBIEC) annealing by 3-MeV-Ne 2 + ion irradiation at 750 ° C, and thermal annealing at 850 ° C in vacuum, by quantitatively comparing the random and channeling yields in Rutherford backscattering (RBS) measurements. In the analysis of RBS spectra, simulated distributions were fitted to the experimentally obtained spectra and we quantitatively identified the occupational sites and ratio of implanted P atoms. In addition, we investigated the diffusion phenomenon of implanted P atoms during annealing processes from the depth profile of scattered He ions. Consequently, the averaged occupational ratio in the substitutional lattice site was stably recorded around 50%. From the depth profile in the RBS spectra, thermal diffusion of implanted P atoms was clearly observed during the last thermal annealing after MeV-IBIEC annealing. We finally discuss the potentials of the MeV-IBIEC irradiation at relatively low temperatures applying to useful and credible annealing methods for electrical activation as well as recovery of damaged crystallinity. [ABSTRACT FROM AUTHOR]

Deposition of Mo S 2 on Au(111) alters the electronic properties of Mo S 2 . In this study, we investigate the free-standing MoS 2 monolayer and the MoS 2 /Au(111) heterostructure, with and without strain, as well as defects of interest in memristive and neuromorphic applications. We focus on the so-called atomristor devices based on monolayer materials that achieve resistive switching characteristics with the adsorption and desorption of metal adatoms. Our study confirms that the formation of midgap states is the primary mechanism behind the resistive switching. Our results show that strain lowers the adsorption/desorption energies of Au+defect structures of interest, leading to more favorable switching energies, but simultaneously reduces the switching ratio between states of differing conductivities. The presence of the Au(111) substrate additionally introduces non-uniform amounts of strain and charge transfer to the MoS 2 monolayer. We propose that the induced strain contributes to the experimentally observed n - to p -type transition and Ohmic to Schottky transition in the MoS 2 monolayer. The charge transfer leads to a permanent polarization at the interface, which can be tuned by strain. Our study has important implications on the role of the electrode as being a source of the observed variability in memristive devices and as an additional degree of freedom for tuning the switching characteristics of the memristor device. [ABSTRACT FROM AUTHOR]

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ∼4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. These results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions. [ABSTRACT FROM AUTHOR]

Alchemical transformations showed that perturbation theory can be applied also to changes in the atomic nuclear charges of a molecule. The alchemical path that connects two different chemical species involves the conceptualization of a non-physical system in which an atom possess a non-integer nuclear charge. A correct quantum mechanical treatment of these systems is limited by the fact that finite size atomic basis sets do not define exponents and contraction coefficients for fractional charge atoms. This paper proposes a solution to this problem and shows that a smooth interpolation of the atomic orbital coefficients and exponents across the periodic table is a convenient way to produce accurate alchemical predictions, even using small size basis sets. [ABSTRACT FROM AUTHOR]

The electronic structures of three highly mismatched alloys (HMAs)—GeC(Sn), Ga(In)NAs, and BGa(In)As—were studied using density functional theory with HSE06 hybrid functionals, with an emphasis on the local environment near the mismatched, highly electronegative atom (B, C, and N). These alloys are known for their counterintuitive reduction in the bandgap when adding the smaller atom, due to a band anticrossing (BAC) or splitting of the conduction band. Surprisingly, the existence of band splitting was found to be completely unrelated to the local displacement of the lattice ions near the mismatched atom. Furthermore, in BGaAs, the reduction in the bandgap due to BAC was weaker than the increase due to the lattice constant, which has not been observed among other HMAs but may explain differences among experimental reports. While local distortion in GeC and GaNAs was not the cause for BAC, it was found to enhance the bandgap reduction due to BAC. This work also found that mere contrast in electronegativity between neighboring atoms does not induce BAC. In fact, surrounding the electronegative atom with elements of even smaller electronegativity than the host (e.g., Sn or In) consistently decreased or even eliminated BAC. For a fixed composition, moving Sn toward C and In toward either N or B was always energetically favorable and increased the bandgap, consistent with experimental annealing results. Such rearrangement also delocalized the conduction band wavefunctions near the mismatched atom to resemble the original host states in unperturbed Ge or GaAs, causing the BAC to progressively weaken. These collective results were consistent whether the mismatched atom was a cation (N), anion (B), or fully covalent (C), varying only with the magnitude of its electronegativity, with B having the least effect. The effects can be explained by charge screening of the mismatched atom's deep electrostatic potential. Together, these results help explain differences in the bandgap and other properties reported for HMAs from different groups and provide insight into the creation of materials with designer properties. [ABSTRACT FROM AUTHOR]

This Tutorial focuses on the physics of layered compounds under high pressure. We have chosen h-BN and III–VI layered materials as representative materials. h-BN layers are strictly two-dimensional. Layers in III–VI compounds are more complex, and subtle details in their structural behavior play an important role in the evolution of high pressure properties. They are also interesting because they contain a different number of layers in their primitive unit cell and/or have a different ionic character. We begin describing the structural evolution. We discuss the experimental challenges encountered as well as the main findings related to intra- and interlayer compressibility, polytype influence, and geometrical modifications induced by pressure inside the layers. We then describe lattice vibrations. The origin of the modes is reviewed, paying attention to the relationships between atom motions in different layers. We discuss the convenience of redefining the Grüneisen parameter and describe the behavior of rigid layer modes, soft modes, and Davidov pairs. The last section is devoted to the electronic properties. We show that the changes observed when passing from a single layer to a three-dimensional BN are qualitatively similar to those induced by high pressure. The pressure behavior of electronic transitions in III–VI layered compounds is very rich, revealing the subtle balance between intra- and inter-layer interactions. Finally, we take advantage of high pressure studies to explain the formation of the Mexican hat type of valence band at ambient conditions in single layers of InSe and GaSe, but not in three-dimensional compounds. [ABSTRACT FROM AUTHOR]

A nitrogen K-edge x-ray absorption near-edge structure (XANES) survey is presented for tetrapyrido[3,2-a:2′,3′-c:3″,2″-h:2‴,3‴-j]phenazine (tpphz)-bridged bimetallic assemblies that couple chromophore and catalyst transition metal complexes for light driven catalysis, as well as their individual molecular constituents. We demonstrate the high N site sensitivity of the N pre-edge XANES features, which are energetically well-separated for the phenazine bridge N atoms and for the individual metal-bound N atoms of the inner coordination sphere ligands. By comparison with the time-dependent density functional theory calculated spectra, we determine the origins of these distinguishable spectral features. We find that metal coordination generates large shifts toward higher energy for the metal-bound N atoms, with increasing shift for 3d < 4d < 5d metal bonding. This is attributed to increasing ligand-to-metal σ donation that increases the effective charge of the bound N atoms and stabilizes the N 1s core electrons. In contrast, the phenazine bridge N pre-edge peak is found at a lower energy due to stabilization of the low energy electron accepting orbital localized on the phenazine motif. While no sensitivity to ground state electronic coupling between the individual molecular subunits was observed, the spectra are sensitive to structural distortions of the tpphz bridge. These results demonstrate N K-edge XANES as a local probe of electronic structure in large bridging ligand motifs, able to distinctly investigate the ligand-centered orbitals involved in metal-to-ligand and ligand-to-ligand electron transfer following light absorption. [ABSTRACT FROM AUTHOR]

The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga–Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn–Sham limit is improved relative to other basis sets where convergence to the Kohn–Sham limit is typically not observed. [ABSTRACT FROM AUTHOR]

We describe an apparatus where many species of ultracold atoms can be simultaneously trapped and overlapped with many species of ions in a Paul trap. Several design innovations are made to increase the versatility of the apparatus while keeping the size and cost reasonable. We demonstrate the operation of a three-dimensional (3D) magneto-optical trap (MOT) of 7Li using a single external cavity diode laser. The 7Li MOT is loaded from an atomic beam, with atoms slowed using a Zeeman slower designed to work simultaneously for Li and Sr. The operation of a 3D MOT of 133Cs, loaded from a 2D MOT, is demonstrated, and provisions for MOTs of Rb and K in the same vacuum manifold exist. We demonstrate the trapping of 7Li+ and 133Cs+ at different settings of the Paul trap and their detection using an integrated time-of-flight mass spectrometer. We present results on low energy neutral–neutral collisions (133Cs–133Cs, 7Li–7Li, and 133Cs–7Li collisions) and charge–neutral collisions (133Cs+–133Cs and 7Li+–7Li collisions). We show evidence of sympathetic cooling of 7Li+ (133Cs+) due to collisions with the ultracold 7Li (133Cs). [ABSTRACT FROM AUTHOR]

The Effective Relativistic Coupling by Asymptotic Representation (ERCAR) approach is a method to generate fully coupled diabatic potential energy surfaces (PESs) including relativistic effects, especially spin–orbit coupling. The spin–orbit coupling of a full molecule is determined only by the atomic states of selected relativistically treated atoms. The full molecular coupling effect is obtained by a diabatization with respect to asymptotic states, resulting in the correct geometry dependence of the spin–orbit effect. The ERCAR approach has been developed over the last decade and initially only for molecules with a single relativistic atom. This work presents its extension to molecules with more than a single relativistic atom using the iodine molecule as a proof-of-principle example. The theory for the general multiple atomic ERCAR approach is given. In this case, the diabatic basis is defined at the asymptote where all relativistic atoms are separated from the remaining molecular fragment. The effective spin–orbit operator is then a sum of spin–orbit operators acting on isolated relativistic atoms. PESs for the iodine molecule are developed within the new approach and it is shown that the resulting fine structure states are in good agreement with spin–orbit ab initio calculations. [ABSTRACT FROM AUTHOR]

For the most precise incorporation of single impurities in silicon, which is utilized to create quantum devices, a monolayer of adatoms on the Si(100) surface and a dopant-containing molecule are used. Here, we studied the interaction of phosphorus tribromide with a chlorine monolayer with mono- and bivacancies using a scanning tunneling microscope (STM) at 77 K. The combination of different halogens in the molecule and the adsorbate layer enabled unambiguous identification of the structures after PBr3 dissociation on Si(100)-Cl. A Cl monolayer was exposed to PBr3 in the STM chamber, which allows us to compare the same surface areas before and after PBr3 adsorption. As a result of this comparison, we detected small changes in the chlorine layer and unraveled the molecular fragments filling mono- and bivacancies. Using density functional theory, we found that the phosphorus atom occupies a bridge position after dissociation of the PBr3 molecule, which primarily bonds with silicon in Cl bivacancies. These findings provide insight into the interaction of a dopant-containing molecule with an adsorbate monolayer on Si(100) and can be applied to improve the process of single impurity incorporation into silicon. [ABSTRACT FROM AUTHOR]

The photophysics of thiobases—nucleobases in which one or more oxygen atoms are replaced with sulfur atoms— vary greatly depending on the location of sulfonation. Not only are direct dynamics of a neutral thiobase impacted, but also the dynamics of excess electron accommodation. In this work, time-resolved photoelectron spectroscopy is used to measure binary anionic clusters of iodide and 4-thiouracil, I− · 4TU. We investigate charge transfer dynamics driven by excitation at 3.88 eV, corresponding to the lowest ππ* transition of the thiouracil, and at 4.16 eV, near the cluster vertical detachment energy. The photoexcited state dynamics are probed by photodetachment with 1.55 and 3.14 eV pulses. Excitation at 3.88 eV leads to a signal from a valence-bound ion only, indicating a charge accommodation mechanism that does not involve a dipole-bound anion as an intermediate. Excitation at 4.16 eV rapidly gives rise to dipole-bound and valence-bound ion signals, with a second rise in the valence-bound signal corresponding to the decay of the dipole-bound signal. The dynamics associated with the low energy ππ* excitation of 4-thiouracil provide a clear experimental proof for the importance of localized excitation and electron backfilling in halide–nucleobase clusters. [ABSTRACT FROM AUTHOR]

Four-center two-electron Coulomb integrals routinely appear in electronic structure algorithms. The resolution-of-the-identity (RI) is a popular technique to reduce the computational cost for the numerical evaluation of these integrals in localized basis-sets codes. Recently, Duchemin and Blase proposed a separable RI scheme [J. Chem. Phys. 150, 174120 (2019)], which preserves the accuracy of the standard global RI method with the Coulomb metric and permits the formulation of cubic-scaling random phase approximation (RPA) and GW approaches. Here, we present the implementation of a separable RI scheme within an all-electron numeric atom-centered orbital framework. We present comprehensive benchmark results using the Thiel and the GW100 test set. Our benchmarks include atomization energies from Hartree–Fock, second-order Møller–Plesset (MP2), coupled-cluster singles and doubles, RPA, and renormalized second-order perturbation theory, as well as quasiparticle energies from GW. We found that the separable RI approach reproduces RI-free HF calculations within 9 meV and MP2 calculations within 1 meV. We have confirmed that the separable RI error is independent of the system size by including disordered carbon clusters up to 116 atoms in our benchmarks. [ABSTRACT FROM AUTHOR]

Bottom-up methods for coarse-grained (CG) molecular modeling are critically needed to establish rigorous links between atomistic reference data and reduced molecular representations. For a target molecule, the ideal reduced CG representation is a function of both the conformational ensemble of the system and the target physical observable(s) to be reproduced at the CG resolution. However, there is an absence of algorithms for selecting CG representations of molecules from which complex properties, including molecular electronic structure, can be accurately modeled. We introduce continuously gated message passing (CGMP), a graph neural network (GNN) method for atomically decomposing molecular electronic structure sampled over conformational ensembles. CGMP integrates 3D-invariant GNNs and a novel gated message passing system to continuously reduce the atomic degrees of freedom accessible for electronic predictions, resulting in a one-shot importance ranking of atoms contributing to a target molecular property. Moreover, CGMP provides the first approach by which to quantify the degeneracy of "good" CG representations conditioned on specific prediction targets, facilitating the development of more transferable CG representations. We further show how CGMP can be used to highlight multiatom correlations, illuminating a path to developing CG electronic Hamiltonians in terms of interpretable collective variables for arbitrarily complex molecules. [ABSTRACT FROM AUTHOR]

In density-functional theory, the exchange–correlation (XC) energy can be defined exactly through the coupling-constant (λ) averaged XC hole n ̄ xc (r , r ′ ) , representing the probability depletion of finding an electron at r′ due to an electron at r. Accurate knowledge of n ̄ xc (r , r ′ ) has been crucial for developing XC energy density-functional approximations and understanding their performance for molecules and materials. However, there are very few systems for which accurate XC holes have been calculated since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of λ while the electron density remains fixed at the physical (λ = 1) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron–electron cusp using finite basis sets. Focusing on Hooke's atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the λ-integrated XC hole. Our results indicate that, for Hooke's atoms, the error introduced by the description of the electron–electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle-, and coupling-constant-averaged XC holes are also calculated and provide a benchmark against which the Perdew–Burke–Ernzerhof and local density approximation XC hole models are assessed. [ABSTRACT FROM AUTHOR]

Although it plays a critical role in the photophysics and catalysis of lanthanides, spin–orbit coupling of electrons on individual lanthanide atoms in small clusters is not well understood. The major objective of this work is to probe such coupling of the praseodymium (Pr) 4f and 6s electrons in Pr2O2 and Pr2O2+. The approach combines mass-analyzed threshold ionization spectroscopy and spin–orbit multiconfiguration second-order quasi-degenerate perturbation theory. The energies of six ionization transitions are precisely measured; the adiabatic ionization energy of the neutral cluster is 38 045 (5) cm−1. Most of the electronic states involved in these transitions are identified as spin–orbit coupled states consisting of two or more electron spins. The electron configurations of these states are 4f46s2 for the neutral cluster and 4f46s for the singly charged cation, both in planar rhombus-type structures. The spin–orbit splitting due to the coupling of the electrons on the separate Pr atoms is on the order of hundreds of wavenumbers. [ABSTRACT FROM AUTHOR]

We numerically study the strong-interaction limit of the exchange–correlation functional for neutral atoms and Bohr atoms as the number of electrons increases. Using a compact representation, we analyze the second-order gradient expansion, comparing it with the one for exchange (weak interaction limit). The two gradient expansions, at strong and weak interaction, turn out to be very similar in magnitude but with opposite signs. We find that the point-charge plus continuum model is surprisingly accurate for the gradient expansion coefficient at strong coupling, while generalized gradient approximations, such as Perdew–Burke–Ernzerhof (PBE) and PBEsol, severely underestimate it. We then use our results to analyze the Lieb–Oxford bound from the point of view of slowly varying densities, clarifying some aspects on the bound at a fixed number of electrons. [ABSTRACT FROM AUTHOR]

In this paper, we explore the molecular basis of combining photodynamic therapy (PDT), a light-triggered targeted anticancer therapy, with the traditional chemotherapeutic properties of the well-known cytotoxic agent gemcitabine. A photosensitizer prerequisite is significant absorption of biocompatible light in the visible/near IR range, ideally between 600 and 1000 nm. We use highly accurate multiconfigurational CASSCF/MS-CASPT2/MM and TD-DFT methodologies to determine the absorption properties of a series of gemcitabine derivatives with the goal of red-shifting the UV absorption band toward the visible region and facilitating triplet state population. The choice of the substitutions and, thus, the rational design is based on important biochemical criteria and on derivatives whose synthesis is reported in the literature. The modifications tackled in this paper consist of: (i) substitution of the oxygen atom at O2 position with heavier atoms (O → S and O → Se) to red shift the absorption band and increase the spin–orbit coupling, (ii) addition of a lipophilic chain at the N7 position to enhance transport into cancer cells and slow down gemcitabine metabolism, and (iii) attachment of aromatic systems at C5 position to enhance red shift further. Results indicate that the combination of these three chemical modifications markedly shifts the absorption spectrum toward the 500 nm region and beyond and drastically increases spin–orbit coupling values, two key PDT requirements. The obtained theoretical predictions encourage biological studies to further develop this anticancer approach. [ABSTRACT FROM AUTHOR]

A much-needed solution for the efficient modeling of strong coupling between matter and optical cavity modes is offered by mean-field mixed quantum–classical dynamics, where a classical cavity field interacts self-consistently with quantum states of matter through Ehrenfest's theorem. We previously introduced a modified mean-field approach, referred to as decoupled mean-field (DC-MF) dynamics, wherein vacuum fluctuations of the cavity field are decoupled from the quantum-mechanical ground state as a means to resolve an unphysical drawing of energy from the vacuum fluctuations by a two-level atom. Here, we generalize DC-MF dynamics for an arbitrary number of (nondegenerate) atomic levels and show that it resolves an unphysical lack of emission from a three-level atom predicted by conventional mean-field dynamics. We furthermore show DC-MF to provide an improved description of reabsorption and (resonant) two-photon emission processes. [ABSTRACT FROM AUTHOR]

Preparation of atomic Bose–Einstein condensate (BEC) with tunable interactions in a quasi-one-dimensional (quasi-1D) optical trap is essential for both the observation of bright matter-wave solitons and the quantum simulation based on the discrete atomic momentum states. However, the quasi-1D BEC has been obtained by a complex process, which includes the creation of a three-dimensional BEC and its adiabatic transform into a quasi-1D trap. Here, we report the direct creation of a quasi-1D BEC of 133 Cs atoms by the fast evaporative cooling of ultracold atoms prepared by the degenerated Raman sideband cooling. We produce the pure BEC of up to 5.5 × 10 4 atoms in a quasi-1D optical trap with an evaporative time of 6 s. We demonstrate the anisotropic expansion of the atomic cloud after the release from the quasi-1D trap and study the dependence of both the phase space density and the temperature on the number of atoms in the trap during the evaporative cooling. Our study facilitates the promising applications of quasi-1D interacting atomic gases. [ABSTRACT FROM AUTHOR]

Periodic density functional theory was employed to investigate the impact of chemical modifications on the properties of π-conjugated acenedithiophene molecular crystals. Here, we highlight the importance of the β-methylthionation effect, the position of the sulfur atoms of the thiacycle group and their size, and the number of central benzene rings in the chemical modification strategy. Our results show that the introduction of the methylthio groups at the β-positions of the thiophene and the additional benzene ring at the center of the BDT crystal structure are a promising strategy to improve the performance of organic semiconductors, as observed experimentally. We found that β-MT-ADT exhibits large charge carrier mobility, which is in good agreement with the experimental results and comparable to that of rubrene. In addition, the electronic and optical properties of these ambipolar materials suggest promising performances with β-MT-ADT > ADT > β -MT-NDT > NDT > BEDT-BDT > β -MT-BDT > BDT. Moreover, functionalization with thiacycle-fused sulfur atoms of different sizes and numbers improve the properties of BDT but is still less efficient than the methylthionation effect. Overall, our findings suggest a promising molecular modification strategy for possibly high performance ambipolar organic semiconducting materials. [ABSTRACT FROM AUTHOR]

Inspired by the single-bonded nitrogen chains stabilized by tetravalent cerium, pentavalent tantalum, and hexavalent tungsten atoms, we explored the possibility of single-bonded nitrogen polymorphs stabilized by trivalent lanthanum ions. To achieve this, we utilized the crystal structure search method on the phase diagram of binary La–N compounds. We identified three novel thermodynamically stable phases, the C2/c LaN3, P-1 LaN4, and P-1 LaN8. Among them, the C2/c phase with infinite helical poly-N6 chains becomes thermodynamically stable above 50 GPa. Each nitrogen atom in the poly-N6 chain acquires one extra electron, and the spiral chain is purely single-bonded. The C2/c phase has an indirect band gap of ∼1.6 eV at 60 GPa. Notably, the band gap exhibits non-monotonic behavior, decreases first and then increases with increasing pressure. This abnormal behavior is attributed to the significant bonding of two La–N bonds at around 35 GPa. Phonon spectrum calculations and AIMD simulations have confirmed that the C2/c phase can be quenched to ambient conditions with slight distortion, and it exhibits excellent detonation properties. Additionally, we also discovered armchair-like nitrogen chains in LaN4 and the armchair and zigzag-like mixed nitrogen chains in LaN8. These results provide valuable insights into the electronic and bonding properties of nitrides under high pressure and may have important implications for the design and development of novel functional materials. [ABSTRACT FROM AUTHOR]

The H atom tunneling dissociation dynamics of the S1 state of meta- or para-cresol has been investigated by using the picosecond time-resolved pump-probe spectroscopy in a state-specific manner. The S1 state lifetime (mainly due to the H atom tunneling reaction) is found to be mode-dependent whereas it quickly converges and remains constant as the rapid intramolecular vibrational energy redistribution (IVR) starts to participate in the S1 state relaxation with the increase of the S1 internal energy (Eint). The IVR rate and its change with increasing Eint have been reflected in the parent ion transients taken by tuning the total energy (hνpump + hνprobe) just above the adiabatic ionization threshold (so that the dissipation of the initial mode-character could be monitored as a function of the reaction time), indicating that the mode randomization rate into the S1 isoenergetic manifolds exceeds the tunneling rate quite early in terms of Eint for m-cresol (≤∼1200 cm−1) or p-cresol (≤∼800 cm−1) compared to the case of phenol (≤∼1800 cm−1). Though the H atom tunneling dynamics of phenol (S1) seems to be little influenced by the methyl substitution on the either m- or p-position, the IVR rate has been found to be strongly accelerated due to the sharply-increasing (S1) density of states with increasing Eint due to the pivotal role of the low-frequency CH3 torsional mode. [ABSTRACT FROM AUTHOR]

Constructing metal-semiconductor interfaces by loading metal atoms onto two-dimensional material to build atomically dispersed single-atom catalysts (SACs) has emerged as a new frontier for improving atom utilization and designing multifunctional electrocatalysts. Nowadays, studies on black phosphorus nanosheets in electrocatalysis have received much attention and the successful preparation of metal nanoparticle/black phosphorus (BP) hybrid electrocatalysts indicates BP nanosheets can serve as a potential support platform for SACs. Herein, by using large-scale ab initio calculations, we explored a large composition space of SACs with transition metal atoms supported on BP monolayer (M-BP) and built a comprehensive picture of activity trend, stability, and electronic origin towards oxygen reduction and evolution reaction (ORR and OER) and hydrogen evolution reaction (HER). The results show that the catalytic activity can be widely tuned by reasonable regulation of metal atoms. Ni-, Pd-, and Pt-BP could effectively balance the binding strength of the target intermediates, thus achieving efficient bifunctional activity for OER and ORR. Favorable bifunctional catalytic performance for OER and HER can be realized on Rh-BP. Especially, Pt-BP exhibits promising trifunctional activity towards OER, ORR, and HER. Multiple-level corrections among overpotential, Gibbs free energy, orbital population, and d-band center reveal that the trend and origin of catalytic activity are intrinsically determined by the d-band center of metal sites. The thermodynamic and dynamic stability simulations demonstrate that the active metal centers are firmly anchored on BP substrate with intact M-P bonds. These findings provide a theoretical basis for the rational design of BP-based SACs toward promising multifunctional activity. [ABSTRACT FROM AUTHOR]

Many low-cost or semiempirical quantum mechanical-based electronic structure methods suffer from the use of unpolarized minimal atomic orbital (AO) basis sets. In this work, we overcome this limitation by a fully DFT variationally optimized, adaptive minimal basis set consistently available for the elements up to radon (Z = 86). The new key feature is to make the linear coefficients of the primitive Gaussians in a contracted AO dependent on the effective atomic charge of the atom in the molecule, i.e., each symmetry-unique atom obtains its "own" specifically adapted basis functions. In this way, the physically important "breathing" of the AOs in a molecule with (a) atomic charge (expansion/contraction for anionic/cationic states) and (b) the number of close-lying bonded neighbor atoms is accounted for. The required atomic charges are obtained from a specially developed extended Hückel type Hamiltonian and the coordination numbers from the molecule geometry. Proper analytical derivatives of the resulting adaptive basis functions can easily be derived. Moreover, the basis functions are electric field-dependent, thus improving the description of, e.g., dipole moments and polarizabilities. The new basis set termed q-vSZP (charge dependent valence single-ζ, polarized) is thoroughly benchmarked for atomic/molecular and thermochemical properties compared to standard minimal and double-ζ basis sets at the DFT level with the accurate ωB97X-D4 functional. It is shown that q-vSZP is clearly superior to existing minimal basis sets, often reaching double-ζ quality or even better results. We expect it to be the optimal choice in future semiempirical quantum mechanical methods. [ABSTRACT FROM AUTHOR]

We introduce ACEpotentials.jl, a Julia-language software package that constructs interatomic potentials from quantum mechanical reference data using the Atomic Cluster Expansion [R. Drautz, Phys. Rev. B 99, 014104 (2019)]. As the latter provides a complete description of atomic environments, including invariance to overall translation and rotation as well as permutation of like atoms, the resulting potentials are systematically improvable and data efficient. Furthermore, the descriptor's expressiveness enables use of a linear model, facilitating rapid evaluation and straightforward application of Bayesian techniques for active learning. We summarize the capabilities of ACEpotentials.jl and demonstrate its strengths (simplicity, interpretability, robustness, performance) on a selection of prototypical atomistic modelling workflows. [ABSTRACT FROM AUTHOR]

The mechanisms of carbon sticking reactions to C36 and C–C80 fullerenes were investigated with molecular dynamics simulations (MD) using the Second-generation Reactive Empirical Bond Order (SREBO) and Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potentials that were specifically optimized for carbon-carbon interactions. Results showed the existence of three possible sticking configurations where the projectile atom can stick either to one, two or three atoms of the target fullerene. They also showed that although the two potentials give similar magnitudes for the sticking cross-sections, they yield fairly different results as far as sticking mechanisms and configurations at thermal collision-energies, i.e., in the range 0.05–0.5 eV, are concerned. While AIREBO, that takes into account the long-range Lennard-Jones interaction, essentially results in a surface-sticking configuration with a single atom of the target fullerene, SREBO potential yields both surface- and two neighbors-sticking (2N-sticking) configurations. The fullerene structure is preserved in the last configuration while it can be recovered by a 2000 K annealing in the former configuration. Results obtained with SREBO eventually showed larger sticking probabilities for C36 as compared with C80. In spite of this, the sticking cross-sections obtained for C80 are similar to or even larger than those obtained for C36 due to the larger size of C80 that compensates for its smaller sticking probabilities. [ABSTRACT FROM AUTHOR]

The quantum state-to-state rotationally inelastic quenching of N2O by colliding with a He atom is studied on an ab initio potential energy surface with N2O lying on its vibrational ground state. The cross sections for collision energies from 10−6–100 cm−1 and rate constants from 10−5–10 K are calculated employing the fully converged quantum close-coupling method for the quenching of the j = 1–6 rotational states of N2O. Numerous van der Waals shapes or Feshbach resonances are observed; the cross sections of different channels are found to follow the Wigner scaling law in the cold threshold regime and may intersect with each other. In order to interpret the mechanism and estimate the cross sections of the rotational energy transfer, we propose a minimal classical model of collision between an asymmetric double-shell ellipsoid and a point particle. The classical model reproduces the quantum scattering results and points out the attractive interactions and the potential asymmetry can affect the collision process. The resulting insights are expected to expand our interpretations of inelastic scattering and energy transfer in molecular collisions. [ABSTRACT FROM AUTHOR]

A new approach for obtaining interacting quantum atoms-defined components of binding energy of intermolecular interactions, which bypasses the use of standard six-dimensional integrals and two-particle reduced density matrix (2-RDM) reconstruction, is proposed. To examine this approach, three datasets calculated within the density functional theory framework using the def2-TZVP basis have been explored. The first two, containing 53 weakly bound bimolecular associates and 13 molecular clusters taken from the crystal, were used in protocol refinement, and the third one containing other 20 bimolecular and three cluster systems served as a validation reference. In addition, to verify the performance of the proposed approach on an exact 2-RDM, calculations within the coupled cluster formalism were performed for part of the first set systems using the cc-pVTZ basis set. The process of optimization of the proposed parametric model is considered, and the role of various energy contributions in the formation of non-covalent interactions is discussed with regard to the obtained trends. [ABSTRACT FROM AUTHOR]

The structure and dynamics of water on solid surfaces critically affect the chemistry of materials in ambient and aqueous environments. Here, we investigate the hydrogen bonding network of water adsorbed on the majority (101) surface of anatase TiO2, a widely used photocatalyst, using polarization- and azimuth-resolved infrared spectroscopy combined with neural network potential molecular dynamics simulations. Our results show that one monolayer of water saturates the undercoordinated titanium (Ti5c) sites, forming one-dimensional chains of molecule hydrogen bonded to surface undercoordinated bridging oxygen (O2c) atoms. As the coverage increases, water adsorption on O2c sites leads to significant restructuring of the water monolayer and the formation of a two-dimensional hydrogen bond network characterized by tightly bound pairs of water molecules on adjacent Ti5c and O2c sites. This structural motif likely persists at ambient conditions, influencing the reactions occurring there. The results reported here provide critical details of the structure of the water–anatase (101) interface that were previously hypothesized but unconfirmed experimentally. [ABSTRACT FROM AUTHOR]

Faraday polarization rotation control of 1529 nm wavelength has been obtained using the 5 S 1 / 2 → 5 P 3 / 2 → 4 D 5 / 2 transition of rubidium. The traditional off-resonant polarization rotation method commonly used in the Faraday effect for direct transitions of atoms is not applicable to the transitions between excited states. In this study, we proposed a near-resonant polarization rotation method using the Faraday effect between excited states of atoms, by which the polarization rotation of the rubidium 5 P 3 / 2 → 4 D 5 / 2 transition can be controlled from 0 ° to 90 ° with low distortion. The 780 nm pump light corresponding to the 5 S 1 / 2 → 5 P 3 / 2 transition can also be used simultaneously as a control source to manipulate the polarization state of the 1529 nm optical signal. [ABSTRACT FROM AUTHOR]

In this work, we have proposed a Continuous Constant Potential Model (CCPM) based on grand canonical density functional theory for describing the electrocatalytic thermodynamics on single atom electrocatalysts dispersed on graphene support. The linearly potential-dependent capacitance is introduced to account for the net charge variation of the electrode surface and to evaluate the free energetics. We have chosen the CO2 electro-reduction reaction on single-copper atom catalysts, dispersed by nitrogen-doped graphene [CuNX@Gra (X = 2, 4)], as an example to show how our model can predict the potential-dependent free energetics. We have demonstrated that the net charges of both catalyst models are quadratically correlated with the applied potentials and, thus, the quantum capacitance is linearly dependent on the applied potentials, which allows us to continuously quantify the potential effect on the free energetics during the carbon dioxide reduction reaction instead of confining it to a specific potential. On the CuN4@Gra model, it is suggested that CO2 adsorption, coupled with an electron transfer, is a potential determining step that is energetically unfavorable even under high overpotentials. Interestingly, the hydrogen adsorption on CuN4@Gra is extremely easy to occur at both the Cu and N sites, which probably results in the reconstruction of the CuN4@Gra catalyst, as reported by many experimental observations. On CuN2@Gra, the CO2RR is found to exhibit a higher activity at the adjacent C site, and the potential determining step is shifted to the *CO formation step at a wide potential range. In general, CCPM provides a simple method for studying the free energetics for the electrocatalytic reactions under constant potential. [ABSTRACT FROM AUTHOR]

The associability of artificial intelligence (AI) as a dual-use technology for nuclear material production (NMP) within the academic and practitioner communities remains widely neglected, and so a widening opportunity for AI to aid in illicit and covert non-peaceful applications exists. To address this emerging gap, this paper investigates the evolving and applicable uses of AI and finds broad evidence of its use to optimize performance, promote innovation, reduce costs, and enhance safety associated with the development and production of nuclear material. AI's use in this arena will, thereby, facilitate broader accessibility of peaceful uses of nuclear science and technology, while at the same time cause concerns that said improvements can aid the illicit development of nuclear weapons. As such, this paper advocates for a three-dimensional solution to manage the evolving dual-use concern of AI that involves advancing states-centric monitoring and regulation, promoting intellectual exchange between the nonproliferation sector and the AI industry, and encouraging AI industrial contributions. [ABSTRACT FROM AUTHOR]

The pseudopotentials and dispersion potentials are applied to a theoretical study of the hyperfine splitting frequencies of the ground-state paramagnetic hydrogen (H) and alkali-metal (Li, Na, K, Rb, and Cs) atoms in noble gases (He, Ne, Ar, Kr, and Xe). Using classical turning points for statistical averages, we find that numerical calculations based on second-order perturbation theory fit the measured frequency shifts well over a wide temperature range. The characteristic energy, pseudopotential height, and electric-dipole polarizability allow us to consistently determine the van der Waals radii and electron scattering lengths of noble-gas atoms. This study shows that the hyperfine splitting frequency of alkali-metal atoms is a good measure for investigating colliding partners. [ABSTRACT FROM AUTHOR]

The nature of an atom in a bonded structure—such as in molecules, in nanoparticles, or in solids, at surfaces or interfaces—depends on its local atomic environment. In atomic-scale modeling and simulation, identifying groups of atoms with equivalent environments is a frequent task, to gain an understanding of the material function, to interpret experimental results, or to simply restrict demanding first-principles calculations. However, while routine, this task can often be challenging for complex molecules or non-ideal materials with breaks in symmetries or long-range order. To automatize this task, we here present a general machine-learning framework to identify groups of (nearly) equivalent atoms. The initial classification rests on the representation of the local atomic environment through a high-dimensional smooth overlap of atomic positions (SOAP) vector. Recognizing that not least thermal vibrations may lead to deviations from ideal positions, we then achieve a fuzzy classification by mean-shift clustering within a low-dimensional embedded representation of the SOAP points as obtained through multidimensional scaling. The performance of this classification framework is demonstrated for simple aromatic molecules and crystalline Pd surface examples. [ABSTRACT FROM AUTHOR]

A new approach for generating Gaussian basis sets is reported and tested for atoms from H to Ne. The basis sets thus calculated, named SIGMA basis sets, range from DZ to QZ sizes and have the same composition per shell as Dunning basis sets but with different treatment of the contractions. The standard SIGMA basis sets and their augmented versions have proven to be very suitable for providing good results in atomic and molecular calculations. The performance of the new basis sets is analyzed in terms of total, correlation, and atomization energies, equilibrium distances, and vibrational frequencies in several molecules, and the results are compared at several computational levels with those obtained with the corresponding Dunning and other basis sets. [ABSTRACT FROM AUTHOR]

We consider a one-dimensional effective quantum electrodynamics (QED) model of the relativistic hydrogen-like atom using delta-potential interactions. We discuss the general exact theory and the Hartree–Fock approximation. The present one-dimensional effective QED model shares the essential physical feature of the three-dimensional theory: the nuclear charge polarizes the vacuum state (creation of electron–positron pairs), which results in a QED Lamb-type shift of the bound-state energy. Yet, this 1D effective QED model eliminates some of the most serious technical difficulties of the three-dimensional theory coming from renormalization. We show how to calculate the vacuum-polarization density at zeroth order in the two-particle interaction and the QED Lamb-type shift of the bound-state energy at first order in the two-particle interaction. The present work may be considered a step toward the development of a quantum-chemistry effective QED theory of atoms and molecules. [ABSTRACT FROM AUTHOR]

Carbon nanotubes can serve as one-dimensional nanoreactors for the in-tube synthesis of various nanostructures. Experimental observations have shown that chains, inner tubes, or nanoribbons can grow by the thermal decomposition of organic/organometallic molecules encapsulated in carbon nanotubes. The result of the process depends on the temperature, the diameter of the nanotube, and the type and amount of material introduced inside the tube. Nanoribbons are particularly promising materials for nanoelectronics. Motivated by recent experimental results observing the formation of carbon nanoribbons inside carbon nanotubes, molecular dynamics calculations were performed with the open source LAMMPS code to investigate the reactions between carbon atoms confined within a single-walled carbon nanotube. Our results show that the interatomic potentials behave differently in quasi-one-dimensional simulations of nanotube-confined space than in three-dimensional simulations. In particular, the Tersoff potential performs better than the widely used Reactive Force Field potential in describing the formation of carbon nanoribbons inside nanotubes. We also found a temperature window where the nanoribbons were formed with the fewest defects, i.e., with the largest flatness and the most hexagons, which is in agreement with the experimental temperature range. [ABSTRACT FROM AUTHOR]

The regulation of intramolecular vibrational energy redistribution (IVR) to influence energy flow within molecular scaffolds provides a way to steer fundamental processes of chemistry, such as chemical reactivity in proteins and design of molecular diodes. Using two-dimensional infrared (2D IR) spectroscopy, changes in the intensity of vibrational cross-peaks are often used to evaluate different energy transfer pathways present in small molecules. Previous 2D IR studies of para-azidobenzonitrile (PAB) demonstrated that several possible energy pathways from the N3 to the cyano-vibrational reporters were modulated by Fermi resonance, followed by energy relaxation into the solvent [Schmitz et al., J. Phys. Chem. A 123, 10571 (2019)]. In this work, the mechanisms of IVR were hindered via the introduction of a heavy atom, selenium, into the molecular scaffold. This effectively eliminated the energy transfer pathway and resulted in the dissipation of the energy into the bath and direct dipole–dipole coupling between the two vibrational reporters. Several structural variations of the aforementioned molecular scaffold were employed to assess how each interrupted the energy transfer pathways, and the evolution of 2D IR cross-peaks was measured to assess the changes in the energy flow. By eliminating the energy transfer pathways through isolation of specific vibrational transitions, through-space vibrational coupling between an azido (N3) and a selenocyanato (SeCN) probe is facilitated and observed for the first time. Thus, the rectification of this molecular circuitry is accomplished through the inhibition of energy flow using heavy atoms to suppress the anharmonic coupling and, instead, favor a vibrational coupling pathway. [ABSTRACT FROM AUTHOR]

Herein, the first-principles technique combined with a coherent potential approximation is utilized to investigate the magnetic anisotropy energy (MAE) of (Fe1-xCox)16(N,C)2 alloys. The alloys contain an ordered α″-phase and partially ordered α′-phase, where α-Fe is distorted to form a bct structure with c/a = 1.1 owing to the random intrusion of (N, C) atoms. For the α″-phase, the MAE is sensitive to Co substitution sites and reaches approximately 3 MJ/m3 at an x of 0.25 when Co atoms occupy the 4d sites. For the α′-phase, the lattice distortion c/a = 1.1 alone cannot provide a positive MAE when x = 0 (FeN0.125 and FeC0.125). Evidently, both Co substitution and further increasing the c/a over 1.1 must be performed to realize a positive MAE in the α'-phase. Finally, the importance of local distortion of the (Fe, Co) lattice in further increasing the c/a beyond 1.1 is discussed as it increases the effective c/a even when the global c/a is fixed at 1.1. [ABSTRACT FROM AUTHOR]

Koopmans spectral functionals are a class of orbital-density-dependent functionals designed to accurately predict spectroscopic properties. They do so markedly better than their Kohn–Sham density-functional theory counterparts, as demonstrated in earlier works on benchmarks of molecules and bulk systems. This work is a complementary study where—instead of comparing against real, many-electron systems—we test Koopmans spectral functionals on Hooke's atom, a toy two-electron system that has analytical solutions for particular strengths of its harmonic confining potential. As these calculations clearly illustrate, Koopmans spectral functionals do an excellent job of describing Hooke's atom across a range of confining potential strengths. This work also provides broader insights into the features and capabilities of Koopmans spectral functionals more generally. [ABSTRACT FROM AUTHOR]

The origin of β-relaxation in metallic glasses is still not fully understood, and the guidance of slow atoms for caged dynamics and β-relaxation is rarely mentioned. Using molecular dynamics simulations, we reveal the bridging role of slow atoms on unusual caged dynamics and β-relaxation. In the stage of unusual caged dynamics, slow atoms are bounded by neighboring atoms. It is difficult for the slow atoms to break the cage, producing more high-frequency vibration, which causes more atoms to jump out of the cage randomly in the next stage. Precisely, the movement of the slow atoms changes from individual atoms vibrating inside the cage and gradually breaking out of the cage into a string-like pattern. The string-like collective atomic jumps cause decay of the cages, inducing β-relaxation. This situation generally exists in binary systems with the large atomic mass difference. This work offers valuable insights for understanding the role of slow atoms in unusual caged dynamics and β-relaxation, complementing studies on the origin of β-relaxation in metallic glasses and their glass-forming liquids. [ABSTRACT FROM AUTHOR]

Ecological traits are functional characteristics measurable at the species level and provide valuable insights into how organisms respond to environmental constraints. Here, we investigated how diatom trait‐groups and individual species respond to environmental variables, and identified indicator species that are particularly sensitive to environmental variation. Diatoms were sampled at 129 sites in the subarctic streams of Norwegian islands and mainland, and were categorised into three trait groups: high‐profile species that live in an erect position, low‐profile species living in low position along the surface, and motile diatoms. Data were analysed using a recently developed method known as Hierarchical Modelling of Species Communities, which is a flexible framework for joint species distribution modelling. We found that diatom trait‐groups responded relatively weakly to measured environmental variables but showed positive or negative relationships with major ion levels of the water (e.g., conductivity, calcium [Ca2+], sodium [Na+] or chlorine [Cl−]). Variance partitioning showed a similar, important contribution of the chemical variables for all of the trait groups, while the contributions of physical variables and especially random (spatial) factors were notably lower for all trait groups. Our findings also highlighted considerable among‐species variation in their relation to environmental variables within the trait groups. Notably, we identified a high number of indicator species within each trait group that were explained by specific environmental factors, mostly chemical variables (conductivity, pH, total nitrogen and phosphorus, Ca2+, Na+, Cl−). Our study suggests that certain diatom species can be considered as useful environmental indicators but the variability in species preferences within the trait group may in some circumstances hamper the use of ecological traits in environmental assessments. Thus, we suggest using species‐level ecology combined with trait information to better track environmental change when using diatoms as indicators. [ABSTRACT FROM AUTHOR]